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High performance, low cost, and relative 


ease of storage and handling make liquid oxygen 
not only the best oxidizer available but guarantee 
it a prominent place among high performance 
oxidizers in the future. 


Liquid oxygen has certain advantages when 
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1. A potential supply unlimited 


2. Production possible at source 


3. Inexpensive 


4, Non-toxic, non-corrosive, smoke-free 
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—_-§. Easily ignited and smooth burning 


6, Few limitations on component ma- 
terials of construction 


7. appreciably by ambient 
emperatures 
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An analysis has been made of rocket engines which con- 
vert heat into mechanical energy by means of an iso- 
The basic characteristics of such en- 


thermal expansion. 
zines have been derived and compared with conventional 
rocket engines utilizing adiabatic expansions to obtain 
thrust. It has been found that, in general, an isothermal 
expansion is a less desirable work cycle than an adiabatic 


expansion for the conversion of heat into mechanical 


energy in jet propulsion devices. 


Nomenclature 


= area of duct 
8 = ratio of limiting isothermal expansion temperature to the 
limiting adiabatic chamber temper: ature 

¢, = specific heat at constant pressure 

€ = area ratio of nozzle 

= enthalpy per unit mass 


AH = heat release associated with the addition of a unit mass of 
material to the gas stream 

= specific impulse 

Tmax= Maximum specific impulse attainable from a given pro- 
pellant combination 

(1:/I.¢= ratio of the isothermal to the adiabatic specific impulse 
for the condition of excess heat available 


« = ratio of the limiting adiabatic temperature to the stagna- 
tion temperature of the isothermal flow en 

l= specifie heat ratio 

= Mach number at any station 

VW, = Mach number at the nozzle throat in an isothermal expan- 
sion 

M, = a particular Mach number defined in the text 

p = static pressure 

Q = heat release per unit mass 

! = mixture ratio, defined in text 

= mass density 

R = ideal gas constant for unit mass of fluid 

7 = absolute static temperature 

7 = isentropic stagnation temperature 

7, = adiabatie static flame temperature based upon | a 
at a fixed reference temperature. a 

7, = isothermal static expansion temperature 

v= velocity of gas 

w = mass rate of flow of gas i ~~ ‘ 


Presented at the Annual Convention of the AMERICAN RocKET 
Socrery, Atlantie City, N. J., November 30, 1951. 
! Research Engineer, Flight Research Laboratory. 
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Superscript 
* = the value of a quantity at the nozzle throat 
Subscripts 
= 
t = the value of the quantity for isothermal expansion 
s = the value of the quantity for isentropic expansion 
1 = the value of a quantity at unity Mach number 
c = the value of a quantity in the combustion chamber 
e = the value of a quantity at the nozzle exit 


T HAS recently been suggested by Zwicky (1)? that a | 


rocket engine operating on an isothermal expansion 


might offer attractive advantages. This kind of 
thermodynamic process is clearly desirable from the 
standpoint of minimizing the temperature which ma- 
terials of construction must withstand. However, it~ 
is necessary also to consider the characteristics of such 
a process with respect to propellant economy and to the 
geometry of the associated engine. Seifert and Altman 
(2) have analyzed one special case of the general iso- 
thermal expansion problem and compared certain of its | 
features with those of a conventional isentropic expan- 
sion. The purpose of this paper is to develop the char- 
acteristics of isothermal rockets in a more general way; 
it is also intended to compare the isothermal and adia- 
batic cycles on a more direct and meaningful basis than 
has heretofore been done. 

An isothermal expansion in a jet engine 
ceived to occur in one of two ways which, for the sake 
of convenience, we shall refer to as Case 1 and Case 2 


may be con- 


in subsequent discussion. 

In Case 1, all of the working fluid in the engine is 
introduced into the combustion chamber, and the rate 
of injection, timing of the propellant flows, and geom- 
etry of the chamber and nozzle are such that partial 
combustion occurs in the chamber and raises the gas 
to a predetermined temperature level just before it 
enters the nozzle. The remainder of the combustion 
process is completed in the nozzle as the gas expands to 
its final pressure. The expansion in the nozzle is 


imagined to be so scheduled that in any section of the 


2 Numbers in parentheses refer to the References on page 287. 
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nozzle the release of heat by combustion is just suffi- 
cient to balance the drop in temperature due to expan- 
sion, thus maintaining the gas at a constant static 
temperature. In terms of a gas-dynamic analogy, Case 
1 corresponds to the continuous addition of heat to a 
gas flowing in a variable area duct. If all of the heat 
to the working fluid is supplied externally through heat 
exchangers, as for example in a jet engine activated by a 
nuclear reactor, the gas-dynamic analogy corresponds 
exactly to the physical model of the engine. 

In Case 2, only a portion of the working fluid is 
introduced into the combustion chamber; the re- 
mainder is injected into the expanding gas in the nozzle 
at just such a rate that its combustion maintains the gas 
at a constant temperature. In this instance, the total 
mass of working fluid in the nozzle cannot be assumed 
constant without the introduction of a serious error. 
Hence, the gas-dynamic equivalent for the Case 2 
rocket engine must be taken to correspond to the con- 
tinuous addition of heat and mass to a gas flowing in a 
variable area duct. 

"Development of Flow Equations > 

The | problems of predicting the combined effect of a 
number of simultaneous influences (such as area change, 
heat addition, and mass addition) on the properties of a 
moving gas stream are very complex and will therefore 
be approached from the one-dimensional point of view. 
This treatment is justified by the satisfactory results 
achieved by the application of one-dimensional gas- 
dynamics to the problem of flow in conventional adia- 
batic nozzles, and should serve to supply a valid com- 
parison between the isothermal and isentropic expan- 
sion cycles. In addition to the assumption of one- 
dimensional flow, it will also be assumed for this analy- 
sis that: (a) There are no discontinuities in the flow, 
and changes in stream properties are continuous. (6) 
The gas obeys the ideal gas law; also, the specific heat 
and molecular weight are the same for products and 
reactants and remain constant during the flow. (c) 
Fluid friction is negligible. (This assumption is more 
optimistic in the present problem than it is for the 
case of an adiabatic expansion. Considerable fluid 
turbulence is created by the process of heat and/or 
mass addition and this energy dissipation would be ex- 
pected to be more significant than that due to wall fric- 
tion alone.) 

Consider a section of infinitesimal length of a nozzle 
in which the processes with which we are concerned are 
occurring. The mass injected at this section of the 
nozzle is assumed to enter as a gas with a negligible 
velocity component parallel to the moving gas stream. 
The external heat addition arises from the chemical 
reaction of the entering gases and is proportional to the 
amount of these gases. The conditions characterizing 
the flow through this section are as follows: 

(a) Continuity 


dw _de dA, dv 


(b) Conservation of Momentum 
dp , kM*dv? | ,,,,dw 
w 
(c) Conservation of Energy 
dQ = dH = + 1/2 dv®............ [3] 
In these equations, dH denotes the net enthalpy change 


associated with combustion. The Mach number is re- 
lated to the other variables by the equation: 


dM? dv? dT 
M2 2? T 


It will be convenient, for our purposes, to make use of 
the isentropic stagnation temperature, which is defined 
as follows: 


The energy equation, in terms of this quantity, becomes 


The above equations can be cast directly into usable 
form by application of the method of influence coeff- 
cients developed by Shapiro and Hawthorne (3). The 
working equations for the flow then become: 

For Case 1 (heat addition only) 


dTy 
(7] 
ant 2 2 
dT, 
_ kM? ™ 


For Case 2 (heat addition + mass addition) 


2 2 2 
~ A — M? 

2 
aT _ 9 (k= jdt 
_ — 1) [1 + kM?) M? dw 11] 
k-1 
2 
p 1—MA 1 — M? To 


Equations [7]-[12] completely define the character- 
istics of the flow for both Case 1 and Case 2. However, 
in order to obtain unique solutions for Case 2, the 
mode of mass addition must be specified in terms of the 
other variables in the flow. 

For this purpose, consider the following model. All 
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of one component (either the oxidizer or the fuel) is 
introduced into the combustion chamber at the iso- 
thermal expansion temperature 7;. The deficient com- 
ponent is added along the duct, as a gas at temperature 
T,, at the proper rate to maintain an isothermal expan- 
sion. The addition of mass is terminated when the 
stoichiometric mixture ratio is attained for the propel- 
lant combination under consideration. Since the pro- 
pellants are assumed to have been preheated, the heat 
of combustion of the propellant combination is cor- 
rected by subtracting the sensible heat and latent heat 
of vaporization of the propellants. Then AH, the ef- 
fective heat of combustion so obtained, will be the 
energy release associated with the reaction of a unit 
mass of the added component with the gas stream. 
\ simplified combustion reaction is assumed for which 
the heat of reaction is strictly proportional to the quan- 
tity of the second component introduced. Considering 
the energy flux into and out of an element of volume, it 
follows that 


[AH — C,(T. — = weAT......... [13] 
where 7p is the total temperature of the stream entering 
the volume element. Furthermore, it follows that 

dw dT» 


[14] 
from the definition of AH, it is evident that 
[15] 


where 7’, is the adiabatic flame temperature for a 
stoichiometric mixture. Hence 

dw dT) 

To 


It will be convenient at this point to introduce another 
parameter, M,, defined as follows: 


M, would be the Mach number of the flow in the duct if 
the isothermal expansion were carried out to the point 
where the total temperature of the stream is equal to 
the adiabatic flame temperature. In terms of this 
parameter 


where M is the Mach number of the flow when the total 

temperature is 7’. It is now possible to solve the flow 

equations for both cases to arrive at integral relations 
between the significant variables. 

For Case 1, operating upon Equations [7], [8], and 

[9] we find that 


exp E (Me — »| [19] 


The subscript 1 appearing in the above equations de- 
notes the value of the quantity so designated at the 
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station for which the Mach number is unity. a 


= 


Similarly for Case 2, the following relations can san 
derived: 


= exp [ - 5 (M? — 1) 


2+ (k — 1)M? 


M,;? — M2] kMP k 
Fos | exp | 5 (M? [24] 


where, as before, the subscript denotes the value of the 
quantity appearing in each of the equations at a Mach 
number of unity. ; 

Since the Mach number at the throat need not neces-_ 
sarily be unity for an isothermal expansion, it will be — 
instructive to determine the conditions at the throat for 
both cases of isothermal expansion. For Case 1, the 
value of the Mach number that yields a minimum value 
for the duct area is found from Equation [19] to be 


7) 


Similarly, from Equation [22] the value of the Mach — : 
number at the throat, M,, is given by 


_ —(kMy? +3) + +3)? + [26] 


M? ok 


We are now in a position to express the flow relations © 
in terms of the familiar rocket design parameters ¢ and — 
p-/Pe. Cast in these terms they become: 

For Case 1 


= in [27] 
) 
T.=T: In [29] 


For Case 2 


27 -k My? 
DelDe = [1 Gr) ] exp [ - [30] 
M; 2 


In the temperature expression above for Case 1, 7’; is the _ 
isentropic stagnation temperature of the flow at the 
nozzle exit and is identical with the adiabatic flame 
temperature. 7’; has been defined previously as the iso- 
thermal expansion temperature. For Case 2, the ratio — 
T,/T, is independent of the pressure ratio and may be © j 
assigned upon the basis of other considerations. Its — 
value will determine the Mach number attainable for a — 
given pressure ratio and the geometric expansion ratio | 
necessary to achieve this pressure ratio. It will be — 
seen, upon inspection of the flow relations for Case 2, — 
that for any temperature ratio selected, the free stream — 
Mach number approaches the value M, as the pressure — 
ratio p,/p. and the expansion ratio approach infinity. 
In other words, these relations state that for the case — 
of an isothermal expansion involving heat and mass — 
addition, the total temperature of the stream can attain 
the adiabatic flame temperature only at the expense of 


an infinite expansion and pressure ratio. 
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FIG. 1 RATIO OF ISOTHERMAL EXPANSION TEMPERATURE TO 
ADIABATIC FLAME TEMPERATURE AGAINST PRESSURE RATIO FOR 
CORRECT EXPANSION. 


The principal attraction of an isothermal rocket lies 
in the possibility of operating it at temperatures lower 
than those prevailing in conventional rockets. In Fig. 
1, the ratio of the isothermal expansion temperature to 
the adiabatic flame temperature is plotted as a function 
of the pressure ratio for conditions of correct expansion 
for the isothermal Case 1 (heat addition alone). Values 
for this plot were calculated from Equation [29]. It is 
interesting to note that in the range of pressure ratios 
normally used in rockets (20-50), the static temperature 
of the isothermal expansion is only about 60 per cent of 
the adiabatic flame temperature. It is well to observe 
also that the effect of pressure ratio diminishes with in- 
creasing pressure ratio, and that further substantial 
reductions in temperature can be achieved only at the 
expense of a vastly increased pressure ratio. Since for 
Case 2 the temperature ratio is independent of the 
pressure ratio, the comparison of the characteristics of 
the two types of isothermal expansions will be made at 
values of 7/7, equal to 0.6 and 0.1 for the Case 2 
expansion. The former value of 7,/7, will yield a 
comparison of the two processes at approximately the 
same expansion temperature; the latter value will 


serve to demonstrate the effect of temperature ratio on 
an isothermal expansion with mass addition. 

The geometries of isothermal — are compared 
with those of isentropic nozzles in Fig. 2, which depicts 
the relations between area expansion ratio and pressure 


R/% 
FIG. 2 AREA EXPANSION RATIO AGAINST PRESSURE RATIO. 


isentropic expansions. 


ratio. Isentropie expansion and both types of iso- 
thermal expansion are shown, at a specific heat ratio 
of 1.3. The exit mass flows for the adiabatic expansion 
and the Case 1 isothermal expansion are equal. The 
exit mass flow for the Case 2 isothermal expansion is a 
function of the pressure ratio and is smaller than the 
exit mass flow for isentropic expansion at any finite 
pressure ratio under the conditions of the present analy- 
sis, which postulates equal mass flows at the time that 
all available chemical energy in the propellant has been 
released to the gas stream. The curves indicate that 
over the entire range of pressure ratios investigated, the 
isentropic nozzle requires a considerably smaller area 
ratio to expand to a given pressure ratio than either of 
the isothermal nozzles. The area ratio for the Case 2 
isothermal expansion (with mass addition) is greater 
than that for the Case 1 process, increasing as the tem- 
perature ratio 7,/T, is increased. 

To define completely the geometry of the nozzle, it is 
necessary to know, in addition to its expansion ratio, 
the throat area required for the passage of a given mass 
flow at any operating pressure ratio. We shall now de- 
rive expressions relating the throat area of an adiabatic 
nozzle to the throat areas of the two types of isothermal 
nozzles for the same total mass flow at the completion 
of chemical reaction. 

Consider first the case of heat addition only. The 
mass of gas passing the throat is the same for both the 
adiabatic and isothermal flows; hence 


where the subscripts t refer to the isothermal expansion, 
the subscripts s refer to the adiabatic expansion, and 
the asterisk denotes the values of the quantities at the 
throat of each nozzle. After insertion of the proper 
values for the density and velocity of the gas stream 
at the throat for both expansions in Equation [32], the 
following expression is obtained after some algebraic 
manipulation: 


k+1 


For the Case 2 isothermal process, with mass addi- 
tion, the total mass flow through the throat is not equal 
to the common value prevailing for the Case 1 and 
Recasting Equation [18] gives 


Lin (p. . [33] 


dw 2MdM 


where w, is the mass flow of propellant entering the 
nozzle from the combustion chamber, w* is the mass 


flow of propellant at the throat, M, is the Mach num- 
ber of the flow at the throat, and M, is the parameter 
defined previously. If r is the weight ratio of the defi- 


cient component to the other component in a stoi- 
chiometric mixture, the ratio of the mass flows at the 
_ throat will be given by: 
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HEAT ADDITION 
R=1.3 
= 0.6 
O4F € MASS ADDITION 
O2 - = 0.1 
fe) 50 00 
R./Pe 
11G. 3. RATIO OF THROAT AREA FOR ISOTHERMAL EXPANSION TO 
'HROAT AREA FOR ADIABATIC EXPANSION AGAINST PRESSURE 
RATIO. 
w* M;* 


(1+r(M? — M2) 


where w is the mass rate of flow through the isentropic 
nozzle. Treating this case in a manner similar to that 
used in the previous instance, the result is 


A,* k+1 2(k—1) 


The throat area ratios for both types of isothermal ex- 
pansion are shown in Fig. 3. For the case of heat addi- 
tion alone, the throat area is approximately 10 per cent 
less than that for the adiabatic expansion over the en- 
tire pressure range plotted and decreases gradually with 
increasing pressure ratio. For the case of both heat and 
mass addition, the throat area ratio is independent of 
the pressure ratio and is a function of the temperature 
ratio T,/T, alone. The curves in Fig. 3 have been 
calculated for a mixture ratio r equal to unity. The 
throat area for the case of heat and mass addition is 
considerably smaller than it is for an adiabatic expan- 
sion and decreases with decreasing temperature ratio. 

In order to achieve a balanced evaluation of the char- 


2\ —kMys? 
exp [-3 . [37] 
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FIG. 4 RATIO OF ISOTHERMAL TO ADIABATIC SPECIFIC IMPULSE 
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acteristics of an isothermal expansion, the possible pen- 
alty in performance incurred as a result of the lowered 
operating temperature must be investigated. This is 
best shown by the relation between the specific impulse 
of an isentropic and an isothermal expansion as a 
function of the operating pressure ratio. The velocity 
of the gas stream at the nozzle exit may be determined 
from the following equation for an isothermal expan- 
sion: = 


vr? = 

in which »,, is the velocity and M, the Mach number of 
the flow at the nozzle exit. The velocity of a gas at the 
end of an adiabatic expansion may, on the other hand, 


be represented by 
k-—1 
(39 
Pe 


2kR 
T 
The ratio of the isothermal to the adiabatic specific 
impulse will be equal to the ratio of the terminal gas 
velocities for the two types of expansions and will be 
given for both Case 1 and Case 2 by 


2/7, (40) 


1 (2: ) k 
Pe 


For the Case 1 isothermal expansion, both the tem- 
perature ratio and the Mach number at the nozzle exit 
are uniquely determined by the operating pressure 
ratio and may be replaced by their equivalent expres- 
sions in terms of pressure ratio with the following 
result: 


k-1 
k In (2 ) 
k | pe) * \ 

The comparison between the performance of an 
isothermal and an isentropic expansion operating at 
the same pressure ratio is presented in Fig. 4 as a func- 
tion of pressure ratio. For the heat and mass addition 
process (Case 2), the specific impulse ratio is a function | 
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of the expansion temperature as well as the pressure 
ratio. This type of isothermal expansion is seen to be 
much inferior to the Case 1 expansion from the stand- 
point of propellant economy; moreover, it becomes pro- 
gressively worse with decreasing temperature ratio over 
the range of values considered. Even the Case 1 proc- 
ess, however, is slightly inferior in performance to the 
isentropic expansion, its relative performance also 
decreasing with increasing pressure ratio for this range 
of values. Actually the specific impulse ratio for the 
Case 1 expansion eventually goes through a minimum 
and begins to increase at higher pressure ratios, ap- 
proaching unity asymptotically. The curves for the 
Case 2 expansion also exhibit a different behavior at 
pressure ratios much higher than those plotted here; 
the temperature ratio curves cross over, and the specific 
impulse ratio increases with decreasing temperature 
ratio at a constant pressure ratio. However, the high 
pressure ratios at which these effects set in are beyond 
the range of practical interest. 

It is instructive to compare the thermodynamic effi- 
ciencies of the energy conversion process in the various 
kinds of expansion. One measure of this efficiency is 
the ratio of the specific impulse attained at any pressure 
ratio to the maximum specific impulse obtainable (at an 
infinite pressure ratio). The expressions for this pa- 
rameter as a function of pressure ratio are readily de- 
duced from Equations [38] and [39] for the gas velocity 


and are as follows: 
(a) Isentropie expansion 28 
4 

I/lmax = (®) [42] 

(b) Isothermal expansion, Case 1 
| De 1/2 re 

I/Imax k 1 Pel [43] 

In 


(c) Isothermal expansion, Case 2 


--- The plot of the thermodynamic efficiency against 


pressure ratio for all the expansions, shown in Fig. 5, 
reveals that the isentropic expansion is the most efficient 
in its utilization of pressure for the conversion of heat 
into mechanical energy. The isothermal expansion 
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involving heat and mass addition is significantly poorer 
in its pressure utilization than either of the other two 
processes and is adversely affected by decreasing tem- 
perature ratio. 

It is possible to conceive of at least one situation in 
which the performance of an isothermal jet engine can 
surpass that of an adiabatic engine, namely the situa- 
tion in which more heat is available to the working 
fluid in the engine than can be utilized in an adiabatic 
process due to structural limitations of the combustion 
chamber and nozzle materials. Such a condition might 
arise in the use of high impulse rocket propellants, 
where adiabatic flame temperatures in the order of 5000 
K are possible, or in the external addition of heat from 
a nuclear reactor which can supply a relatively limitless 
amount of energy. We shall complete our comparison 
of the isothermal and adiabatic expansion cycles with 
the analysis of this case. This will be carried out only 
for the Case 1 isothermal expansion, since the perform- 
ance of the Case 2 expansion is too poor to warrant 
further consideration. 

Let us now define the following quantities: 


T,, the limiting adiabatic chamber temperature 
(static and total are equal). 

8, the ratio of the limiting isothermal expansion 
static temperature to the limiting adiabatic cham- 
ber temperature. 

T.x, the stagnation temperature of the isothermal 
flow at the nozzle exit for a static temperature of 
BT. 

o, the ratio of the limiting adiabatic chamber tem- 
perature to the stagnation temperature of the iso- 
thermal flow. 

The value of 8 can never exceed unity and will be 
fixed by the maximum permissible rate of heat transfer 
in the nozzle, which in turn will be proportional to the 
stream total temperature, all other things being equal. 
Hence the limitations of the nozzle material will set 
a maximum tolerable stream total temperature. Pre- 
sumably this temperature will be lower than the adia- 
batic chamber temperature due to the greater value of 
the heat transfer coefficient in the nozzle as compared 
with the combustion chamber. 

The parameter o is in effect a measure of the relative 
amounts of heat utilized in the adiabatic and isothermal 
expansions, respectively, expressed as the ratio of the 
two quantities. It is related to the parameter 6 and 
the operating pressure ratio in the following way: 


The relation between the stagnation temperature of the 
isothermal flow and the adiabatic limiting temperature 
is given by: 


From the preceding two equations the exhaust velocity 
of the isothermal flow is found to be 


= [2c,T(1 — [47] 
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while for the adiabatic flow the exhaust velocity is 


k—1) 


The ratio of the isothermal to the adiabatic specific im- 
pulse for the case of excess heat available is then 


In Fig. 6 families of curves have been drawn which 
depict the variation of the specific impulse ratio with 
pressure ratio for a series of values of the parameter 8. 
Lines of constant o have been cross-plotted on the same 
co-ordinate system to form a network which completely 
defines the performance of an isothermal rocket relative 
to an adiabatic rocket when excess heat is available. 

It is evident that the performance increase obtained 
with an isothermal expansion, in the case of an un- 
limited energy source, is a function of the pressure ratio 
alone and may be read directly from the plots for the 
proper values of 8. In the case where the total energy 
available is limited, the curves, entered at appropriate 
values of o and 8, yield the pressure ratio necessary to 
utilize this amount of energy in an isothermal expansion 
and the associated increase in specific impulse. To 
illustrate, suppose that a certain propellant combina- 
tion possesses an adiabatic flame temperature of 4400 K 
in the pressure range of interest and that the maximum 
chamber gas temperature the proposed rocket can 
tolerate is 3500 K. Let us further suppose that the 
maximum temperature that can be tolerated in the 
nozzle gas is 0.8 times this value, and that the specific 
heat ratio for the gases concerned is 1.3. In this ex- 
ample, o and 8 are both equal to 0.8. Entering the 
network at these values, the pressure ratio required to 
utilize all the available chemical energy is found to be 
12, and the increase in specific impulse of the 
isothermal expansion over the adiabatic expansion is 
approximately 2'/. per cent. If a higher pressure ratio 
is available for the expansion, the isothermal expansion 
temperature may be reduced, but the gain in specific 
impulse will not change. A,” 


Conclusions 


The preceding analysis has demonstrated that, in 
general, an isothermal expansion is a less desirable proc- 
ess than an adiabatic expansion for the conversion. of 
chemical energy into mechanical energy in a jet propul- 
sion device. The thermodynamic efficiency and per- 
formance of an isothermal expansion cycle are lower than 
those of an adiabatic cycle, and larger area ratios are re- 
quired in the nozzle for the same pressure expansion 
ratio. These disadvantages, coupled with the formida- 
ble practical problems inherent in the attainment of an 
isothermal expansion, ordinarily outweigh any benefit 
of reduced operating temperatures. 
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The influence of local variations in mixture ratio in a 
rocket-motor combustion chamber on the motor perform- 
ance is analyzed with regard to arbitrary modes of such 
variation.‘ The results of this analysis show that the ef- 
fect in question can be evaluated by referring to a curve of 
theoretical specific impulse vs. weight fraction. In partic- 
ular, if such a curve has a negative second derivative 
everywhere, local variations in mixture ratio can only 
decrease performance. If, on the other hand, there is a 
region with positive second derivative, then local variations 
in mixture ratio may improve performance; and if these 
variations could be controlled, the maximum performance 
so obtainable would be on the double tangent near the over- 
all mixture-ratio point. 


= exhaust velocity 
c* = characteristic velocity 


Cy = specific heat at constant pressure 

f = cross-sectional area in nozzle 
F = thrust 

I = specific impulse a 17 
k = constant (cf. Equation [17]) } 7 
m = mass flow rate 
M = molecular weight of gases 

Pp = pressure 

r = mixture ratio (wt. oxidizer/wt. fuel) 


R = gas constant 
S = entropy 


7 = temperature 
Xi = weight fraction of total propellant in 7th een nt 
y(p) = po distribution which maximizes specific impulse 
n = thermodynamic efficiency age 
= Lagrange multiplier (cf. Equation [15]) 
p = weight fraction of oxidizer in propellant - 

g proj 
Subscripts 
c = chamber 
e = exit 
i = combustion element 
t = throat 
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Nomenclature 


(a) Consideration of Thermodynamic Limitations 


HERMODYNAMICALLY there does not seem 

to be any objection to the attainment of an in- 
creased or decreased performance if there is nonuniform 
mixing of the fuel and oxidizer streams. Even if the 
total enthalpy of the combustion gas, relative to the 
elements at some base temperature, is a constant (in an 
adiabatic system) regardless of mixing, the distribution 
of energy among its various forms can be considerably 
changed, and hence the entropy and thermodynamic 
efficiency of the system altered. Depending upon which 


I Thermodynamic Treatment 


_ additional energy states are available, the entropy can 


be shifted up or down by mixing variations, and hence 
_ the work potential of the combustion gas is increased or 
decreased, respectively. Two simple illustrations will 
_ suffice to demonstrate this effect. 
EXAMPLE 1: Consider a rocket motor A operating 
with helium gas heated to 2000 K, and an adjacent 
motor B identical in all respects except that half of the 
gas is at 1500 K and the other half at 2500 K. For the 
sake of argument, the difficulties in maintaining the 
sharp temperature boundary in the gas in motor B may 
be disregarded since this boundary has no relation to 
the principle concerned. Now the specific impulse J 
and the entropy S are given by the following equations 
(for a nondissociating gas with constant specific heat): 


| 


S = + C, In 300 ~ [2] 


7; = chamber temperature 

Pc = chamber pressure 

Pe = exit pressure 

Cy = specific heat at constant pressure 7 
n = thermodynamic efficiency 


Syoo® = standard entropy at 300 K 


As is shown in Section I (b), the J of motor B is deter- 


mined by taking the weighted average of the J values 
corresponding to the two temperatures. Thus 

1 

oy n(V1500 + 1/2500) .......... [3] 


¥ is now evident that J(A) > I(B) since 1/2000 > 
'/5(4/ 1500+ +/ 2500) and 7 is the same for both motors. 
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Similarly, it is found that oy | 


C,(In 2000 — In 1500 — In 2500) 
2000 
= Co In x 2500)", > ° 


Sa — Sz = 


Thus, for the case of a monatomic gas with constant 
specific heat in which only the translational energy 
states are activated, a nonuniform energy distribution 
will cause a lowering of the entropy and hence of per- 
formance. 

EXAMPLE 2: An examination will now be made of a 
svstem in which other energy states besides transla- 
tional are available. Consider a rocket motor A con- 
tuining hydrogen at 20.4 atm and 6000 K. At this 
temperature, equilibrium-constant data on H». show 
that it is about 92.5 per cent dissociated into atoms. 
(‘onsider also a rocket motor B identical with A in all 
respects except that half of the gas has 12 kcal more heat 
per gram and the other half has 12 keal less heat per 
gram. Equilibrium calculations show that the hotter 
half will attain a chamber temperature of 8000 K (H, 
is 100 per cent dissociated), an increase of 2000 K, 
whereas the cooler half will be 5450 K (Hz is 75.5 per 
cent dissociated), showing a drop of only 550 K. This 
effect, which may be contrasted with that of the motor 
containing helium, results from the fact that a lowering 
of the temperature causes H atoms to recombine, thus 
liberating heat. A performance calculation carried out 
for the three temperatures 5450, 6000, and 8000 K 
vields the IJ values (frozen expansion assumed for 
simplicity) of 1238, 1325, and 1550 sec, respectively. 
Thus the net J observed for rocket B would be !/, 
(1238 + 1550) = 1394 sec, which represents about a 5 
per cent increase in J although the total heat content 
per gram of hydrogen in the chamber is the same in both 
cases. If equilibrium were assumed, the magnitude of 
the increase would have been diminished, but the gen- 
eral conclusion would be the same. The entropy differ- 
ence between motors A and B is not significant in this 
case because the thermodynamic efficiencies of the two 


systems are different. 

These two examples show that, when temperature 
variations are imposed on an adiabatic system, the 
propellant performance may be shifted up or down de- 
pending upon the availability of the various internal and 
chemical energy states of the molecules. In a rocket 
system such temperature variations may be accom- 
plished through deviations in mixture ratio, but it would 
be necessary to know the equilibrium composition of 
ach mixture-ratio element before conclusions could be 
drawn concerning whether performance was increased 
or decreased. Since the character of the J-r plot de- 
pends on the thermodynamic distribution of energy in 
the combustion gas, conclusions concerning observed 
performance when inhomogeneous mixing occurs can be 
be drawn directly from such a curve. 


Performance of Rocket Containing Several Distinct 
Mixture Ratios in Combustion Gas 


(b) 


It will be demonstrated here that the parameters 
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E+ 


c*, c, and I for the whole motor can be evaluated as 
weighted averages of the corresponding local parame- 


ters® In carrying out this analysis, the following 
assumptions were made: (1) Each element at a given 
mixture ratio r of oxidizer to fuel remains intact 


throughout expansion. (2) Friction and thermal con- 
duction between filaments are negligible. (3) Each 
element reaches its minimum cross section at the throat 
and thus operates under the same constant pressure 
ratio. (4) Each element is in thermodynamic chemical 
equilibrium corresponding to its own composition. 


Let 


x; = weight fraction of total propellant in ith element 
fi = throat area occupied by ith element 

m; = mass flow rate of 7th element 

fi = total throat area 


The basic equations are 


m; 

c* Defi 


* is ob- 


The 
tained through the condition that 


relation between c* and the quantities c 


m 


fi = > fi = : = mic*; = 
i i 


Equations [6] and [7] can now be combined to yield 


. 


Similarly it can be shown that the exhaust velocities — 
and hence J should be averaged in the same way. If | 
F; = m,c; is the thrust contribution of the 7th element, 
then 


{9} 


F = = = n>, ric; = me 
i i i 


The last equality of Equation [9] shows that 
c= 


Since it has been shown that the parameters c*, c, and 
I are evaluated by weighted averages, the problem re-_ 
mains to determine what arbitrary weighting shall be — 
placed on the various mixture ratios in order to maxi-— 
This problem is treated in Section IT. 


[10] 


mize the J. 


Mathematical Analysis _ 


(a) Mathematical Formulation 


The problem of maximizing the specific impulse of a — 


5 “Influence of Non-Uniform Mixing of Fuel and Oxygen on 
Characteristic Data of a Rocket Motor,” by Carl Wagner, Fort 
Bliss, Texas, TR No. 22, Aug. 1947. In this report the influence 
of statistical deviations from uniform mixing is discussed, and 
the conclusion is reached that the effect on performance is least 
in the neighborhood of flat regions of the impetus-mixture ratio 


curve. This conclusion is implied as a special case of the present 


results. 
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rocket motor, as has been described earlier in this 
paper, can be interpreted mathematically as follows: 
a ot A: Let I(p) be a given function of p, single-valued, bounded, 
and nonnegative for p on the interval 0 to 1. Then for 
2 each value of p (for example, p = po) it is required to 


find a function y(p) which satisfies the conditions 


y(0) = Oand y(1) = 
y(p) increases monotonically with p.......... 


dy(p) = po where 0 < po 
" a and which maximizes the integral 


1:(p0) = I(p) dy(e) 


The function I(p) corresponds, of course, to the curve 
of specific impulse vs. weight fraction for the given pro- 
__ pellant system. The distribution function y(p) is that 
fraction of the total propellant which is burned at mix- 
ture ratios corresponding to values of the weight frac- 
tion less than or equal to p. Thus, condition [11] 
states that none of the propellant is burned at a value 
of p less than zero, whereas all of it is burned at values 
of p less than or equal to unity. Condition [12] asserts 
that the fraction burned at a given mixture ratio is not 
negative. Condition [13] is the statement of conserva- 
tion of mass. The integral to be maximized, J;(po), is 
the net specific impulse of the motor written as the 
7 . weighted average of local specific impulses (cf. Section 


Solution 


Impossibility of Continuous and Differentiable 


the problem involves maximizing an integral 
___-with regard to a variable function, the methods of the 

of variations are However, as 
will be seen next, the Euler equation of the calculus of 
_—- variations does not yield a useful solution. In fact, 
_ there is in general no continuous and differentiable 
solution of the type that satisfies the Euler equation. 
R= If there were such a solution, Euler’s equation would 
take the form 


dfo fe) 

wit] -2 +m) = 0... [15] 

where ) is a Lagrange i.e., a constant. 

_ Equation [15] reduces to 

[16] 


whence I’ is a constant. Thus a solution exists only if 
4 J is linear in p; that is, of the form 


In this case, Equation [14] becomes 


= + &) dy (0) ( 

=dpo + k = I(0) 

is, the maximized specific impulse function J, 

coincides with the given functionZ. Thus the only case 

for which the Euler equation yields a solution proves to 

be uninteresting. 
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(c) Geometric Derivation of Solution 


A direct method of attack yields better results. In 
fact, the general analytic solution to the problem as 
stated in Section II (a) can be given. The analytic 
solution, which involves considerable detail not relevant 
to the physical understanding of the problem, is given 
elsewhere. The following geometric interpretation is, 
perhaps, more enlightening. 

Consider the specific impulse function I(p) shown in 
the accompanying sketch. It is apparent that the value 
of I(p) at point C is exceeded by the weighted average 


of the values of I(p) at points A and B. This weighted 
average, in fact, lies on the straight line connecting A 
and B. On the other hand, the value of Z(p) at point B 
cannot be exceeded by any weighted average because 
no chord of the curve lies above B. 

Extrapolating these results gives the following cri- 
terion for increasing J(p) by means of weighted over- 
ages. If, for any value of p (for example, p = po), the 
curve I(p) has a chord passing above I(po) and has its 
ends on either side of po, the weighted average of the 
values of I(p) at the ends of the chord exceeds I(pp). 
If no such chord can be found, I(po) cannot be exceeded 
by weighted averages. It follows that the maximum 
weighted average is determined by the chord which is 
also tangent to I(p) at its ends or, more generally, which 
lies above I(p) over the whole range. 

Thus the maximum specific impulse is obtained as a 
weighted average of the specific impulse at only two 
points. Averaging over more than two points cannot 
give further improvement. 


III Application to Several Propellant Systems 


The results of Section II show that if the Z(p) curve is 
continuous everywhere, as is generally true for propel- 
lant systems investigated thus far, then the existence 
of a line touching at least two points of the Z(p) curve 
and lying above the curve defines a region in which 
I(p:) > I(po). This section is devoted to the examina- 
tion of two typical I(p) curves (calculated on basis of 
equilibrium expansion) in order to determine whether 
or not such tangent lines can be formed. 


(a) RFNA-Aniline 


Fig. 1 is a plot of the specific impulse I vs. the weight 
fraction p of RFNA. At the p values of 0.5 and 0.4, 


6 David Altman and Jack Lorell, Jet Propulsion Laboratory, 
Pasadena, Calif., PR No. 9-51, Jan. 15, 1951 (Restricted). 
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the chamber temperatures are 1345 and 1230 K, re- 
spectively, allowing for the formation of carbon and 
methane. Other hydrocarbons were not considered 
because of the complexity of the calculations, although 
it is not believed that their omission will appreciably 
alter the results. 

Since the J curve should go approximately to zero at 
p = 0, a reasonable extrapolation shows a point of in- 
flection near p = 0.3. This point permits the drawing 
of a tangent line between p = 0.3 and 0.7 maximizing I 
in that interval. Thus, in the mixture-ratio range of 
0.43 to 2.33, it is possible to observe performance above 
the theoretical as a result of inhomogeneous mixing. 


(b) NeHy-F.(D) 


Fig. 2 is a graph of the specific impulse of N2H, and 
F.(1) over a considerable range of p. The value of J 
at p = 0 has the value for N2H, as a monopropellant, 
corresponding to either of the reactions NsH, = Nz + 
2H, or 3NesH, = NH; + No. The entire portion of this 
curve from p = 0 to 0.83 shows no region where d?I/dp? 
> 0, and hence the maximum J coincides with the curve. 


IV Conclusions 


The results of this study have yielded the following 
practical conclusions: 


T 3 4NH,+ 


100 - 
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= Wt of / TOTAL wt 
FIG 2. SPECIFIC IMPULSE OF F,(/) vs. WEIGHT FRACTION 


or F, 


1 Uniform equilibrium combustion throughout the 
cross section of a rocket motor does not necessarily 
yield the highest performance attainable at that mix-_ 
ture ratio unless it corresponds to the over-all maximum > 
performance. 

2 Nonhomogeneous mixing of propellant streams 
can either raise or lower J depending on the shape of — 
the I-p curve. If d?J/dp* > 0 for some value of p, non- — 


uniform mixing will generally improve performance if 


the nonuniformities are not too large. 

3 The region of possible improvement of perform- 
ance by nonuniform mixing can be recognized quickly 
by examining an J-p curve. If the curve is such that a_ 
straight line can be drawn touching the curve at two > 
points but lying above the curve between these two- 
points, then an improvement in performance is possible 
for values of p between these two points. In fact, the 
new I-p curve follows this straight line. 

4 The facts presented in this paper point to a pos-_ 
sible explanation of a common observation that combus- 
tion efficiencies in a given motor and injector are gener-_ 
ally higher on the under-oxidized portion of the I-p | 
curve than near stoichiometric. Since all injectors — 
have some degree of nonuniformity, the shape of the — 
performance curve will influence the observed perform- m- 
ance at a given mixture ratio, 
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+ Let m»(t) be the mass rate of generation of hot gas by 


This paper shows that the combustion in the rocket 
motor can be stabilized against any value of time lag in 
combustion by a feedback servo link from a chamber 
pressure pickup, through an appropriately designed ampli- 
fier, to a control capacitance on the propellant feed line. 
The technique of stability analysis is based upon a com- 
bination of the Satche diagram and the Nyquist diagram. 
For simplicity of calculation, only low-frequency oscilla- 
tions in monopropellant rocket motors are considered. 


However, the concept of servo-stabilization and method 


of analysis are believed to be generally applicable to other 
cases. 


HE phenomenon of rough burning in liquid-pro- 

pellant rocket motor has been interpreted as the 
instability of the coupled system of propellant feed and 
combustion chamber by D. F. Gunder and D. R. Friant 
(1),2 M. Yachter (2), M. Summerfield (3), and L. Crocco 
(4). The essential feature of these theories is the time 
lag between the instant of injection of the propellant 
and the instant when the propellant is burned into hot 
gas. Crocco has further improved on this concept by 
considering the time lag as an integrated effect of con- 
secutive stages, each of which is controlled by the prevail- 
ing pressure in the combustion chamber. As a result 
of this new concept, Crocco showed the possibility of in- 
trinsic instability with constant injection rate not in- 
fluenced by the chamber pressure. 

The present paper will first give a slightly more gen- 
eral formulation of Crocco’s concept of time lag, allow- 
ing arbitrary pressure dependence of lag. Then the 
problem of intrinsic stability is discussed by applying a 
method suggested by M. Satche (5). This method is 
based upon a modification of the Nyquist diagram and is 
particularly useful for systems having time lag. For 
easy reference, this new diagram will be called the 
Satche diagram. The later sections of the paper will 
show the possibility of stabilizing the combustion by 
means of a feedback servo for all values of time lag. 
Such possibility of servo-stabilization was first men- 
tioned by W. Bollay in his admirable paper (6) on the 
application of servomechanisms to aeronautics. The 


present study definitely shows the power of this idea. 


a 


Time Lag in Combustion 


Consider, for simplicity, 
Then the mass rate of injec- 


combustion at time instant ¢. 

a monopropellant motor. 
Received February 22, 1952. 
1 Robert H. Goddard Professor of Jet Propulsion. 
2 Numbers in parentheses refer to the References on page 268. 
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tion at ¢ can be denoted by m;,(t). Let 7(¢) be the time 
lag for that parcel of propellant which is burned at the 
instant ¢. Then the mass burned during the interval! 
from ¢ to t + dt must be equal to the mass injected dur- 
ing the time from — rtot—7+d(t— 7). Thus 


mi(t)dt = mi(t — r)d(t — r)............ [1 


The mass of hot gas generated is either used to fill the 
combustion chamber by raising its pressure p(é), or is 
discharged through the rocket nozzle. If the frequency 
of the possible oscillations in the chamber is small, then 
the pressure in the chamber can be considered as uni- 
form, and as a first approximation (7) the rate of flow 
through the nozzle can be taken as proportional to the 
instantaneous chamber pressure p(/). Thus if m is the 
steady mass rate flow through the system, M, is the av- 
erage mass of hot gas in the chamber, and if the volume 
occupied by the unburned liquid propellant is neglected 


(7) dt +d (i, [2] 


where # is the steady state pressure in the combustion 
chamber. 

By following Crocco, the nondimensional variables 
for the chamber pressure and the rate of injection are 
defined as 


g and yw are then the fractional deviation of pressure and 
injection rate from the average. With Equation [3], 
m, can be eliminated from Equations [1] and [2], and 


dr 
4) {u(t —7r) + 1]..... [4] 


To calculate the quantity dr/dt, Crocco’s concept of 
pressure dependence of time lag has to be introduced. 
If the rate at which the liquid propellant is prepared for 
the final rapid transformation into hot gas is a function 
(p), then the lag r is determined by 


t 


By differentiating Equation [5] with respect to ¢, 


const [5! 


The concept of small perturbation from the steady state 
will now be explicitly introduced: Assume that the de- 
viation of the pressure p from the steady state value p 
is small. Then f(p) at the instant ¢ and f(p) at the in- 
stant ¢ — 7 can be expanded as Taylor’s series around p. 
By taking only the first order terms, 
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Here is the lag at the average pressure a 


[6] 


Then 


dlogf 
1+ (7 


now. 


— At — r)]... 


By combining Equations [4] and [6], the following _ 


equation is obtained 


d dp 
dz 


+ = — 8) + n[o(z) — oz — [7] 


where 


log 
~ \d log p/ r= 


[9] 


and 
7 = t/6,, 


If n is a constant independent of /, then f(p) is propor- 
tional to p". This is the form of f(p) assumed by 
Crocco. The present formulation of the problem is 
slightly more general in that f(p) is arbitrary and the 
value of n is to be computed by using Equation [8], and 
isa function of p. 6, is, of course, the gas transit time. 


Intrinsic Instability 


Crocco called the instability of combustion with 
constant rate of injection the intrinsic instability. If 
the injection rate is constant and not influenced by the 
chamber pressure p, then » = 0. Therefore the sta- 
bility problem is controlled by the following simple 
equation obtained from Equation [7], 


+ (1 — n)g(z) + ne(z — 6) = 0........ [10] 


e 


Now let 


Then 


s+(1 — n) + ne-& = 


This is the equation for the exponent s. 

Crocco determined the value of the complex number s 
by studying the set of two equations for the real and the 
imaginary parts of Equation [11]. However, if the 
point of interest is whether the system is stable or not, 
one can use the well-known Cauchy theorem w with ad- 
vantage. Let 


G(s) = 


Then the question of stability is determined by whether 
G(s) has zeros in the right half of the complex s-plane. 
This question itself can be in turn answered by watching 
the argument of G(s) when s traces a contour enclosing 
the right half s-plane. Specifically, let s trace clockwise 
the contour consisting of the imaginary axis and a large 
half circle to the right of the imaginary axis (Fig. 1). 
If the vector G(s) makes a number of complete clockwise 
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od 


CONTOUR TRACED BY THE VARIABLE S$ FOR THE SATCHE 
DIAGRAM OR THE NYQUIST DIAGRAM 


FIG. 1. 


revolutions, then that number is, according to Cauchy’s 
theorem, the difference between the number of zeros and 
the number of poles of G(s) in the right half s-plane. 
Since G(s) evidently has no poles in the s-plane, the 
number of revolutions of G(s) is the number of zeros. 
Hence for stability, the vector G(s) must not make any 
complete revolutions, as s traces the specified contour. 
Therefore the stability question can be answered by 
plotting graphically G(s) on the complex plane. This 
graph is, of course, the well-known Nyquist diagram. 

A direct application of this method to G(s) given by 
Equation [12] is, however, inconvenient for the com- 
plication caused by lag term e~* (8). M. Satche (5), 
however, proposed a very elegant and ingenious method 
of treating such a system with time lag: Instead of 


G(s), break it into two parts, 
G(s) = gi(s) — gels)... [13] 
where 
g(s) = 


The vector G(s) is thus a vector with vertex in g;(s) and 
its tail on go(s). The graph of g;(s) is the unit circle for 
s on the imaginary axis. For s on the large half circle, 
gi(s) is within the unit circle. The graph of g2(s) 
is the straight line (Fig. 2) paralleled to the imag- 
inary axis when s is on the imaginary axis. When 
s is on the large half circle, go(s) is a half of a large 
circle closing the contour on the left. A moment’s re- 
flection will show that in order for the vector G(s) not 
to make complete revolutions for any value of 6, the 


STABLE SATCHE DIAGRAM FOR INTRINSIC OSCILLATIONS; 
O<n< 


FIG. 2. 
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g2(s) contour must lie completely out of the g:(s) con- 
tour. That is, for unconditional intrinsic stability 


[15] 
n 2 


When n > '/s, the g;(s) contour and the g2(s) contour 
intersect. Stability is still possible, however, if for g2(s) 


FIG. 3. UNSTABLE SATCHE DIAGRAM FOR 
INTRINSIC OSCILLATIONS; n > 1/2 


within the unit circle (Fig. 3), gi(s) is to the right of 


g2(s). This condition is satisfied if ee 
(8\/2n — 1) > 

Or if 

where 

= 

V2n—1 


G(iw*) = 0. Therefore when 6 = 6*, ¢ has the oscilla- 
tory solution with the angular frequency w*. 

These results on intrinsic stability were obtained by 
Crocco. The present discussion with the Satche dia- 
gram, however, seems to be simpler. For the more 
complicated stability problem treated below with feed 
system and servo control, the solution is hardly practical 
without the Satche diagram. 


System Dynamics with Servo Control 


as Consider now a system including the propellant feed 
and a servo control represented by Fig. 4. In order to 
approximate the elasticity of the feed line, a spring load 
capacitance is put at the midway point between the 
propellant pump and the injector. The spring constant 
is to be computed from the feed-line dimensions.* 
Near the injector there is another capacitance controlled 


3 See the Appendix for details. 
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sure feed, a@ is zero. 


by the servo. The servo receives its signal from the 
chamber pressure pickup through an amplifier. If the 
feed system and the motor design are fixed by the de- 
signer, the question is whether it is possible to design 
an appropriate amplifier so that the whole system will 
be stable. Because there is no accurate information on 
the time lag of combustion, a practical design should 
specify unconditional stability, i.e., stability for any 
value of 6. 

Let mi be the instantaneous mass flow rate out of the 
propellant pump, and pp be the instantaneous pressure 
at the outlet of pump. The average flow rate must be 
m. The average pressure is fp. The pump character- 
istics can be represented by the following equation, 

[18] 

Po m 

If the time rate of change of mass flow is small, a is sim- 
ply related to the slope of the head-volume curve of the 
pump at constant speed near the steady-state operating 
point. For constant pressure pump or the simple pres- 
For conventional centrifugal 
pumps, ais approximately 1. For displacement pumps, 
a is very large. 

Let i, be the instantaneous mass rate of flow after 
the spring loaded capacitance, x the spring constant 
of the capacitance, and p, the instantaneous pressure at 
the capacitance. Then 


™m — m = 


where p is the density of the propellant, a constant. 

In the following calculation, the pressure drop in the 
line by frictional forces will be neglected. Then the 
pressure difference pp — p; is due to the acceleration of 
the flow only. That is 


fi = 


where A is the cross-sectional area of the feed line, a 
constant, and lis the total length of the feed line. Simi- 


larly, if pe is the instantaneous pressure at the control 
capacitance. 


m mi [22] 
t wb 
LINE CAPACITANCE 


Nin 


ate 
= 


SERVO-CONTROLLED LIQUID MONOPROPELLANT ROCKET 


MOTOR 


FIG. 4. 
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Since the control capacitance is very close to the injec- 
tor, the inertia of the mass of propellant between the 
control capacitance and the injector is negligible. Then 


where A, is the effective orifice area of the injector. A‘ 
can be eliminated from the calculation by noting that 
at steady state, the difference of pressures fp and A, 
or Ap is 


Equations [18] to [24] describe the dynamics of the 
feed system. By a straightforward process of elimina- 
tion of variables, a relation between m;, p, and C is ob- 
tained. To express this relation in nondimensional 
form, the following quantities are introduced, following 
the notation of Crocco: 


~ 2ap E= px, J= [25] 
and 


where 6, is the gas transit time given by Equation [9]. — f 
Then the nondimensional equation relating ¢, u, and x is: 


(P + 5) & JE d*) 
a) dz + 2 + 
(P+5)i + jen (P +3) 44 He 

2)\ 2 \ dz 

aJ E l JE} ad? J°E d3 
d ad E 1\ d3 


where z is the nondimensional time variable defined by 

[9]. 

l 


( 9) 
e+ (P+5) +4 


¢ = ae*, p = be*?, [29] 


By substituting Equation [29] into Equations [7,] [27], 
and [28], three homogeneous equations for a, b, and c 
are obtained. In order for a, b, c to be nonzero, the 


determinant formed by their coefficients must vanish. 


s? + 


This condition can be written as follows: 
J?E JE 
[s + (1 — n)) 


> 
+5){ s+ 
Jak (P +5) s+ + 
/ 1 
2 


aJ E (P+ 


= 0... [30] 


This is the equation for determining the exponent s. 
F(s) is now recognized as the over-all transfer function 
of the servo-control link. The complete system stabil- 
ity depends upon whether Equation [30] gives roots 
that have positive real parts 


Instability Without Servo Control oa 


The system characteristics without the servo control 
can be simply obtained from the basic Equation [30] 
by setting F(s) = 0. Let it be assumed that the poly- 
nomial multiplied into e~* has no zero in the positive 
half s-plane, as is usually the case. Then Equation 
[30] can be divided by that polynomial without intro- 
ducing poles in the positive half s-plane into the result- 
ant function. That is, for the Satche diagram, one has 
again 

G(s) = gi(s) — gs), gi(s) = 
gi(s) is thus again the “unit c circle.” 


ge(s) is now much 
more complicated : 


~ 


g(s) = — 


41 


(3 


The dynamics of the servo control are specified by 
the composite of the instrument characteristics of the 
pressure pickup, the response of the amplifier, and the 
properties of the servo. Since it is not the purpose of 
the present paper to discuss the detailed design of the 
servo control, the over-all dynamics of the servo control 
are represented by the following operator equation: 


where F is the ratio of two polynomials with the denomi- 
nator of higher order than the numerator. 

Equations [7], [27], and [28] are the three equations 
for the three variables ¢g, u, and x. Since they are equa- 
tions with constant coefficients, the appropriate forms 
for the variables are 
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The intercept of go(s), when s is pure imaginary, is given 
by setting s = 0 in Equation [31], i.e., 


1+a(P +5) 
[32] 
1+a(P+5)+ 


n 


Since all the parameters n, a, P are positive, the mag- 
nitude of g2(0) is now smaller than the magnitude of 
g(0) given by Equation [14] for the intrinsic stability 
problem. Thus the effect of the feed system is to move 
the go(s) curve toward the unit circle of g:(s) in the 
Satche diagram. For instance, for n = '/2, go(s) is just 
tangent to the unit circle for the intrinsic system with- 
out considering the propellant feed. But with the 
propellant feed system, go(s) contour will intersect the 
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g(0) = 


If the polynomial H(s) 


unit circle and the system will become unstable, for time 
lag 6 exceeds a certain finite value. The influence of 
the feed system is thus always destabilizing. This is 
further confirmed by considering the asymptote of g2(s) 
for large imaginary s, obtained from Equation [31]. 


That is 


gx(8) ~ — + 
Therefore, for large imaginary s, g2(s) approaches 
asymptotically a line parallel to the i aanginary axis ata 


. [33] 
n 


>> 1. 


distance 
n Jn? 


to the left of the imaginary axis. The effect of the 
feed system is again to move g2(s) toward the unit cir- 
cle. 

It is thus evident that for the parameter n near '/2 
or larger than '/», it would be impossible to design the 
system for unconditional stability. In the Satche dia- 
gram, g;(s) contour and g»(s) will always intersect with- 
out a servo control. 


Complete Stability with Servo Control 


]e+ 


Law 
aJE (P + 5) +5) . [34] 


2 
which multiplies into e~* in Equation [30], has no poles 
and zeros in the right half s-plane, then the occurrence 
zeros of the expression in Equation [30] in the right 
half s-plane can be determined from the Satche diagram 


with 


gi(s) 
and 
an 


8 1- J°E JE 


\ak ( 


H(s) = 


pa! 


As s traces the contour of Fig. 1, g:(s) is again a unit cir- 
cle. Therefore, if simultaneously the go(s) contour is 
completely outside the unit circle, there can be no root 
of Equation [30] in the right half s-plane. In other 
words, if the transfer function F(s) of the servo-control 
link is so designed as to place the g2(s) contour com- 
pletely out of the unit circle (Fig. 5), then the system is 
stabilized for alltimelags. 

As an example, take 


H(s).. . [35] 


n 


1+a(P+$) 
140 (P#$)+ 


WITH SERVO— 


~ 
~ 


WITHOUT SERVO — 


2P 


/- 


Jn? 


_ FIG. 5. SATCHE DIAGRAM FOR THE ORIGINAL AND FOR 
THE SERVO-STABILIZED SYSTEM 
Then without the servo control, the go(s) is 
1 (2s + 1) (2s + 38% + 9s + 6) 
2 s? + 3s? + 68 + 6 


Of primary interest is the behavior of go(s) when s is a 


g(s) = — 


pure imaginary number iw, w real. Thus 
gliw) = Ss 
1 (6 — 2lw? + (6 — 3w*) + — — e 
~ § (6 — 3w?)? + — w?)? 
1. (21 — 8w?)(6 — 3w?) — (6 — 2lw? + 4w!)(6 — 
(6 — 302)? + w6 — w?)? 

This contour for w > 0 is plotted in Fig. 6. It is evident 
that for sufficiently large values of time lag, the system V 
_ will be unstable. On the other hand, if the go(s) con- ” 
tour can be changed by the servo control to, say, 7 

g({s) = — 2 + 3X0 + 3) 3) 

Then, as plotted in Fig. 6, the new g2 contour is com- 
pletely’ outside of the unit circle of g,(s). Therefore T 
the system is now unconditionally stable. A straight- tl 
forward calculation from Equations [31] and [35] shows “a 
that the required transfer -F(s) for the servo te 
F(s) = la 
— 4.875 (s +.1.0528)(s? + 0.7164s + 2.6304) in 
s(s + 2)(s + 3)(s + 0.5332)(s? + 0.4668s + 3.7511) ur 
The servo link has thus the character of an integrating m 
circuit. If, with given response of the chamber pres- al 
sure pickup and of the servo for the control capacitance, T 
an amplifier could be designed to give an over-all trans- cr 
fer function close to that specified above, the com- W! 
bustion can be stabilized by such a servo control. sn 
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SATCHE DIAGRAM FOR THE ORIGINAL AND 
FOR THE SERVO-STABILIZED SYSTEM 
(g2(tw) without servo intersects the unit circle; 
g2(iw) with servo is outside the unit circle. Num- 
bers beside points are the value of w.) 


FIG. 6. 


As the second example, take 


1 3 


4 
Since a = 0, the feed pressure pp is thus constant with 
even variable flow of propellant. The case then corre- 
sponds to that of a simple pressure feed. Without the 
servo control, 
_ 1 (28 + 1) (288 + s* + 88 + 2) 

2 s3 + 2s? + 4s + 4 


g(s) = 


When s is pure imaginary, 
g{tw) = 
_ 1 (4 — 2w*)(2 — 17w? + 4w*) + w(4 — w?)(12 — 40?) 
2 (4 — 2w?)? + w%(4 — 
(4 — (12 — 4w*) —(4 — w*)(2 — 170? + 
(4 — Qu?)? + w4 — w?)? 


Ni— 


This contour of gs is plotted in Fig. 7. It is evident 
that without servo control the combustion will be un- 
stable for sufficiently long time lag. In fact, the sys- 
tem is even less stable than the system considered in the 
first example: It will become unstable at shorter time 
lag. The part of the g. contour near w = 2 is of special 
interest. Near w = 2, the contour comes so close to the 
unit circle of g; that if the value of time lag 6 is such as to 
make g; and g» for w ~ 2 very close to each other, then 
an almost undamped oscillation at # ~ 2 can occur. 
This critical value of 6 is evidently smaller than the 
critical 6 determined from the true intersection of go 
with the unit circle at w ~ 0.65. Such near instability at 
smaller values of time lag can be easily overlooked in 
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the analytic treatment of the stability condition by 
Crocco, and yet such possible instability should not be 
dismissed. This, perhaps, indicates the superiority of 
the present graphical method. 

For unconditional stability, gz should be displaced 
out of the unit circle, to, say, the same “‘stable’’ contour 
as in the first example. The required transfer function 
F(s) is calculated to be 

, ~- (s + 0.8126)(s? — 0.04337s + 2.6506) 

M(s) = 4.875 + 26 + + 4) 
The required servo link must then have the character 
of double integrating circuit. Furthermore, the trans- 
fer function has two purely imaginary poles at +27. 
This unrealistic requirement on the amplifier comes 
from the original feed-system dynamics and is due to the 
neglect of frictional damping in the feed line. In 
any actual system, the frictional damping in the feed 
line will remove these purely imaginary poles of the 
required transfer function F(s) and replace them by 
two complex conjugate poles. 


Stability Criteria 


In the preceding discussion of servo-stabilization, 
it is assumed that the polynomial H(s), Equation 
[34], has no pole or zero in the right half s-plane. This 
is, however, not necessarily the case. In general then, 
one should first investigate the number of zeros and 
poles of H(s) in the right half s-plane. To do this, it 
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SATCHE DIAGRAM FOR THE ORIGINAL AND 
FOR THE SERVO-STABILIZED SYSTEM 
P= %),,J=4, E = ‘/,a=0 

(g2(tw) without servo intersects the unit circle; 
-g:(iw) with servo is outside the unit circle. Num- 
; bers beside points are the value of w.) 


FIG. 7. 
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should be recognized that the polynomial in Equation 
[34] before the factor F(s) usually does not have zeros 
in the right half s-plane. Therefore instead of studying 
H(s), one can study the ratio of H(s) and that poly- 
nomial. That is, the number of zeros and poles of 
H(s) in the right half s-plane is the same as the number 


of zeros and poles of the following function < om 


order to obtain g,(s) and ge(s) as given by Equation 
[35], g zeros and r poles are introduced in the right half 
s-plane. The g poles of K(s) must come from F(s), 
since the polynomial in the denominator of Equation 
[37 | has no zero in the right half s-plane. Therefore the 
original expression in Equation [30] also has g poles in 


H(s) 


1+ K(s) = 


where 


- [36] 


K(s) 


According to the Nyquist criterion, the number of poles 
and zeros for 1 + K(s) in the right half s-plane can be 
found by plotting the Nyquist diagram of 1 + K(s) 
with s tracing the contour of Fig. 1. In fact, if 1 + 


K(s) or H(s) has r zeros and gq poles in right half s- 
plane then K(s) will carry out r—q clockwise revolutions 
around the point —1, as s traces the contour of Fig. 1. 
Hence the necessary information on H(s) can be ob- 
tained by plotting the Nyquist diagram of K(s). 

When one divides the Equation [30] by H(s) in 


Fig. 8. FUL 


CURVE FOR POSITIVE w; DOTTED CURVE FOR 
NEGATIVE 

(a) Nyquist diagram for K(s), with two zeros for 1 + K(s) 


in right half s-plane. (6) Corresponding stable Satche diagram. il (Continued on page 268) 
Ans 


[37] 


the right half s-plane. Hence in order for the original 
expression in Equation [30] to have no zero in the right 
half s-plane, g2(s) must make —q + (q — r) = -r 
clockwise revolutions around the unit circle. In order 
for stability to be unconditional, i.e., stable for all time 
lag, the ge(s) contour should never intersect the unit 
circle. Therefore the general unconditional stability 
criteria are, first, g2(s) contour completely outside of the 
unit circle; and, second, g2(s) making r counterclock- 
wise revolutions around the unit circle as s traces the 
conventional contour enclosing the right half s-plane. 
These are the criteria for stability with the Satche 
diagram. To determine r, one has to use the Nyquist 
diagram of K[s], Equation [37]. Thus the stability 
problem for the general case requires both the Satche 


diagram and the Nyquist diagram (Fig. 8). 


Concluding Remarks 


In the previous sections of this paper, the theoreti- 


eal possibility of completely stabilizing the combustion 


for any value of time lag by servo control is demon- 
strated. The great flexibility of electronic amplifier 
seems to indicate that this theoretical possibility 
can be always realized. On the other hand, with- 
out the servo link, unconditional stability is shown 
to be generally impossible. Therefore the concept 
of feedback servo is indeed a _ powerful tool in 
controlling the behavior of a time-lag system. It 
is to be realized, of course, that the proposed scheme 
is but one among many. No attempt is made here to 
give an exhaustive treatment of all possible schemes. 


The best scheme is certainly to be determined by de- 


tailed considerations on all aspects of the engineering 


- problem, such as the possibility of high-frequency 
acoustic oscillations which are not considered here. 


The main purpose here is to give a general discussion 
of the concept together with a suggested general method 


of analyzing the stability by the Satche diagram. 
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An analysis is presented for the determination of the 
natural frequency in longitudinal vibration of a static gas 
mass at ambient pressure contained in a rocket thrust 
chamber. The technique used is the Rayleigh-Ritz 
method of equating the kinetic energy to the work done in 
a conservative system, and determining the natural fre- 
quencies from the solution of the differential equations 
defining the energy relationship of the system. A com- 
parison of the expression derived from the theory and that 
of the one-quarter wave length oscillator is made and the 
relationship indicated. Tests were made to determine the 
resonant frequency of thrust chambers filled with air at 
normal temperature and pressure. Results from these 
tests with thrust chambers of three different configura- 
tions are presented in substantiation of the theoretical 
results. 


Introduction 


HE idea has been advanced that unstable combus- 

tion in the rocket motor is a manifestation of one or 
more of the possible modes of vibration of the gas mass 
contained in the chamber being excited into resonance. 
This theory, in order to be pursued, requires two pieces 
of information: (1) The frequencies associated with the 
various possible modes of vibration, which are best in- 
vestigated under nonrocket-operating conditions, i.e., 
with ambient air in the system; and (2) the average 
temperature and molecular weight of the gases through- 
out the combustion chamber under combustion condi- 
tions. Knowing the natural frequencies under am- 
bient conditions, they can then be corrected for the 
acoustic velocity related to the gas properties at the 
conditions existing in the combustion chamber. A com- 
parison of these calculated frequencies with those 
measured by experiment will then serve as a check of 
the theory. 

The normal modes of vibration of any enclosure are 
determined by the geometry which describes that en- 
closure. For the case of the rocket motor thrust cham- 
ber, its geometry lends itself best to the use of cylindri- 
cal co-ordinates for a proper definition. Consequently, 
there are three types of modes of vibration possible; 
namely those in the radial, circumferential, and longi- 
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tudinal direction; or symbolically, 7, 6, and y. Those 
vibrations in the r and @ direction are quite adequately 
covered by the classical acoustic theory describing the 
oscillations in a cylinder (1).4 However, the analysis 
for the longitudinal modes is not strictly applicable in 
view of the deviation of the rocket motor geometry from 
that of a constant area channel. 

The literature on acoustics is quite meager on the 
subject of resonances in variable area ducts. A. T. 
Jones (2) has treated the case of longitudinal resonances 
in a channel consisting of two cylindrical tubes of dif- 
ferent diameters connected by a conical section. Rich- 
ardson (3) has performed an analysis, assuming a cavity 
as a resonator, and has derived an expression for the 
fundamental mode of oscillation as a function of the 
cavity volume and neck dimensions. This is an exten- 
sion of the classical analysis (4) contained in acoustic 
literature pertaining to resonators characterized by 
having a throat with diameter small compared to the 
‘avity diameter. Richardson’s work is broader in that 
it is not limited by this requirement on the throat-to- 
cavity ratio. However, it was felt that none of these 
analyses could be applied satisfactorily to the rocket 
motor thrust chamber. 

At first it seemed possible to apply the theory de- 
veloped for resonances in a conical horn to this prob- 
lem. The rocket motor can be broken up into two sec- 
tions: (1) One section consisting of two horns with 
their throats adjacent to each other, thus accounting 
for the converging, diverging portion of the chamber, 
and (2) a second section consisting of a cylinder corre- 
sponding to the combustion chamber proper. The na- 
tural frequency of the entire system can then be deter- 
mined by matching impedances at the two interior 
boundaries in accordance with transmission line theory. 
However, it was considered more desirable to analyze 
the problem directly and develop a theory based on 
fundamental physical principles rather than adopt 
secondarily derived expressions. The theory as de- 
veloped is presented below.® 


Theory 


The Rayleigh-Ritz method (5) of determining the 
gravest mode of oscillation of a conservative system es- 
sentially consists of four steps: (1) Assume a form of 

‘ Numbers in parentheses refer to the References on page 268. 


* Based on an approach suggested by Dr. ‘Theodore von Kar- 
man. 
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elastic displacement of the system. (2) Compute the 
potential energy of the system at the maximum dis- 
placement. (3) Compute the kinetic energy of the sys- 
tem at the neutral position, i.e., position of zero poten- 
tial energy. (4) Equate the two expressions for energy. 
If the system being analyzed contains potential 
energy at the time of maximum kinetic energy or kinetic 
energy at the time of maximum potential energy, then 
it is necessary to use differences of potential and kinetic 
energy rather than absolute values of these two proper- 
ties. For example, if this technique were being used to 
analyze the behavior of a pendulum, then the difference 
of the potential energies would be the difference in ele- 
vation of the pendulum bob from the extreme of the 
swing to the neutral position, rather than the elevation 
of the bob above some arbitrary plane. In a similar 
manner, when considering the energies in a vibrating air 
column, one must use the difference in pressures be- 
tween the maximum and minimum compression and 
the difference in volume between the same points. The 
potential energy, U, is thus equal to the work done by a 
mass of gas undergoing expansion from its maximum to 
minimum pressure. This is then one half the change in 
pressure times the change in volume, which can be writ- 
ten as follows (7): 


and since Ap / Ap = a®, where a = 
sumed constant) 


Consider the expansion of a mass of gas, dm, through 


a duct of variable area; then the relations depicted in 
Fig. 1 hold where 


= displacement of the given mass of gas Rea! 


Volume before expansion = SAz. Volume after expansion 


(s + )(ar + ar) 


OS 
or Oz 


ll 


Sar + Sax + ar 98 + gar 


Rearranging the above expression and ns the 
higher order terms, we obtain 


--f; apa (5 )am ite 


which becomes, after performing the indicated dif- 


acoustic velocity (as- 


S = cross-sectional area at position zx 7 

Az = thickness of slab of the given mass of gas a — 
p = gas density 

V = volume of given mass of gas 


1 o(éS) 


~Saz [1 


Then 
1 0 (éS) 
AV _ ae + — SAz 
SAz 
Then substituting into [1] - 
The density p at any position z in the duct is a 
p= m+ =m (1-247) 
Po 
” Po S oz 
Then substituting into [2] _ 
q 


The second term is dropped due to its higher order and 
consequent insignificance relative to the first term. 
(This is equivalent to assuming nearly incompressible 
flow, that is, small oscillations. ) The final expression for 
the potential energy then becomes 


The kinetic energy of the system is 


l 
12, 258, 


ferentiation, and using the relationship t 
oe _av —_ Assuming plane wave motion through the duct, and 
— assuming that the velocity will vary sinusoidally with 
1 ap (AV time; then at any x 
2 Ap ) - fi Ap v = usin wt = 


for any normal mode. The displacement is: 


Vo 
= — — al 
w 


In establishing the limits of integration, we note that 


the kinetic energy will be a maximum at time t = -S 2 


2 w 


i.e., when the velocity is the greatest. It is this inte- 
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FIG. 1 VARIABLE AREA DUCT 
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grated maximum 7’ which is equated to the potential 
energy U. It follows, then, that 


Vo a/2w 
Em = — — COS wl —— 
\0 w 
Um? = (wm)? 


4 
The expression for the maximum kinetic energy be- 
comes 


The variation in p may again be neglected as shown 
in the development for the potential energy expression. 
Summarizing then, the expressions for the maximum 
kinetic and potential energies are 


0 


2 


where the limits of integration extend over the length 
of the cavity for the geometry shown in Fig. 2 


FIG. 2. THRUST CHAMBER GEOMETRIES 


It remains now to relate S and £,, to x in order to per- 
form the above indicated integrations. 

Since both S and é,, are functions of x, the equation 
becomes complicated mathematically. Therefore, as- 
sume a series solution 


S SA == 8 
tm = ~ n 2 i) [ ] 
n=1 


where the n’s define the characteristic values as deter- 
mined by the boundary conditions of the particular 
configuration. The geometry under consideration has 
an open end,® and consequently only odd values of n 
are valid. Considering only the first and third har- 
monics, 


6 6 The effect of this assumption is discussed later on in the paper. 
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The functions of Sé,, shown in Equations [6] and [7] 
are 


(Stn)* = [ sin (57) + sin (5 5) sin + 
— 
lz - EB A; cos + 
A,A; cos (5 ;) cos i) + 


Then, substituting [9] into [6] and introducing the 
constants S, for the minimum duct cross-sectional area, 


and L for the length of the duct 


2 
op gem, + + 11] 
2 St 


So (7) (7) 
f i) (= i) d (7) [12] 
Sv (27) 4 (7) 

and substituting [10] into [7] 


where 


The integrals K, .... Ke are dimensionless and can be 
evaluated by graphical integration or by numerical 
methods for the particular geometry being considered. 
The system is conservative and therefore 7 = U. 
Since 7 — U = 0 must be true for all values of the two 
A, and Ag, it follows: 


ar -U)_xT-U)_ 
oA; 


Then, conducting the above indicated operations, we 


i 


Ke 


obtain 
fe 
Aik, + 3 (57) [15] 
L 
= pow? [A, Ke + A3K3] 


2 2 
[3 Ai Ks; + 9 AK | [16] 
t 


From [15] and [16], obtain 


Ped 
] 
oi 
> 
Cc 
. 
sin 
2 
2 
65 
: 


The ratio of the two constants, A;/A3, can thus be ob- 
tained from [17] or [18]. It is seen that two different 
values of A,/A; are obtained. As will be seen later, two 
frequencies are obtained. The equation » = % sin wt 
which was utilized in the derivation to establish the 
kinetic energy implies a single frequency along the chan- 
nel for any normal mode. However, there will be dif- 
ferent frequencies for the different normal modes. Each 
value of A,/A; describes a different normal mode. (The 
obtaining of multiple frequencies is discussed in greater 
detail in Reference (8)). The calculated values of 
A,/A; differ by large amounts. One value of Ai/A; 
may be considered to define fi, while the other defines /3. 

At this point, it is convenient to introduce the sub- 
stitution, (ra/2wL)* = 


Then forming the fraction [17]/(18] 


Ky — 0K, 


30Ks — ¢ 
K, — 30K; 


A quadratic expression in @ is obtained which upon 
solution becomes 


18(KiKe — K;?) 


V(6K2K; — — — 36(K.Ks — Ks*)(K3K, — 
18(K,K, — K;?) 


Two values of @ are obtained, 6, and 62, which thus de- 
fine two values of w. Then, recalling the substitution 
made earlier, (ra/2wl)? = 6, and solving for w (and neg- 
lecting the negative value), obtain 


we ra 


a 
or since, w = 
a a 
= fs = ——........... 21] 
fi 4LV/0,  4LV/ 


Note that @ is a constant which is determined solely 
by the geometry being considered, When @ = 1, the 
expression for the '/,\ oscillator is obtained. Therefore, 
+/@ can be interpreted as a correction factor to the 
chamber length, L, which will make the '/, \ formula 
applicable. 

Furthermore, the value of +/@ is independent of L, 
provided the dimensionless function S/S, = F(z/L) is 
the same for the configurations being considered. This 
function may be used, therefore, as a definition of 
acoustic similarity. 


Discussion 


1. Theoretical 


Calculations based on the analysis developed above 
were made to predict the natural frequency at ambient 
conditions of three thrust chambers of different con- 
figuration. Two are of the conventional rocket motor 


design consisting of a cylindrical section followed by a 
converging-diverging nozzle. A third is characterized 
by a conical section to the throat of the nozzle, followed 
by a diverging section. For purposes of this discussion, 
these chambers will be designated as Chambers A, B, 
and C, and they are shown schematically in Fig. 2. Re- 
sults of the calculations are presented in Table 1. 


TABLE 1 
Chamber A B C 


fi Calculated by Equation [21] 127 96 102 
fs Calculated by Equation [21] 576 8359 
fi: Calculated by Equation [23] 142 
fs Calculated by Equation [23] 573 
f, Experimental 125 98 102 


These calculations were carried through with only two 
terms in the series expansion, as shown previously, thus 
producing values for the fundamental and first overtone. 
It is of interest to note that the overtone is not a har- 
monic of the fundamental. This is in contrast to the 
condition for constant area pipes, where the natural 
frequencies in longitudinal resonance form a harmonic 
series. 

As can be seen from the analysis, a specific applica- 
tion of the method involves the evaluation of six in- 
tegrals when two terms are used in the series expansion. 
If three terms are used, there are twelve integrals de- 
veloped, or generally n(n + 1) expressions need to be 
evaluated where n is the number of terms considered in 
the series expansion. Each of these expressions re- 
quires either graphical or numerical calculation meth- 
ods for their solution—procedures which are admittedly 
tedious and time-consuming. Attempts to obtain a 
solution which is completely analytical have been un- 
successful to date. 

In an attempt to simplify the method, an analysis has 
been conducted with only one term in the expansion for 
the particular tone required, i.e., 


Resulting expressions for the frequencies are thus ob- 
tained which require only two integrals to be evaluated, 
or specifically, 


hea VRE [23] 


where the subscripts are those identifying the integrals 
in [12] and [14]. Calculations carried through with 
this method give values for the fundamental which are 
12 per cent higher than when both terms are used in the 
expansion. There is negligible deviation for the third 
harmonic. The results obtained from experimentation 
(which will be discussed later and are shown in Table 1) 
show satisfactory agreement with the frequencies cal- 
culated using two terms in the series expansion. There- 
fore, to obtain a satisfactory degree of accuracy, it is 
recommended that the lengthier set of calculations be 
performed. 

Attention must be directed to one point. In the 
analysis, an assumption is made that plane wave motion 
exists throughout the chamber. Obviously, this as- 
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sumption is not rigorously correct since the sloping 
walls of the converging-diverging nozzle must distort 
the plane wave. However, it seems reasonable that 
the effect of the converging section on the plane wave 
front emanating from the cylinder will be offset by the 
effect of the diverging section on the resultant spherical 
wave. These two effects conceivably may combine to 
make valid a plane wave analysis for the over-all wave 
propagation. Certainly, the agreement obtained be- 
tween the experimental and theoretical values show 
that the assumption of plane wave motion does not in- 
troduce enough error to make this approach invalid. 


2. Experimental Results 


Experimental results obtained from testing with the 
three afore-mentioned chambers are also included in 
Table 1. Values for the third harmonic cannot be pre- 
sented due to the inconclusiveness of the experimental 
data obtained. The experimental technique consisted 
of vibrating the air mass contained in the chamber by 
means of a high impedance sound source driven by an 
audio oscillator. Resonances were detected by noting 
the increased output of a microphone located in the 
chamber as the driving frequency was swept through 
the audio range. Care was taken that the microphone 
be sufficiently small to avoid distorting the sound field 
present in the chamber. As can be seen by comparison 
with the calculated values, satisfactory agreement exists 
between experimental and theoretical values. 

_ Comparison with Organ Pipe Theory 7 

It is interesting to compare the theory developed for 
the organ pipe with that developed in this paper for 
variable area ducts. One difference has already been in- 
dicated earlier, namely, the higher frequencies for the 
variable area duct are not integral values of the funda- 
mental, as is the case for the constant area duct. It has 
been shown in Equation [21] that the final expression 
for the frequency is related to the '/, wave-length ex- 
pression by a geometrical factor, +/@. Then, taking the 
chamber as a constant area duct with the same length 
leads to a frequency of 152 eps for Chamber A and 108 
cps for Chambers B and C. A comparison of these fre- 
quencies with those calculated from the theory results 
in values for +/@ = 1.2, 1.13, and 1.06 for Chambers A, 
B, and C, respectively. It is obvious that, although 
Chambers B and C are of equal length, thus possessing 
the same frequencies when the '/, \ expression is ap- 
plied, different values for +/6 would be obtained as a re- 
sult of their considerably different internal geometries. 

The theory pertaining to oscillations in a tube, from 
which the !/, \ formula is derived, assumes that a true 
antinode exists at the open end. This is not strictly 
correct, however, in that the air beyond the open end is 
also in vibration, which consequently makes the effec- 
tive length greater than the actual length. Experi- 
menters in acoustics have observed that end corrections 
to the physical length are required in order to reconcile 
the '/, \ formula with test results. Values used vary 
from .6R to .8R (where R is the radius of the tube) de- 


SEPTEMBER-OCTOBER 1952 


pending upon the nature of the flange at the exit. 
Therefore, for a proper test of the applicability of the 
1/,\ expression to the longitudinal resonances in rocket 
motor thrust chambers, it is necessary that end correc- 
tions be applied. This has been done for the three 
chambers, using an intermediate value of 0.7 of the 
radius of the exit, as the length to be added to the geo- 
metrical length of each chamber. The resultant fre- 
quencies obtained are 135 cps for Chamber A and 90 
cps for Chambers B and C. For comparison, Table 2 
shows the values discussed above. _ 


TaBLE2 
Chamber A B C 
fi Calculated by 1/, \ formula 152 108 108 
V6 1.20 1.13 1.06 
f; Calculated by 1/, \ formula 
with end correction 135 90 90 


It may be argued that end corrections should be ap- 
plied as well to the expression for the frequency derived 
in this paper. However, a value for the end correction 
for a flared tube has not been firmly established. Two 
experimenters in acoustics, Mokhtar and Messih (6), 
have shown recently in their work on conical pipes that 
the end correction is less than that of a cylindrical pipe 
having the same radius as the mouth, without advanc- 
ing any quantitative results. A frequency obtained by 
the Rayleigh-Ritz method is higher than the true fre- 
quency, the degree of error depending on the lack of 
similarity between the assumed and true forms of 
elastic displacement. Neglect of the end correction also 
increases the calculated frequency. It would appear, 
therefore, that the calculated frequency of Equation 
[21] should be higher than the measured frequencies. 
However, the error in the Rayleigh-Ritz method is 
usually small. Den Hartog (5) analyzes the frequency 
of a vibrating string by this method, utilizing three as- 
sumed forms of elastic displacement. The resulting rela- 
tive frequencies are as follows: 


RELATIVE 
Form or DISPLACEMENT FREQUENCY 
Sinusoidal (true form) at? 1.00 
Parabolic 1.006 
Triangular 1.102 


It is believed that the assumption of sinusoidal motion 
for gas vibrations cannot be far in error. The agree- 
ment between the experimentally obtained and theoret- 
ically calculated results indicates that end corrections 
for the three chambers considered would be small. A 
possibility does exist that the lack of end correction and 
the assumption of plane waves in the convergent and 
divergent sections of the duct give rise to compensating 
errors. 


Application to Operating Rocket Thrust 
Chambers 


It remains now to relate the preceding theory which 
was developed for the nonoperating ambient condition 
with the operating thrust chamber. There are dif- 
ferences between the two conditions which must be 
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considered further, the most apparent being as fol- 
lows: 

(1) A nonplanar sonic boundary exists in the throat 
of the rocket nozzle. 

(2) The gas contents of the thrust chamber are 
flowing at high velocity, whereas in the derivation no 
through-flow velocity was assumed. 

(3) The distribution of the RT function throughout 
the chamber is not constant and therefore the acoustic 
velocity is not constant. 

(4) Supersonic flow exists near and downstream of 
the throat. This will demand consideration of compres- 
sibility effects. 

(5) The oscillations of chemical energy release and 
entropy with time as a result of combustion under oscil- 
lating pressure will affect the conservation of the poten- 
tial and kinetic energy within the system, one of the 
basic assumptions of the derivations. 

An extension of the theory developed should include 
the factors listed above. It is hoped to extend the above 
theory and to present data on the experimental results 
of firing tests on the thrust chambers, the dimensions of 
which were used to calculate the values presented in 


1,8 

1 “Vibration and Sound,” by P. M. Morse, McGraw-Hill 
Book Company, Inc., New York, Toronto, London, 2nd edition, 
1948, pp. 398-401. 

2 “Resonance in Certain Non-Uniform Tubes,” by A. T. 
Jones, Journal of the Acoustic Society of America, vol. 10, 1939, 
pp. 167-172. 

3 “The Amplitude of Sound Waves in Resonators,” by E. G. 
Richardson, Proc. Physical Society of London, vol. 40, December 
1927 to August 1928, pp. 206-220. 

4 “The Theory of Sound,” by Rayleigh, Dover Publications, 
New York, Ist American edition, 1945, vol. II, pp. 170-172. 

5 ‘Mechanical Vibrations,’ by J. P. Den Hartog, McGraw- 
Hill Book Company, Inc., New York and London, 3rd edition, 
1947, pp. 178-185. 

6 “The Acoustic Characteristics of Conical Pipes,” by M. 
Mokhtar and G. Messih, Physical Society of London, vol. 62, 1949, 
pp. 793-799. 

7 “The Theory of Sound,” by Rayleigh, Dover Publications, 
New York, Ist American edition, 1945, vol. II, p. 17. 

8 “Vibration Problems in Engineering,” by S. Timoshenko, 
D. Van Nostrand, Inec., New York, 2nd edition, 1947, pp. 370- 
376. 


Servo-Stabilization of Combustion in Rocket Motors 
(Continued from page 262) 


It is of interest to point out that stabilization by 
servo control is only one phase of the general concept of 
feedback link. The opposite case of destabilization 
could be of importance also. For instance, consider 
the so-called valveless pulsejet. It is not always 
possible to operate the engine with the desired pulsa- 
tion. With a feedback servo linking the combustion 
chamber pressure pickup through an amplifier to the 
fuel line, the system can be destabilized at the desired 
operating frequency and thus operate the engine at that 


frequency of pulsation. This application of servo- 
destabilization gives the valveless pulsejet a new 
flexibility and an extended range of operation. There- 
fore it seems worth while to explore carefully all possible 
applications of feedback control to systems with time 


lag. 


APPENDIX 
, Calculation of Parameters J and E 


If L* and c* are the characteristic length and the 
characteristic velocity of the motor, and if 7, is the 
chamber temperature, R the gas constant, the transit 
time is 0, = L*c*/RT.. 

To calculate J and EF defined by Equation [25], it is 
more convenient to use the average propellant velocity 


v in the feed line. Thus m = pAv. 
rhus, according to Equation [25] 
1 l 
ger 
A consistent set of units would be p in slugs per cubic 
foot, v in feet per second, / in feet, 8, in seconds, and 
Ap in pounds per square foot. 

If dis the diameter of the feed line, h its thickness, and 
FE’ Young’s modulus of the tube material, then x, the 
change in volume of the feed line per unit rise in pres- 
sure, is 

d\? , 
x = (5) a/ 
a 


Therefore Equation [25] gives 


A consistent set of units would be Ap in pounds per 
square inch, E’ in pounds per square inch, / in feet, 6, 
in seconds, and v in feet per second. 
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A functional border between atmosphere and space is 
defined as a level at which the atmosphere fails as a sup- 
porting medium, and space-equivalent conditions begin. 
Depending upon a particular kind of function the corre- 
sponding limit is located at a certain altitude. The major 
functions of the atmosphere for man and craft are the fol- 
lowing: Contributing to respiration, preventing boiling of 
body fluids, sustaining combustion of fuel, absorbing 
heavy primaries of cosmic radiation, absorbing solar UV- 
radiation, supplying aerodynamic lift, supplying diffuse 
daylight, absorbing meteors, interacting thermally with 
the craft, and interfering by air drag over long periods of 


time (permanence of satellite orbits). Depending upon 


the nature of a particular function, the functional borders 


so defined are more or less extended regions. The various 


functional borders of space lie in the region between 10 and 
120 miles of altitude. The significance of the above men- 
tioned factors for manned rocket flight is discussed with 
special emphasis upon problems of aero-medical and 
space-medical nature. The use of the term ‘‘aeropause”’ 
for the border region between atmosphere and space is pro- 
posed. 


HE upper boundary of the atmosphere is com- 

monly identified with that region of the exosphere 
where the upper most geophysical phenomena, namely 
the highest aurora, are occasionally observed. In 
terms of this concept, the limit of the terrestrial at- 
mosphere is located at about 600 miles above the sur- 
face of the earth. The peak of the highest rocket 
trajectories attained so far—250 miles—lies yet within 
the boundaries of the atmosphere. For all practical 
purposes of rocket engineering, however, the at- 
mosphere ceases to exist at an altitude of 110 — 
miles. Unmanned rocket craft are routinely r reaching 
beyond the physically effective regions of the at-— 
mosphere, and manned flights in the border region’ 
of the atmosphere are being made in experimental 
rocket airplanes. Flights at very high altitudes, i.e., 
in excess of 10 miles, are still short in duration; never- 
theless, the technical means available in aviation of 
today make it necessary to investigate the medical 
and psychological problems peculiar to flight in the 
border area between the terrestrial atmosphere and 
space. 

Since the geophysical concept of the borders of the 
atmosphere is of little use in aviation medicine and 
space medicine, a new concept must be found. This 
new concept must be adapted to the uses of rocketry 
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Manned Flight at the Borders of Space 


The Factor of Manned Ro Rocket t Flight 


and aviation medicine, in other words, it must be based 
on the functions which the atmosphere fulfills for craft 
In a series of papers, Haber (1, 2)? and 
Strughold et al. (3) attempted to develop such a new 
concept by introducing the “functional borders be- 


and men. 


tween atmosphere and space.”’ A functional border 
is defined as a certain level above the ground at which 
the atmosphere fails as a supporting medium. De- 
pending on the nature of the various functions, the 
borders arelocated at different altitudes. They are found 
inthearea between 10 and 120 miles of altitude. In these 
regions of the atmosphere the conditions of conventional 
aviation gradually blend into those of actual space 
flight. 

From this point of view the following functional 
borders can be listed: 


TABLE | 
Funetion Altitudes, miles 
1 Contributing to respiration 10 
2 Preventing boiling of body fluids — - 12 
3 Sustaining combustion of fuel 13-15 
4 Absorbing heavy primaries of cosmic ra- 
diation 13-23 
5 Absorbing solar ultraviolet radiation be- 
ot tween 210 and 300my (Hartley band of 
Os) 22-28 
6 Supplying aerodynamic lift 50-60 


7 Supplying diffuse daylight 
8 Absorbing meteors 

9 Interacting thermally with the craft 
10 Interfering by air drag over long periods 
of time (permanence of satellite orbit) 120 


In addition to these data it may be mentioned that 
the presence of ozone above the 8-mile level can re- 
sult in toxic concentrations of this gas in the cabin 
air, if the pressurization of the cabin is maintained by 
compressing ambient air. Fortunately, ozone is easily 
intercepted by means of simple filters. 

Of course, the functional borders so defined are more 
or less extended regions. Especially the functions men- 
tioned under 6, 9, and 10 are dependent upon the 
velocity of the craft, and the altitude data given are 
related to a velocity of the order of 5 miles/second. 
This velocity must be attained in order to establish 
a craft in a permanent satellite orbit around the planet. 

In the following, those functions of the atmosphere 
that have a bearing on the human factor of manned 
rocket flight will be discussed briefly. 


3 Numbers in parentheses refer to the Bibliography on page 276. 
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A The Function of Contributing to Respiration 


Quantitative information as to the extent to which 
the atmosphere can maintain the oxygen supply of the 
body is obtained through experiments on explosive 
decompression. In experiments of this kind, a test 
subject is transposed, within fractions of a second, 
from an environment of normal or slightly reduced 
oxygen pressure to a state of acute oxygen need such 
as would be found at higher altitudes. If a subject 
is brought abruptly from normal oxygen pressure to 
that found at 5 miles the first psychophysiological dis- 
turbances can be observed after about two minutes; 
after another minute the subject becomes entirely 
helpless, and loses consciousness. The interval during 
which he is still capable of acting is called the ‘time of 
useful consciousness” after Armstrong (4), or ‘time 
reserve” after Strughold (5). This characteristic 
time span is reduced to about 80 seconds at 6 miles, and 
it is further reduced with increasing height, and finally 
reaches an asymptotic value of about 15 seconds at a 
height of 9 to 10 miles. Fig. 1 shows the time of useful 


40 80 120 


SECONDS 
TIME OF USEFUL CONSCIOUSNESS AT VARIOUS LEVELS OF 
ALTITUDE 


Fig. 1 


consciousness as a function of altitude according to 
measurements carried out by Luft, Clamann, and Opitz 
(6). In experiments made in the low-pressure chamber 
at simulated altitudes up to 12 miles, the time of useful 
consciousness remained constant at an average of 15 
seconds. There is no reason to assume that the mini- 
mum time span of 15 seconds will become any shorter 
at still greater heights. The height of the critical level 
in the atmosphere—about 10 miles where we first ob- 
serve the minimum time reserve—can be derived from 
an analysis of the alveolar gases. The constitution 


of the air in the lungs differs from ambient air. In 
average atmospheric air of 760 mm Hg, about 580 mm 
are provided by Ne, 150 mm by O:, 23 mm by water 
vapor, and 7 mm by A. The rare constituents of free 


air, such as the other inert gases and COs, are omitted 
in this listing. The approximate values for alveolar 
air at sea level are: 560 mm No, 106 mm O:, 47 mm 
water vapor, 40 mm CO., and 7 mm A. At higher 
altitudes the partial pressures of the individual gases 
in free air fall off in exactly the same rate as the total 
pressure, since the atmosphere is being kept mixed 
thoroughly throughout. (The gravitational separation 
of atmospheric gases existing in the 22-30 miles layer 
and above, Ref. 7, can be omitted in this discussion.) 
Alveolar air, however, behaves differently. If air of a 
lowered total pressure is breathed, the partial pressures 
of nitrogen and oxygen in the lungs are likewise lowered 
proportionally; the partial pressures of water vapor 
and carbon dioxide, however, remain practically 
constant at about 80 to 87 mm Hg. This is because 
water vapor and carbon dioxide inside the lungs origi- 
nate from the blood through continuous evaporation and 
exhalation. At increasing heights, water vapor and 
carbon dioxide occupy an ever-increasing percentage 
of the alveolar air. Obviously, the air inside the lungs 
will consist exclusively of these two components as soon 
as the pressure of the ambient air falls to 80 to 87 mm, 


a value which is reached at an altitude of about 9 to 10 


_ water vapor and carbon dioxide. 
atmosphere to contribute to man’s respiratory needs 


miles. At this altitude no additional oxygen can enter 
the lungs because their capacity is exclusively claimed by 
The function of the 


comes to an end at an altitude of about 10 miles. 


_ B The Function of Preventing Boiling of Body Fluids 


At a somewhat greater altitude the barometric pres- 
sure of the atmosphere can no longer prevent the body 
fluids from boiling. The boiling point of body fluids 
is reached at an altitude where the barometric pressure 
becomes equal to the water vapor pressure at body 
temperature, i.e., 47 mm Hg at 98.6 F. The critical 
level where boiling of body fluids sets in, lies at about 
12 miles of altitude. Armstrong (8) has studied the 
boiling of body fluids in a low pressure chamber 
by exposing animals to a barometric pressure below 


47 mm Hg. The first manifestations of boiling of 


body fluids appeared as rapid bubbling of the saliva 
around the mouths of the animals. After a delay of 
several seconds, ballooning of the animal skin was 
observed which is attributed to the formation of water 
vapor bubbles in the soft tissue. Boiling of blood in 
the lungs must also be expected upon lowering the baro- 
metric pressure below the critical value of 47 mm Hg. 


C The Function of Absorbing Heavy Primaries of Cosmic 

Radiation 

Originally, the cosmic ray particles were thought 
to consist exclusively of protons, i.e., the nuclei of 
hydrogen atoms. In 1948, Freier, Lofgren, Ney, and 
Oppenheimer (9) presented evidence to the fact that 
primary cosmic radiation contains also nuclei of heavier 
atoms. Bradt and Peters (10) determined the abun- 
dance ratio of different nuclear species in primary 
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cosmic radiation. Their results indicate that the 
following ratios exist: H:He:Heavier nuclei = 
79:20:1. This frequency distribution can be expected 
to correspond approximately to the cosmic abundance 
of the chemical elements. Owing to the tremendous 
energy of the cosmic primaries (~ 2 BeV per nucleon), 
they are characterized by an enormous energy dissi- 
pation in passing through matter. This energy dis- 
sipation, which is dependent upon velocity and nuclear 
charge of the individual particles, determines the pene- 
trating power of the particles. Consequently, the 
lower limits to which a certain percentage of the pri- 
maries can penetrate exhibit a certain stratification 
according to atomic number and energy of the particles. 

About one half of the primary protons penetrate 
down to an altitude of approximately 12.5 miles. 
The lower limit of the heavy primaries is found at about 

12 miles. Above this level their number increases 
considerably reaching a saturation at an altitude of 
about 23 miles. The air blanket above this latter 
level has no appreciable effect on the density and nature 
of the incoming particles of cosmic radiation. 

Because cosmic rays are a highly ionizing type of 
radiation, their possible hazard to man exposed to 
them at greater altitudes must be investigated. The 
biological effectiveness of cosmic radiation was studied 
theoretically by Krebs (11), Schaefer (12, 13, 14), 
and Muller (15). Krebs gave special emphasis to the 
possible effects of the so-called star phenomenon. A 
primary particle entering the atmosphere from outer 
space carries such a high energy that specific hitherto 
unknown phenomena are produced in the processes of 
collision between a cosmic ray particle and a target 
nucleus. The cosmic ray particles create an explosion 
or “evaporation” of nuclei resulting in a complete 
disintegration of nuclei. The fragments shooting away 
from the site of the nuclear event imprint themselves 
on the emulsion of a photographic plate producing star- 
like formations, A typical star is shown in Fig. 2. 
According to Hornbostel and Salant (16), the frequency 
of stars with 4 or more prongs is about 2000 per cubic 
centimeter per day, at an altitude of 17 miles. Krebs 
estimates that the total of the prongs produced in the 
tissue of an individual exposed to cosmic radiation 
outside the air blanket is equivalent to a concentration 
of 6.2 X 10~-!* grams radium elements per cubic centi- 
meter. Comparing this figure with the natural occur- 
rence of radioactive materials found in human tissue, 
Krebs concludes that “the radiation delivered by the 
cosmic ray stars at the top of the atmosphere comes 
very close to the amounts of radioactive energies which 
today are considered to be in no way harmless to tissue 
... the possibility for biological effects of cosmic radia- 
tion becomes, to a certain degree, reality.” 

Hermann J. Schaefer has evaluated the present-day 
knowledge of cosmic radiation in terms of the hazard 
to health. He was the first to point out that the heavy 
particles of cosmic radiation must be expected to be 
particularly effective, even though their number is 
small compared to that of primary protons. A detailed 
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analysis by Dr. Schaefer of the biological aspects of 
the heavy primaries of cosmic radiation is given else- 
where in the JouRNAL (see pages 277—283 in this issue). 


D The Function of Absorbing Ultraviolet Solar Radia- 
tion 

The existence of ozone in the atmosphere in quan- 
tities sufficient to entirely block the passage of ultra- 
violet radiation from the sun was first assumed by Hart- 
ley (17), the discoverer of the absorption band of ozone 
between 210 and 300 mu, that bears his name. The 
absorption of ozone begins at about 300 mu, and soon 
becomes so effective that the solar spectrum is cut off 
abruptly at 286.3 my, as reported by Goetz (18). 
The maximum of ozone absorption lies at 255 my, 
where the absorption coefficient ai) was found to be as 
high as 123/cm. A layer of ozone at standard pressure 
and temperature, having a thickness of 0.002 cm, 
reduces the flux density of an incident beam of ultra- 
violet light to about one half. This shows that ozone 
is a more effective absorber than all metals. 

The bulk of atmospheric ozone is found in the layer 
between 9 and 25 miles of altitude. The mean total 
amount of ozone corresponds to a layer of 3-mm thick- 
ness at STP. The existence of this small amount 
of ozone in the atmosphere is nevertheless of great 
biological importance. It provides an extremely effec- 
tive blanket against those parts of ultraviolet radiation 
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FIG. 2. TYPICAL COSMIC-RAY STAR 


(From H. L. Bradt and B. Peters, Ref. 10. Courtesy of Dr. M. F. 
Kaplon, Dept. of Physics, The University of Rochester, Rochester, 
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_ that have a pronounced effect upon the skin and the 

eyes by causing erythema of the skin and conjunc- 
tivitis. Fig. 3 shows the erythema producing effec- 
tiveness in arbitrary units as a function of wave 
length. The values given were adopted by the In- 
ternational Illumination Commission in 1937 (19). 
On the basis of these data, and by use of recent results 
of ozone investigations, Buettner (20) has estimated 
that erythema of the skin within and above the ozone 
layer is produced 10 to 50 times as fast as at sea level. 
Protection against ultraviolet radiation, therefore, is 
mandatory at altitudes in excess of 15 miles. Fortu- 
nately, window materials capable of absorbing these 
kinds of radiation are available. Upon massive irradi- 
ation by ultraviolet rays, however, most materials be- 
come turbid and discolored. Quartz, for instance, be- 
comes brownish if exposed to large amounts of UV-radi- 
ation of a wave length smaller than 250 mu. 


E The Function of Supplying Aerodynamic Lift and Drag 


It appears as though the final disappearance of aero- 
dynamic lift at greater altitudes is of interest only to 
the engineer. It can be shown, however, that the 
absence of aerodynamic lift and drag will be responsible 
for a medical problem of far-reaching consequences. 
For, in the absence of aerodynamic and propelling 
forces, the crew of a high-flying rocket craft will be 
transposed into the state of weightlessness. 

Following the reasoning of Haber (21), weight of a 
body can be identified with the vectorial sum of all 
external forces acting on the body. The same defini- 
tion can be applied to the concept of “apparent weight”’ 
or “g-forces’” used in aviation medicine which are 
usually defined as the sum of the force of gravity and the 
forces of inertia inherent to accelerated motion, accord- 
ing to the equation 


where W, is the apparent weight of the body, F, the 
force of gravity, and F; the forces of inertia. Using 
the principle of d’Alembert 

[2] 


3. BIOLOGICAL EFFECTIVENESS OF UV RADIATION 


where F, is the vectorial sum of all external forces, 
the apparent weight or g-forces can simply be written as 

The external forces comprise friction, aerodynamic 
lift and drag, propelling forces from an engine, and all 
elastic forces involved in the situation of a body being 
fully or partially supported. As can be derived im- 
mediately from Equation [3], a body becomes weight- 
less if it coasts freely outside the range of the mechani- 
cally effective parts of the atmosphere. 

Owing to the decreasing support which the atmos- 
phere provides at increasing altitudes, states of reduced 
weight and weightlessness become increasingly more 
dominant in modern aviation. There can be no doubt 
that the phenomenon weightlessness is about to become 
an outstanding aeromedical problem as aviation pro- 
gresses. So far, only the problems of increased body 
weight were of aeromedical import. Now we have 
reached a point where investigations of external forces 
acting on the body must extend into the range of 
reduced weight down to the complete lack of weight. 

The problems related to increased body weight 
(g-forces greater than 1) concern primarily the inter- 
ference of additional forces with the mechanical factors 
of circulation, respiration, dislocation of organs, body 
movements and their consequences. All these dis- 
turbances are, in one way or another related to elastic 
deformations of the body and its parts. These disturb- 
ances, up to a certain limit, vary approximately 
in proportion to the increase of body weight. Con- 
sequently, the elimination of body weight will probably 
have no decisive or even lethal effects upon the physio- 
logical functions mentioned previously. The transition 
from the condition of normal weight into the state of 
weightlessness removes the normal elastic stress exerted 
upon the body by its support, and there are no major 
indications of a possible failure of the previously enu- 
merated body functions in the absence of this stress. 
The correctness of this kind of reasoning was proved 
by experiments with monkeys and mice in V-2 and 
Aerobee rockets, carried out by Henry and Ballinger 
(22). In several experiments of this kind, test animals 
were carried aloft in research rockets and subjected 
to the state of weightlessness during the free-flight 
period, lasting between 2 to 3 minutes. The rate of 
pulse and breathing, and breathing movements were 
recorded during these tests, and only slight disturb- 
ances were observed. 

In contrast to the functions of purely mechanical 
nature, the mechanical sense organs of the body do 
not react to alterations of the stimulating forces in a 
linear relationship. They follow a logarithmic rela- 
tionship which is described by the Weber-Fechner law. 
Therefore, we must conclude that changes of mechanical 
excitation in the absence of a basic stress will probably 
have pronounced effects on man’s ability to cope with 
the mechanical factors of his environment in the state 
of weightlessness. According to Gauer and Haber (23), 
and Haber and Gerathewohl (24), we must expect 
disturbances in the function of orientation, in the proper 
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execution of body movements, and in the harmony 
of the coordination between the various components of 
man’s mechano-receptor system. The disturbing oc- 
currence of optical illusions, akin to the oculo-gravic 
illusion (33), cannot be excluded. Our empirical knowl- 
edge of the possible physiological and psychophysical 
consequences of weightlessness is rather fragmentary— 
in fact, it is practically nil. This fact constitutes a 
challenge to aeromedical, or better, to space-medical 
research. 


F The Function of Supplying Diffuse Daylight 


The brightness of the daylight sky at a particular 
point of observation depends upon the zenith distance 
of the sun, the difference in azimuth between the sun 
and the point of observation, the mean reflectivity 
of the earth, and—which is of particular interest 
here—the sky brightness depends on altitude. Haber 
(25) has computed the brightness of the daylight sky 
as a function of altitude by using the theory of 
Tousey and Hulburt (26). In Fig. 4 is depicted the 
zenith brightness in millilamberts as a function of 
altitude. As can be seen from the graph, the bright- 
ness of the sky at a height of 20 miles is about '/s0 of its 
value at sea level. Fig. 5 shows the extension of 
zenith brightnesses up to altitudes of 90 miles. Al- 
though the theory of Tousey and Hulbert was not 
designed to stand an extrapolation to such great heights, 
the data given can be expected to be a fairly good 
approximation. According to Fig. 5, the light of the 
sky vanishes completely only at very great altitudes. 
At heights in excess of 85 miles the zenith brightness 
approaches that of a moonless night at sea level. 
From there on the blackness of space is complete. 

The absence of diffuse illumination provided by the 
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atmosphere poses a visual problem. When looking 
down from a high-flying craft, the eye meets bright- 
nesses of varying degrees, depending upon the features 
of the landscape and meteorological conditions. Look- 
ing up, the daylight sky appears as the only source 
of light. Most features of the ground, the clouds and 
other sunlit objects are rather bright. Their bright- 
nesses lie in the range between 300 and 10,000 milli- 
lamberts. In contrast to these values, the bright- 
ness of the sky is rather low even at moderate alti- 
tudes, and decreases further at higher levels. Owing 
to these conditions strong contrasts of simultaneously 
perceived brightnesses are produced. 

The brightness contrast is expressed by the following 
ratio: 


where B, is the brightress of the target, and B, 
the brightness of the background or surround- 
ing. Guth (27) has established contrast levels 
which produce an initial momentary sensation judged 
to be at the borderline between comfort and discomfort 
(BCD). The critical contrast values depend upon 
the brightness of the target, the solid angle (Q) sub- 
tended by the target, and the brightness of the back- 
ground in the following fashion: 


BCD = — 1.28)...... 


whereby the brightnesses are given in footlamberts 
(1 footlambert = 1.0764 millilambert). From the 
data of Guth it can be concluded that viewing of 
objects at higher altitudes will be anything but com- 
fortable. Even at moderate altitudes, objects having 
a reflectance as low as 20 or 30% surpass the threshold 
of comfort, if they are illuminated by the sun when 
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viewed before the background of the sky. For an 
object the size of the wing of an aircraft, the comfort 
threshold is exceeded at a height of 10,000 feet, even 
if the reflectance of the wing surface is only 25%. 
At a height of 20 miles, sunlit objects as small as the 
disk of the moon appear too bright for comfort, if 
their reflectance is larger than 10%. 

Far more discomfortable brightness contrasts will 
be encountered inside of the cabin of a high flying ship 
where sunlit patches will be seen adjacent to deep 
shadows. Owing to the small illumination of the wan- 
ing light of the sky, the interior of a cabin will be rather 
dark. Reading of instruments will be impossible un- 
less they are illuminated artificially. These conditions 
will undoubtedly cause visual discomfort, fatigue, 
or even disabling glare. Cibis (28) has suggested the 
use of light-diffusing panels covering a certain section 
of the window area. In this way, the interior of the 
ship would be shielded against direct sunlight. At 
the same time, such a device would provide a sufficient 
amount of comfortably diffuse illumination from the 
panel that substitutes for the loss of diffuse illumination 
normally afforded by the light of the sky. 


G The Function of Absorbing Meteors 


The hazards arising from possible impacts with 
meteors is of particular space-medical interest. Even 
moderately large meteors can easily puncture the hull 
of a high-flying ship—an event which would result in a 
very rapid explosive decompression of the crew. At 
altitudes accessible to man with present-day means, 
no hazard from meteors exist, since meteors are vapor- 
ized and annihilated at great altitudes. The great 
majority of meteors ranging between fine meteoric 
dust and particles weighing a gram or more, are ab- 
sorbed before reaching the 60-mile level. Above 90 
miles of altitude, however, a ship would run the same 
risk of being hit by meteors as in interplanetary space, 
except for the protection afforded by the shielding bulk 
of the earth. 


H The Function of Interacting Thermally with the 

Craft 

One of the most intriguing problems of modern 
aviation is the generation of an envelope of hot air 
around a fast-flying aircraft owing to aerodynamic 
heating. This is not the place to discuss the many 
aspects of the “heat barrier” which our high-speed 
aircraft encounter. Relative to our subject we are 
primarily interested in the quesion at what altitude 
aerodynamic or friction heating can be expected to 
cease. This question has been answered by Saenger 
(29), who has studied the stability of a satellite vehicle 
at various levels of altitude. He has found that, at a 
height of 80 miles, a satellite would lose over 6 miles 
of altitude per revolution as a consequence of air drag. 
The altitude loss per revolution would be only 3 feet, 
if the satellite would circle at a height of 112 miles. 
Somewhere in the range between 80 and 110 miles of 
of altitude, air drag and, simultaneously, skin heating 


274 


for velocities of the order of 5 miles per second will 
become negligible. 

The temperature problems of rocket flight, however, 
persist even at greater altitudes. Beyond the border 
of thermal interaction between atmosphere and craft, 
the skin temperature will be determined by the ex- 
change of radiation between the exterior of the craft 
on the one side, and the sun, the earth, and the cosmos 
on the other. Buettner (30, 31) has studied this 
problem in detail by calculating equilibrium tempera- 
tures of various metals and other substances. The 
equilibrium temperatures are calculated by equating 
the rate of absorbed radiant energy to the radiation 
losses of the body in space. The latter is proportional 
to the fourth power of temperature (Stefan-Boltsmann 
law) and to the emissivity of the body’s surface in the 
infrared. The surprising result obtained by Buettner 
was the fact that all metals become excessively hot 
under conditions of radiative equilibrium in space, 
even if they are highly polished and reflecting. The 
high temperatures are caused by the poor emissivity 
of the metals in the infrared. A few of Buettner’s 
data are given in Table 2. 


TABLE 2 EQUILIBRIUM TEMPERATURES OF VARIOUS MATERIALS 
IN SPACE, °F 


Material Illuminated Illuminated 
by sun by “full earth” _ 
Black (soot) 252 154 
White(MgO) -60 9 
Aluminum 802 563 


The figures show that only gleaming white materials 
such as magnesium oxide can serve as surface coatings, 
if a manned craft is exposed to the field of thermal 
radiation in space. 

From the standpoint of aviation medicine and space 
medicine the temperature problem has a further im- 
portant aspect. It is the problem of survival in cases of 
failure of cooling systems—in other words: What is 
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the time reserve of an individual more or less suddenly 
exposed to a hot environment? By “heat time re- 
serve” one understands the time span that elapses 
between the beginning of the exposure to heat and the 
moment when heat collapse must be anticipated. 
Naturally, one must expect the heat time reserve to 
be dependent upon many factors such as type of cloth- 
ing, profuseness of sweating, precooling of the subject, 
and the like. In one of his recent papers, Buettner 
(31) has compiled some average data which are given 
in Fig. 6. 


I The Function of Interfering by Air Drag Over Long 


Periods of Time 


: _ Saenger has shown that at an altitude of 120 miles 


air drag has completely vanished. Beginning at this 
height, the establishment of a permanent satellite orbit 
would be possible, although a somewhat greater height 
would be indicated for practical purposes to allow for 
a certain eccentricity of the satellite orbit. Never- 
theless, the 120-mile level can be considered the me- 
chanical border of space. This line is of space medical 
interest insofar as above this level the state of weight- 
lessness can be maintained for any length of time. 
Above an altitude of 120 miles—which, incidentally, 
is only about 3% of the earth’s radius—the environment 
of space is complete. 

Above an altitude of 120 miles there are only three 
factors of terrestrial origin that make the environment 
of the craft different from that found at any other 
point of interplanetary space: (1) The bulk of the 
earth which shields off half the number of meteors and 
cosmic ray particles. (2) The magnetic field of the 
earth which deflects cosmic ray particles below a certain 
magnetic rigidity, if they approach the earth in or near 
the equatorial plane. (3) The solar radiation re- 
flected by the earth and its atmosphere, and the in- 
frared radiation emitted by the earth proper. 

The problems which arise in the operation of manned 
vehicles at very high altitudes and eventually in free 
space are of an extremely diverse and complex nature. 
Their solution requires contributions from meteorology, 
geophysics, astronomy and astrophysics, cosmic ray 
physics, aerodynamics, radiobiology, physiology, avia- 
tion medicine, space medicine, bioclimatology, and hu- 
man engineering. Owing to these many different 
fields, difficulties in semantics are frequently encount- 
ered. Particularly, the term “upper atmosphere” 
is inadequate and misleading, since it conveys different 
meanings in the various fields such as meteorology, 
geophysics, and aviation medicine. For the common 
benefit it appears expedient to coin a new term for 
designating those regions of the atmosphere where, in 
terms of manned rocket flight, the conditions of con- 
ventional aviation blend into those of actual space 
flight. To this end, Buettner (32) proposed the term 
“aeropause.” The aeropause is defined as that re- 
gion of the atmosphere where its various functions 
for man and craft begin to cease and space-equivalent 
conditions are gradually. approached. The concept 
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of the aeropause appears to be quite useful in modern 
aviation and rocketry. It circumscribes the area 
characterized by certain factors of environment that 
are distinctly different from those found in the area of 
conventional aviation or of space. The aeropause 
encompasses the region between the 10- and 120-mile 
levels. The main characteristics of the aeropause 
are depicted in Fig. 7. 
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Recordings in the lower stratosphere have furnished, in controversial issue ever since. This discourse is not 
recent years, direct pictorial evidence for the so-called intended to deal with this problem in its general mean- 
heavy y nuclei component of the primary cosmic radiation. 
ing, but is limited to the question of what possible effects 


This discovery places the problem of a possible exposure 
hazard for humans in the region of the primary radiation 
on a new basis. 

The heavy primaries are characterized by an extremely upper end of and outside the terrestrial atmosphere 
high specific ionization and by a large radial spread of the where the radiation intensity is about 150 times higher 
ionization dosage around the particle tracks. Both quali- . 
ties are apt to endow them with a greatly increased biologi- 
cal effectiveness which cannot be measured adequately in 
terms of ordinary biological dosage units. 

The geomagnetic field exerts strong, deflecting forces 
on the charged primaries. This influence is strongest in 
the equatorial plane and weakest along the polar axes. 


on the human organism might result from an exposure 
in the region of the primary cosmic radiation at the 


than at sea level. 


2 Nature of the Primary Cosmic Radiation 


A revolutionary change in the concepts of the primary 
cosmic radiation has taken place since 1948. In that 


The minimum energy of the primary spectrum, therefore, year, P. Freier and collaborators (1)* could present the 

shows a marked latitude dependence and the exposure first pictorial evidence for the so-called heavy nuclei 

hazard can be expected to be substantially higher in the . odd — 
component of the primary radiation. This new com- 


polar cap than in the equatorial belt. This holds especially 
since the low-energy particles which represent the most 
harmful type occur only at high latitudes. 


ponent has since been studied systematically in several 
investigations (2, 3, 4, 5), and rather accurate informa- 


For interplanetary space a conclusive answer with regard tion on the composition of the primary radiation is now 
to the exposure hazard cannot yet be given since no experi- available 
mental data on the genuinely “‘primary”’ spectrum outside 1h / . 
The old concept (before 1948) that the primary 


any shielding magnetic field are available, and since the 
existence of a heliomagnetic field which would screen off radiation consisted mainly if not exclusively of high 


low-energy primaries from the full region inside the energy protons had to be abandoned on the basis of the 
earth’s orbit is also a doubtful issue. =e new high-altitude recordings. It is now well established 
that the primary radiation consists of atom nuclei of 

many of the elements of the Periodic System, and their 
1 Introduction relative abundance seems rather similar to the abun- 
dances of the elements in the universe as astro-spectros- 
copy and the study of meteorites reveal them. 

The most startling quality of the heavy nuclei is their 
extremely high kinetic energy. It exceeds by far the 
energies which can be obtained from any natural or 
artificial terrestrial radiation source. Due to this ex- 
cessively high energy these nuclei become stripped of all 
their orbital electrons as soon as they penetrate any ab- 
sorber, i.e., they turn into multiply charged ions. This 
means since nuclei as heavy as iron (atomic number 
Z-26) or even tin (Z-50) have been identified in the 
primary radiation, that 26-fold or even 50-fold ionized 
nuclei occur in the heavy component Now, since the 
ionizing power of a high energy particle increases 
with the square of its charge, it is seen that the heavy 
primaries of cosmic radiation must have an extremely 
high energy dissipation along their paths in any ab- 
sorber. It is this quality which represents the specific 
_ hazard of an exposure in the region of the primary 
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ROM outer space a very penetrating ionizing radia- 
tion whose origin is completely unknown is coming 
to the earth. Ever since this radiation was discovered 
and its extraterrestrial origin clearly recognized in 1912, 
cosmic ray research has been presenting atomic physics 
with a never-ending series of startling new phenomena. 
Some of these discoveries, of paramount importance for 
our understanding of the binding forces of the atom 
nucleus, have been made in most recent years, and to 
predict that more such discoveries are still to come 
seems rather trite especially if one considers that, at 
present, complete ignorance must be admitted with re- 
gard to any reasonable explanation of the origin of cos- 
mic radiation. 
Rather early in this development it had been con- 
jectured that cosmic radiation might have some influ- 
ence on living matter, and this problem has been a 
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In talking of the extremely high energies of the cos- 
mic ray primaries one must clearly distinguish between 
the over-all energy of the total incoming cosmic ray flux 
and the energy of the individual particle. It is only the 
latter which exceeds by many ten powers all other radia- 
tions available. The former is very small even in the 
lower stratosphere, where no attenuation due to absorp- 
tion has occurred. That means more precisely, the 
number of incoming particles per unit time and per unit 
cross section is very small. Due to this circumstance, 
cosmic ray research faces the unfortunate double diffi- 
culty that, while experimentation with the primary 
radiation can only be performed in the gondolas of pilot 
balloons or in the warheads of rockets, prolonged ex- 
posure times are required to produce measurable effects. 


3 The Influence of the Geo- and Heliomagnetic 
Fields 


Before proceeding to the presentation of quantitative 
data on the mass and energy spectrum of the primary 
radiation, another phenomenon has to be briefly dis- 
cussed, namely, the influence of the geomagnetic field 
on the primaries. As is well known, the Earth is a mag- 
net and its magnetic dipole field extends far beyond the 
ionosphere into empty space. The cosmic ray pri- 
maries, being positively charged particles, suffer from 
deflection in this field, and it requires a certain minimum 
kinetic energy for a cosmic ray particle to reach the 
Earth at all. Moreover, since the deflecting influence of 
the geomagnetic field is strongest for particles traveling 
in the equatorial plane and weakest along the magnetic 
axis, this minimum energy for reaching the Earth can be 
expected to show a marked latitude dependence. The 
experimental data give indeed a full, quantitative cor- 
roboration of this theoretical conclusion for the lati- 
tude range from zero (equator) to about 55°N. Table 1 
gives the minimum energy necessary for completely 
ionized nuclei to reach the Earth at vertical incidence 
for different latitudes. 


TABLE 1 GEOMAGNETIC CUT-OFF VALUES FOR CHARGED PRIMARIES 
APPROACHING THE EARTH FROM OUTER SPACE* 


——— Geomagnetic latitude——— 


Particle type 55° 41° 30° 0° 
A = Z (Protons) 1.0 BeV 4.0 13.8 


7.3 
A = 2Z (He and heavier) 0.3 1.6 3.3 6.5 


“ The table gives the minimum kinetic energy in billion elec- 
tron-volts required for atom nuclei to reach the earth at vertical 
incidence for different latitudes. A, Atomic weight; Z, Atomic 
number. 


It is seen from these data that, due to the screening 
influences of the geomagnetic field, low-energy particles 
are entirely excluded from reaching the Earth. It re- 
mains under these circumstances an open question 
whether such low-energy particles are present or not in 
the genuinely ‘“‘primary”’ spectrum of cosmic radiation 
in free space far away from any heavenly body with its 
shielding magnetic field. 

Closely connected with this question is the controver- 
sial issue of the existence of a magnetic field of the Sun 
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(7, 8). It has been mentioned previously that the lati- 
tude dependence of the cosmic ray intensity follows the 
prediction of the geomagnetic theory only to about 
55°N. Beyond that latitude a further strong increase 
should be expected theoretically. Actually, however, a 
knee is observed in the curve and the intensity remains 
constant toward higher latitudes, at least, as far as the 
recordings available at present are conclusive. This 
behavior of the intensity curve has been tentatively 
attributed to the influence of the heliomagnetic field 
whose field strength at the Earth’s orbit is assumed, 
according to this hypothesis, to equal in screening 
strength the geomagnetic field at 55°N and take over 
the shielding function from this point to higher latitudes. 

It might be emphasized once more that this hypo- 
thesis is far from being well established. No clear ex- 
perimental decision can be given at present whether the 
Sun has a magnetic field and what its intensity and con- 
figuration in particular might be. Consequently, we 
also do not know whether the absence of the low-energy 
cosmic ray primaries is a secondary effect due to mag- 
netic shielding or an intrinsic quality of the cosmic 
radiation itself. More detailed high-altitude recordings 
close to the geomagnetic poles are needed for the solu- 
tion of this problem. 

For the radiobiological consequences, this uncer- 
tainty is especially unfortunate since the main threat of 
an exposure hazard rests upon the low-energy part of 
the spectrum, as will be discussed in detail later. 


4 The Energy and Mass Spectrum of the 
Primary Radiation 


The full energy spectrum of the primaries is a con- 
tinuum which extends from the geomagnetic cut-off 
energies given above to extremely high values. One of 
the highest energies so far directly measured has been 
found in a primary alpha particle recorded by M. F. 
Kaplon and collaborators (9). Its kinetic energy can 
be analyzed rather accurately because this particle has 
released, in the stack of photographic plates, a nuclear 
disintegration comprising 72 fragments which can be 
followed through many subsequent plates so that the 
total energy balance can be thoroughly computed. 
The kinetic energy of this particle equals 8 xX 10" 
electron-volts, i.e., 8000 billion electron-volts. 

Such rare giants in the primary radiation are the 
starting point of the so-called large cosmic ray showers. 
These are avalanches of many thousands, or even mil- 
lions, of simultaneously occurring particles sometimes 
covering, at sea level, areas of a mile or more in diame- 
ter. The showers gradually develop in cascades of 
successive nuclear disintegrations from such _pri- 
maries of extremely high energy. Though the com- 
plicated mechanism of shower production has been 
thoroughly investigated and very well clarified in its 
different phases, the main question of where the pri- 
mary particles collect their tremendous energy stays 
completely unsolved. 

Besides the energy spectrum, the mass spectrum of 
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the cosmic ray primaries is of interest. Table 2 gives 
the quantitative data on it as they follow from the re- 
cordings. It should be mentioned that one must clearly 
distinguish between the particle number with which the 
different components are represented in the spectrum 
and the energy which they carry. It is a well-estab- 
lished fact, which is so far not yet understood in its 
causation, that the heavy primaries carry equal 
amounts of energy per nucleon; i.e., the average energy 
of a heavy particle is larger than the one of a lighter 
particle by the ratio of the atomic weight. Thus, the 
percentage of the total energy carried by the heavy 
components is considerably larger than the percentage 
of the particle number, as is shown from Table 2. This 
remarkable phenomenon might help some day to shed 
light on the acceleration me chanism and the origin of 
the cosmic ray primaries. 


TABLE 2 MASS SPECTRUM OF THE PRIMARY COSMIC RADIATION 
ON TOP OF THE ATMOSPHERE AT DIFFERENT LATITUDES” 


———Geomagnetic latitude———— 


Particle type 55° 41° 30° 0° 
Protons ; 12,600 4560 2510 1350 
He a 1,820 723 327 163 
C, N,O = 91 38 18 8.8 
Z>10 27 15 6.2 3.1 


° The table gives the number of particles per second traversing 
a horizontal surface of 1 square meter under any angle, i.e., in- 

coming from the full hemisphere. _ 


5 The Transition Effect 


Though the discussion here is limited to the situation 
on top of and outside the atmosphere, it might be per- 
mitted to outline briefly the situation in the lower 
stratosphere and in the troposphere. The high-energy 
primaries when entering the air blanket undergo differ- 
ent types of collision processes in which their energy is 
transmitted to secondary particles. As a consequence, 
the composition of the cosmic radiation in the lower 
regions is completely different from the primary radia- 
tion. At 70,000 ft, practically all primaries are ex- 
tinguished and the secondaries (mesons, neutrons, pho- 
tons, electrons, positrons, and secondary protons and 
alpha particles) are the exclusive carriers of the energy 
flux. In each collision of a primary particle several to 
many such secondaries are generated. Thus, the num- 
ber of particles increases substantially in the initial 
layers of the air blanket until finally at greater depths 
the absorption processes outweigh this multiplication 
and the particle numbers decrease. It is seen that at a 
certain altitude a maximum particle density must occur. 
This happens at about 70,000 ft. 

Parallel to the maximum number of particles the rate 
of energy dissipation per unit mass of absorber passes 
also through a maximum. In radiobiological terms, 
the ionization dosage in living tissue also goes through a 

2 Nucleon is a general term for proton and neutron. Since 
modern theory of the intra-nuclear forces assumes that protons 
and neutrons exchange their identities in the nuclear bond, it 
was useful to introduce a common name for both particles. 
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FIG. 1 IONIZATION DOSAGE IN LIVING TISSUE FROM COSMIC RADI- 
ATION AT NORTHERN LATITUDES 
maximum. Fig. 1 gives the numerical data for this re- 
lationship. 


The phenomenon that the rate of energy dissipation 
of a radiation in an absorber passes through a maximum 
is called transition effect. The transition effect of the 
cosmic radiation plays an important role in the evalua- 
tion of the actual exposure hazard of human beings in 
a ship. The values of the total dosage in Fig. 1 are 
correct only for an infinitesimally small sample of 
living tissue far apart from any additional masses of a 
greater density. They probably still will be correct for 
an object of the size of a mouse. But as soon as a 
larger mass is involved, for instance a rocket ship, 
local transition effects will develop and the ionization 
dosage in the surroundings will be much higher. Fig. 1 
shows that the dosage multiplication factor due to the 
transition effect in the rarefied air beyond 75,000 ft 
equals about 2. It can be expected that for denser 
materials this value will be markedly exceeded. These 
relationships quite generally have to be kept in mind 
as soon as one starts calculating actual figures for the 
ionization dosage in any kind of ship. It would be a 
gross mistake to assume the numerical values of Fig. 1. 
They merely represent the lower limit for very small 


i 6 The Radiobiological Aspects 
“The cru 


cial point in the problem of a possible expo- 
sure hazard at high altitudes, however, is not connected 
to this phenomenon but is of a very different nature. 
As has been seen from the brief analysis of the transition 
processes which occur when the primary radiation inter- 
acts with the air blanket, a continuous change in the 
types of particles which carry the energy flux and pro- 
duce the ionization takes place from the top of the at- 
mosphere to deeper regions. As far as the heavy com- 
ponent of the primary radiation is concerned, Fig. 2 
gives the numerical data for their gradual extinction. 
The primaries transfer their energy mainly to mesons, 
alpha particles, neutrons, protons, and electrons. This 
change in the type of ionizing particle has decisive im- 
portance for the biological effectiveness of the cosmic 
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FIG. 2. DEPTH OF PENETRATION FOR HEAVY NUCLEI OF PRIMARY 
COSMIC RADIATION 


ray energy flux because these secondaries differ greatly 
in specific ionization from the heavy primaries. 

The specific ionization of a radiation is a measure of 
the density of ionization along the path of the individual 
ionizing particle. It is usually given by the number of 
ion pairs per unit path length. The specific ionization 
varies greatly for different types of particles. More- 
over, it varies greatly for the same particle for different 
kinetic energies. A fair first approximation is given by 
the general rule that the specific ionization is propor- 
tional to the square of the charge of the particle and in- 
versely proportional to the square of its velocity. 

It is a well-established fact in radiobiology that the 
specific ionization is a factor which greatly influences the 
quantity of effect of a radiation exposure. The same 
total ionization dosage, say 500 roentgen units, ad- 
ministered once with gamma rays (low specific ioniza- 
tion) and once with alpha rays (high specific ionization), 
will produce a very different degree of radiation damage. 
For the so-called long-term low-dosage damage with 
which we deal in prolonged or often repeated exposures 
to very low intensities, this difference amounts to a 
factor of about 20 for alpha radiation, meaning that 
alpha radiation is twenty times as effective as gamma 
radiation. To take into account this difference in the 
biological effectiveness of the same ionization dosage of 
different types of radiations, the original definition of 
the roentgen unit has been amended by introducing the 
differentiation in two units, the roentgen-equivalent- 
physical and the roentgen-equivalent-man, the former 
giving the ionization dosage, the latter giving the 
biologically effective dosage. 

The specific ionization of the heavy primaries of the 
cosmic radiation covers a wide range from values below 
that for alpha rays to values far above it. Under these 
aspects it is seen that the right-hand part of the curve 
of the total ionization in Fig. 1 represents, with regard 
to the biological effect, something completely different 
from the left-hand part. On the right side of the 
ionization maximum of 70,000 ft, the heavy component 
is participating in the production of the total ioniza- 
tion dosage. The actual biological dosage in this 


region, therefore, can be expected. to be considerably 


characteristic feature of the heavy primaries. 


tracks. Fig. 3 might demonstrate this more clearly. 
It shows a heavy nucleus and an alpha track recorded in 
the same nuclear emulsion plate and microphotographed 
at the same magnification. Though the photo emul- 
sion is not suited to furnish accurate quantitative infor- 
mation, Fig. 3 shows clearly the much larger spread of 
the ionization events in the eross section of the heavy 
nucleus track. The reason for this phenomenon is the 
high kinetic energy of the heavy nucleus which results 
in correspondingly higher kinetic energies transmitted 
to the secondary electrons which the heavy nucleus par- 
ticle produces in collisions with absorber atoms. This 
large radial spread of the ionization dosage around the 
heavy nuclei tracks is a novel phenomenon, since no 
other artificial or natural radiation source furnishes 
these excessively high kinetic energies. It is, therefore, 
not yet possible to appraise the biological significance in 
all its consequences. 

As has been mentioned already, the specific ionization 
varies greatly with the velocity of the particle obeying 
roughly an inverse square law. It is due to this fact 
that the low-energy part of a heavy track represents the 
most harmful section with regard to the exposure haz- 
ard. Fig. 4 demonstrates such a terminal section of a 
heavy track in nuclear emulsion. The left-hand part 
shows the track section of a heavy nucleus of Z 20 at 
15 billion e-volts kinetic energy. The upper part of 
the central section shows the same track at 1.25 billion 
e-volts kinetic energy. The track here is in its terminal 
phase and passes through the so-called ionization peak 
and thin-down part before coming to a final stop. The 
scale at the upper left gives the magnification; 1 divi- 
sion equals 10 microns. It has to be kept in mind, how- 


FIG. 3 NUCLEAR EMULSION TRACKS OF A HEAVY NUCLEUS AND 
AN ORDINARY ALPHA TRACK 


Both pictures were taken from the same plate and at the same 
magnification. Vertical length of the visual field: 82 microns. 
(Recorded by P. Freier and E. P. Ney, University of Minnesota.) 
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higher than in the corresponding region below 70,000 ft. 

The high values of specific ionization are not the only 
In ad- 
dition, profound differences exist with regard to the 
spatial distribution of the ionizations along the particle 
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ever, that Fig. 4 shows the track in nuclear emulsion. 
If it had traversed living tissue, the length of the track, 
because of the considerably lower absorptive power of 
living tissue, would have to be increased by about a 
factor 4. Fig. 4 demonstrates strikingly that the 
ionization peak and thin-down part are the specifically 
harmful phases of the entire track. 

Besides the energy dissipation by ionization due to 
collisions with orbital electrons, heavy nuclei can un- 
dergo another type of collision process, namely, with the 
nuclei of absorber atoms. This process often is called 
catastrophic absorption because, in such a collision, 
usually the total energy of the projectile nucleus is trans- 
formed in one elementary act and a total disruption of 
both the projectile and the target nucleus occurs. 

If such a nuclear collision takes place in an atom of 
heavy atomic weight (silver or bromine) in photo- 
graphic emulsion, the fragments of the nuclear explo- 
sion form a starlike pattern. Microphotographs of 
such disintegration stars in emulsions give an impres- 
sive pictorial demonstration of the large energies in- 
volved. Fig. 5 gives an example of such a high-energy 


NUCLEAR EMULSION TRACK OF A HEAVY NUCLEUS 
OF 2-20 

The left section shows the track at 15 billion e-volts kinetic 
energy. The center and right-hand sections (travel direction 
downward) show the terminal part. The kinetic energy at the 
uppermost point of the center section is 1.25 billion e-volts. 
Seale at upper left: 1 division—10 microns. (Pilot balloon 
recording of P. Freier and E. P. Ney, University of Minne- 
sota.) 


FIG. 4 
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event recorded in the lower stratosphere by H. Yagoda 
(4). One can distinguish in the picture a narrow cone 
consisting of many mesons, a somewhat wider cone of 
similar composition, and several noncollimated heavier 
fragments. 

In a discussion of possible effects on living matter 
one should clearly remember, however, that the special 
emulsions used for such recordings contain more than 
80 per cent in weight of silver bromide, whereas living 
tissue contains heavier atoms only to a negligible per- 
centage. Disintegration stars released in tissue, there- 
fore, can consist only of a small number of prongs. 
Columns of such an extremely dense ionization as is 
seen in the initial part of the narrow meson cone in 
Fig. 5 are not likely to occur in living matter from 
disintegration stars. 

The probability for a nuclear encounter quite gen- 
erally is incomparably smaller than for a collision with 
some orbital electrons resulting in ordinary ionization. 
This ratio is of the order of 10°. That means a heavy 
nucleus before undergoing a catastrophic absorption 
produces on the average about 10° ionizations before 
being annihilated in a nuclear collision. It is seen 
from this relationship, even without entering into 
mathematical detail, that high-energy heavy nuclei 
passing through an absorber have a higher chance to 
be eliminated by catastrophic absorption before reach- 
ing the ionization peak and thin-down phase in the 
terminal part of their track than low-energy nuclei. 
This consequence has important bearing for the ex- 
As has been seen above, the geo- 
rably with lati- 


posure hazard. 
magnetic cut-off values vary conside 


NUCLEAR DISINTEGRATION OF A SILVER OR BROMINE NU- 
CLEUS OF PHOTOGRAPHIC EMULSION RELEASED BY A HIGH ENERGY 


FIG. 5 


PARTICLE OF PRIMARY COSMIC RADIATION 
(Recording of H. Yagoda, Laboratory of Physical Biology; National 


Institute of Health) 
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a The numerical values for this relationship 


happen to be just in the critical range so that, for the 
lowest energy at which heavy nuclei can reach the earth 
in the equatorial belt up to about-30°, the probability 


for catastrophic collision is practically unity, whereas 
p p 


= equatorial belt reach the ionization peak. They are 


at about 50° N this probability has decreased to about 
0.85. That means that practically no particles in the 


all eliminated beforehand by catastrophic absorption. 
At 50° N, however, about 15 per cent of all incoming 
heavy nuclei escape nuclear collisions and pass through 
the ionization peak. The detailed latitude curve for 
this relationship has yet to be investigated more thor- 
oughly, but it seems safe to conclude on the basis of 
the present experimental data that the most harmful 
type of the heavy primaries, those which pass through 
the ionization peak with its excessively high specific 
 jonization, seems to be limited to northern latitudes 
and the polar cap. 
_ Whether this greatly increased hazard due to the low- 
energy heavy primaries is fairly constant in the entire 
region from 55°N to the pole, or whether the heavy com- 
_- ponent at extreme altitudes does not follow the general 
_ leveling-off of the total intensity, remains an open 
- question. It can only be hoped that the modern tools 
_ for high-altitude sounding will be more extensively used 
in the near future at the geomagnetic north pole itself to 
clarify this unknown part of the problem and the con- 
_ troversial issue of the Sun’s magnetic field. 


_ 7 The Exposure Hazard in the Interplanetary 
Region 


Still more vague and uncertain is the situation with 

regard to the interplanetary regions at greater distances 
tf rom the Earth. The highest direct Geiger counter tele- 
scope recordings extend to an altitude at about 150 kilo- 
meters. They indicate that in the range between 40 
_ and 150 kilometers the total intensity of the primary 
~ cosmic radiation is constant. This information, how- 
ever, does not contain any details in regard to the com- 
_ position of the radiation. 

_ For still larger distances from the earth the intensity 


of all components can be expected to increase slowly to 
twice their stratospheric values, due to the fact that the 
influence of the shadow of the earth gradually vanishes. 
For an observer in the upper stratosphere, the Earth 
still covers very nearly the full hemisphere and, there- 
fore, screens off half of the cosmic ray intensity. This 
simple geometric relationship is plotted in Fig. 6. The 
actual intensity can be expected to increase somewhat 
steeper and to higher values than those in the figure 
because at increasing distances from the earth the 
shielding influence of the —— field gradually 


8 Conclusions 

All these conclusions concerning the situation beyond 
the altitude of 150 kilometers which represents at the 
moment the limit of the zone directly accessible for 
measurements are, of course, guesswork. The un- 
known magnitude is the lower part of the energy spec- 
trum of the primary cosmic radiation, and we have seen 
that it is just this part which is most important for the 
problem of a possible exposure hazard. Thus, very 
disappointingly, the present state of knowledge can 
only be summarized in three “‘if’’ statements: 

(1) If there is no such thing as a lower end in the 
genuinely primary spectrum and only the geomagnetic 
field produces artificially the cut-off edges in the spec- 
trum which we actually observe, it can be assumed that 
radiation hazards in regions outside the shielding in- 
fluence of the geomagnetic field will be prohibitively 
great. 

(2) If the hypothesis of the heliomagnetic field 
taking over the shielding function at 50°N latitude 
should prove true, exposure of humans seems feasible— 
though with a rather marked deterioration of the mean 
life expectancy due to the radiation damage—in that 
part of interplanetary space in which the protection due 
to the heliomagnetic field is effective. This would cer- 
tainly be true for the full space inside the Earth’s orbit, 
but it would be rather questionable whether the “‘safe’’ 
zone would extend outward to Mars’ orbit. 

(3) If already in the genuinely primary spectrum 

the low-energy heavy nuclei are absent, ex- 
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posure of humans seems practicable everywhere 
in space with the restriction concerning chroni- 
cal radiation damage mentioned before. What 
the actual figure for the abbreviation of the 
lifetime or the other possible types of radiation 
injury from exposure of humans to low-energy 
heavy nuclei will be cannot yet be sharply de- 
fined. The heavy nuclei component of the 
primary cosmic radiation is a new phenomenon 
not only for the cosmic ray physicist but all the 
more for the radiobiologist. Its exploration is 
for the latter unequally more difficult since 
he has to deal with reactions in living mat- 
ter which, for full development and manifesta- 
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tion, require periods of months and _ large- 


scale statistical experimentation. To perform 


ARS JourRNAL 


Sl 


rr 


t 
n 
t 
ti 
b 
p 
‘ 
N 
I 
n 
= I 
| 
ly 
FIG. 6 IONIZATION DOSAGE IN LI) 
= 


the 
es, 


such experimentation in the lower stratosphere is both 
expensive and time-consuming. 
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Jet Propulsion News 


C. F. WARNER, Purdue University, Associate Editor — 
with the assistance of W. G. BOHL 


Rockets and Missiles 


: a A ‘CORDING to an article by G. L. Christian in 
Aviation Week, a scale model of one of the Navy’s 
tailless delta-wing aircraft is being flight-tested by 
means of rocket power. The model, equipped with 
i telemetering apparatus, is attached near the tail of a 
standard 5-in. military rocket. The extensive tele- 
metering equipment built by Bendix Aviation Corpora- 
tion records continuously such in-flight information as 
ram air pressure, longitudinal acceleration and decelera- 
tion, yaw and pitch acceleration, and amplitude and 
frequency of wing vibrations. The use of the relatively 
inexpensive rocket-powered model eliminates the need 
for wind tunnel tests. The model, approximately 12 ft 
7 long with a span of 4 ft, reaches a flight Mach number of 
over 1 in approximately two seconds and has a flight 
time of twenty seconds. 
+ + + 


THE Navy has announced that its Douglas Sky- 
rocket established a new world altitude record for 
piloted aircraft of higher than 79,000 ft during its June 
15, 1951, flight. The flight was checked and certified by 
the National Advisory Committee for Aeronautics. 
A speed of 1238 mph has also been confirmed. Sa — 


IT has been reported that the Borg-Warner Corpora- 
tion has obtained a patent for the design of a rocket- 
ramjet missile motor. The use of nitroaliphatic rocket 
propellants by supplying additional oxygen are said to 
contribute to a more stable flame front. The patent 
claims that the use of these oxidizers increases the 
thrust and efficiency of the ramjet. 


+ + 


THE Northrop Decelerator, designed to ‘‘determine 
the tolerances of the human body’s resistance to deceler- 
ative forces so that better equipment can be provided 
for emergency and crash landings,’’ employs from one to 
four Aerojet 5 AS-1000 solid rockets, each rated at 
5000 Ib-sec of impulse, to reach its maximum speed. 
The Decelerator consists of a sled riding on a 2000-ft 
long track. At 150 mph, brakes are applied which slow 
the sled to 70 mph in a distance of 25 ft in a mere 0.16 
sec. Both dummies and human volunteers have been 
used in the tests. The decelerative force is applied in 
a 50-ft section of the 2000-ft track, where 45 brakes, 
each with four brake shoes, are remotely applied by in- 
struments developed especially for use with the sled. 


FIG. | AEROJET AEROBEE SOUNDING ROCKET LEAVING THE 
LAUNCHING TOWER UNDER TREMENDOUS THRUST FROM ITS SOLID 
PROPELLANT BOOSTER ROCKET 


Much of the equipment tested has already proved its 
usefulness in saving lives. 
+ + 


A CONTRACT for production of missile motor com- 
ponents of a type recently developed by Aerojet Engi- 
neering Corporation has been awarded by Aerojet to 
Ryan Aeronautical Company. This is the second Aerojet 
rocket project for which Ryan has been selected as a pro- 
ductive source. Previously Ryan built major compo- 
nents of the Aerobee high-altitude sounding rocket. 
Currently, Ryan is also building complete rocket motors 
for surface-to-surfaee missiles under a production con- 
tract from Firestone Tire and Rubber Company. Ryan 
also supplied this rocket motor under previous experi- 
mental contracts to Douglas Aircraft. The motor was 
developed by the Caltech Jet Propulsion Laboratory. 


* The information reported in this Section has been selected from approved news releases originating with the Department of De- 
fense, private manufacturers, universities, etc., and from published news accounts in journals and newspapers. The reports are con- 


sidered generally reliable, although no attempt has been made to verify them. 
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IT has been reported that, under Britain’s large un- 
publicized guided missile program, missiles are being 
built by six major aircraft or engine manufacturers. 
Several Australian-built target drones powered by 
British Armstrong Siddeley Adder turbojets have been 
flown successfully. Several models of ramjets and ex- 
pendable engines are to be tested soon at the Australian 
Woomera range. Since test flights are restricted by 
short ranges and weather in England, a joint British- 
Australian proving ground was started in Australia in 
1947. This area, known as the Long Range Weapons 
Establishment, extends 1200 miles into the open desert. 

+ 

McDONNELL Aircraft Company, in a report on 17 
airplane, helicopter, and guided missile models, revealed 
that it has been manufacturing for the U. 8. Navy 
a self-propelled armor-piercing radio-controlled pilotless 
bomb designated the RTV-2 Gargoyle. 

+ + 


A NEW altitude record of 80 miles was recently estab- 
lished for the Aerojet AEROBEE free-flight sounding 
rocket in launchings at the proving grounds in New 
Mexico. By making it possible to investigate the 
upper atmosphere for conditions of temperature, com- 
position, density, magnetic-field, cosmic ray intensity, 
x-ray intensity, winds aloft, etc., the AEROBEE has 
been the medium for obtaining valuable scientific data 
heretofore unavailable. Since initial launchings in 1947, 
it has been the most widely used and least expensive 
research sounding rocket available for transporting 
scientific instrumentation on nearly vertical flight. 

Carrying a nominal pay load of 150 lb, the AEROBEE 
vehicle is powered by a bipropellant liquid rocket and 
is initially accelerated by a large JATO booster rocket, 
which detaches and falls away upon cessation of booster 
thrust (Fig. 1). After reaching a high velocity during 
powered flight, the AEROBEE coasts to the zenith alti- 
tude, at which time nose ejection may be accomplished 
for parachute recovery of the nose section. Launched 
from a tower 150 ft in height, the AEROBEE is 
stabilized with fixed fins. The vehicle is 20 ft long 
and 15 in. in diam, the booster rocket adding another 
6 ft to the over-all length. 


Aircraft 


THE USAF has announced the existence of three new 
jet fighters, the North American XF-100, the McDon- 
nell XF-101 Voodoo, and the Convair XF-102, and a 
light bomber, the Douglas RB-66. 

+ + + 


AMERICA’S latest answer to air invasion, an almost 
automatic jet warplane armed with radar and rockets to 
outshoot enemy bombers, is now in production at Lock- 
heed Aircraft Corporation. The new airplane, the 
Starfire F-94C, built for USAF, carries twenty-four 
2.75-in. rockets housed in a ring of firing tubes around 
the nose. Different views of this fighter-intercepter are 
shown in Figs. 2, 3, and 4. as 7 
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FIG. 2.) THE ALL-ROCKET-ARMED LOCKHEED F-94C STARFIRE 


Radar and an electrical ‘‘copilot’’ enable the Starfire 
to spot the enemy miles away, lock onto the target, 
track, close, aim, and open fire without pilot assistance. 
The main duties of the pilot and radar operator are to 
take the plane off, maneuver to the general target area 
as guided by ground radar watchers, switch on the 
“electronic crew” at the proper time, monitor operation 
of the piloting and rocket-control apparatus during 
attack, and then land. The automatic pilot with un- 
limited maneuverability was developed by the Westing- 
house Electric Corporation. With the autopilot at the 
controls, the plane will be able to take the steepest dives, 
turns, rolls, and loops with ease. 

The Starfire, with its 1200 lb of electronic gear, is 
powered by a Pratt and Whitney J-48-P-5, 6250-lb 
thrust jet engine equipped with an afterburner. A rib- 
bon parachute, stowed in a tail compartment, can be 
released just as the ship lands, to halt it in a short space. 
The take-off weight is more than 20,000 Ib; it has a 
length of 41 ft and a span of 37'/. feet. 

+ + + 


A RADICALLY new flying-boat hull bottom, affixed 
to the wing and full crown of a World War II vintage 
seaplane, is currently being taxi- and flight-tested on 
and over the waters of the Chesapeake Bay. Desig- 
nated as research model M-270, the hull embodies in 
full scale the result of years of extensive towing tank and 
wind tunnel tests, according to a joint announcement by 


FIG. 3. THE F-94C WITH ROCKET DOORS OPEN 
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Fic. 4 THE STARFIRE IN FLIGHT 


the U. S. Navy and the Glenn L. Martin Company. 

For years, flying boat engineers have envisioned larger 
slimer hulls to provide superior aerodynamic seaplane 
shapes as well as more efficient structural bodies. These 
shapes are likewise dictated through the necessity of 
having more efficient hulls in which to use jet engines, as 
these will give seaplanes of the future much greater 
speed. Fig. 5 shows a model of the M-270 in the towing 
tank at Stevens Institute of Technology. 

+ + + 


ACCORDING to Jnteravia, the Russians are flight- 
testing a two-jet swept-wing bomber designated the 
Type 150. The new bomber, designed by a group of 
German engineers under the leadership of Dr. B. Baade, 
formerly of the Junkers works, is said to have a span of 
102 ft, a length of 82 ft, and to be powered by two 
11,000-lb thrust Lulkov axial-flow turbojet engines. 
The bomber is reported to have a speed of 620 mph. _ 


+ + + 


IT is reported that a French process has been devel- 
oped for building supersonic aircraft wings from pre- 
stressed concrete. An experimentally developed wing 
by Bregault is cheaper and easier to build than a con- 
ventional wing. The new wing is also reported to have 
greater strength. 


Turbojet Engines 

+ 

BRISTOL has developed a twin-spool compressor jet 
engine—the Olympus—in the 10,000-lb thrust class. It 
is claimed to be the world’s most economical and most 
powerful jet engine at this stage of development. The 
engine is able to accelerate from idling to full power in 
two seconds and is said to have a fuel rate of 0.766 lb/Ib 
thrust/hr. The Olympus has a frontal diameter of only 
40 in., a length of 124 in., and weighs 3250 lb. Under 
its agreement with Bristol, Wright Aeronautical has had 
an Olympus under test at Wood-Ridge, N. J. Wright 
Aero has been able to obtain a static thrust of 9800 Ib 
from the engine. 


A NEW wind tunnel has been placed in operation at 
the Polytechnic Institute of Brooklyn capable of attain- 
ing Mach 5 in its 4-in. X 4-in. test section. This blow- 
down tunnel, designed by Dr. Paul A. Libby and Dr. 
Antonio Ferri, is supplied with 100 cu ft of air at 3000 
psia from Navy-loaned torpedo air bottles. Other 
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nozzles are available for operating the tunnel at Mach 2 
and 3. Schlieren and shadowgraph apparatus may be 
used for photographic studies. 

~ + + 


A TRANSONIC wind tunnel at the Navy’s David 
Taylor Model Basin, when completed late in 1953, will 
permit testing of models at speeds up to 900 mph. The 
tunnel, based on a design of a tunnel partly built by the 
Germans near Munich, will be of closed circuit type, 574 
ft long. Construction will be of reinforced concrete, 
with a steel test section. 

+ + 


THE world’s largest and most powerful “human 
centrifuge,” a machine to subject pilots to extreme 
gravitational conditions encountered in sonic-speed air- 
craft, has been put into operation at the Naval Aviation 
Medical Acceleration Laboratory. Built by the Mc- 
Kiernan-Terry Corporation, the centrifuge is capable of 
accelerating from a dead stop to 173 mph in less than 
seven seconds, and will speed up from zero to approxi- 
mately 90 mph in one and one half seconds. 

The giant ‘‘whirligig”’ is capable of exerting a force 
equal to 40 times the gravitational pull of the earth on a 
pilot seated in a gondola at the end of a 50-ft rotating 
arm. It will be used by the Navy to evaluate human 
tolerance to acceleration and to study the physiological 
systems which limit the tolerance. Since the human 
“guinea pig’ is visible to researchers only when the 
centrifuge is motionless, high-speed television, x-ray, 
and motion picture cameras are mounted in the gondola. 


Personalities 


W. L. ROGERS, who for the past year has been man- 
aging Aerojet’s Washington office, has been appointed 
assistant general manager. The announcement is made 
by W. E. ZISCH, General Manager. 


+ + + 


RAYMOND W. YOUNG, president of Reaction 


Motors, Inc., of Rockaway, N. J., has announced the 


FIG. 5 MODEL OF M-270 IN TOWING TANK AT STEVENS INSTITUTE 
OF TECHNOLOGY 
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appointment of DR. LUIGI CROCCO, R. H. Goddard 
Professor at the Daniel and Florence Guggenheim Jet 
Propulsion Center, Princeton University, to the Advi- 
sory Board of the Company. In accepting the appoint- 
ment to the Advisory Board of Reaction Motors, Dr. 
Crocco is joining two other outstanding scientists, 
DR. HUGH 8. TAYLOR of Princeton University, and 
DR. WILLIAM R. SEARS of Cornell University. 
+ + + 


ROY E. MARQUARDT, president and chairman of 
the Board of Marquardt Aircraft Co., has announced 
the election of two new directors at the annual stock- 
holders meeting held today at the Company’s Van Nuys 
headquarters. Former board members were re-elected. 

The new directors are DR. CLARK B. MILLIKAN 
of California Institute of Technology, and CLARENCE 
E. UNTERBERG of C. E. Unterberg & Co.,a New York 
securities underwriting concern. Dr. Millikan, long 
associated with the aircraft industry, is chairman of the 
Guided Missiles Committee of the Research and Devel- 
opment Board, and is director of Southern California 
Cooperative Wind Tunnel and of the Guggenheim 
Aeronautical Laboratory at the California Institute of 
Technology. He is also a member of the Air Force 
Scientific Advisory Board Committee and chairman of 
the Subcommittee on Fluid Mechanics of the National 
Advisory Committee for Aeronautics. Mr. Unterberg 
founded the New York investment banking firm which 
bears his name. 


In addition to Marquardt’s development and pro- 
duction contracts with the U. 8. Armed Services for sub- 
sonic and supersonic ramjet engines for many of their 
guided missiles projects, the Marquardt Aircraft Com- 
pany is producing turbine-driven auxiliary power units, 
fuel and power controls for missile and other applica- 
tions. The Company also has engineering and manu- 
facturing subcontracts with such companies as Boeing, 
Douglas, Grumman, Westinghouse, and others. 


[sothermal Combustion Under Flow Conditions 


(Continued from page 251) 


However, the performance of the isothermal expan- 
sion cycle can exceed that of an adiabatic expansion 
cycle in the event that thermal limitations of structural 
materials preclude the adiabatic utilization of the total 


energy available to the working fluid. 
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ey New Mexico-West Texas 


Section Elects 


CAUSE of increased interest in 

ARS activities among the personnel 
of the various organizations using the 
White Sands Proving Ground testing 
facilities, the New Mexico-West Texas 
Section is concentrating on an enlarged 
membership and has elected new officers to 
serve until December 13, 1952, or until 
1953, the tenure of office to be decided 
later. 

The election was completed on May 23, 
1952, and the following officers were 
chosen: President, Edward T. Munnell; 
Vice-President, Robert B. Bolles; Secre- 
tary-Treasurer, Edward E. Francisco, Jr. 

The Section also held a business meeting 
on July 31 to elect a board of directors and 
appoint several committees. The results 
of this meeting will be given in the next 
issue of the JouRNAL. 


A. C. Clarke Speaks at Open Meeting 


On July 1, 1952, the New Mexico-West 
Texas Section held an open meeting at the 
Branigan Library, Las Cruces, N. M. 
Guests and members attended to hear 
Arthur C Clarke, chairman of the British 
Interplanetary Society and author of “The 
Exploration of Space,” which was recently 
published by Harper & Bros. and chosen 
for the Book-of-the-Month Club July 
selection. 

Mr. Clarke discussed various aspects of 
interplanetary flight, pointing out the 
various difficulties—engineering, medical, 
and economic—that must be overcome. 
He concluded that each of these problems 
was susceptible of eventual solution, and 
that extraterrestrial flight would some 
day be a reality. 

The lecture was followed by a question- 
and-answer period in which a number of 
technical problems were discussed at con- 


siderable length. 


Alabama Section Hears 
Discussion on Rocket and 


Guided Missile Materials 


HE ARS Alabama Section held its 
monthly meeting on July 8, 1952, in 
the auditorium of the Huntsville Electric 
System, Huntsville, Ala. Hans Heuter, 


vice-president of the Section, presided in 
the absence of the president, Martin B. 
Schilling. 

Following a brief business meeting, a 


New Officers a” 


discussion of “Materials in the Develop- 
ment of Rockets and Guided Missiles” was 
presented by the guest speaker, W. H 
Steurer, chief of the Material Test Lab 
Station, Guided Missile Development 
Group, Redstone Arsenal. 

Dr. Steurer outlined the difference be- 
tween rocket engineering and general en- 
gineering practices with regard to material 
requirements. He pointed out that the 
primary consideration in the former in- 
stance is development of materials to 
withstand wide extremes of temperatures 
for brief periods of time, whereas the lat- 
ter case involves development of materials 
for application in nominal temperature 
ranges over long periods of time. Many 
properties of materials which heretofore 
have not been investigated, he said, are 
now being examined to determine new as- 
pects of their high temperature behavior 
and their basic metallurgical background. 

The speaker continued with an analysis 
of the properties required of metallic and 
nonmetallic materials during actual flight 
conditions, such as high temperatures and 
various atmospheric conditions, and the 
methods of testing materials to determine 
their acceptability. He compared prop- 
erties of various materials to illustrate the 
criteria by which materials are selected for 
application in rockets and guided missiles 
with regard to fabrication and performance 
in flight. Some insight was given regard- 
ing the development of new materials, 
such as titanium and titanium alloys, and 
their acceptance in rocket engineering. 

Dr. Steurer concluded by discussing the 
need for and methods of evaluating missile 
remnants after destruction (i.e., after 
firing) by metallurgical investigation to 


reconstruct events which occurred during 
flight. He illustrated his talk with movies 
and slides of test procedures, 


Arthur C. Clarke Guest 
Lecturer at Washington- 
Baltimore Section Meeting 


ANOTHE R successful lecture meeting 

of the ARS Washington-Baltimore 
Section was held on June 6, 1952, in the 
auditorium of the Baltimore Polytechnic 
Institute, Baltimore, Md. Approximately 
250 members and guests were present to 
hear the British scientist and chairman of 
the British Interplanetary Society, Arthur 
C. Clarke, deliver a talk on his recently 
published book, ‘The Exploration of 
Space.” 

The meeting was called to order by Joel 
M. Jacobson of Aircraft Armaments, Inc., 
who substituted in the absence of the 
president, Harry J. Archer, Jr. Mr. 
Jacobson then introduced Andrew G. 
Haley, ARS general counsel, who spoke a 
few words about Mr. Clarke and his book 
before presenting the guest speaker to the 
audience. 

Mr. Clarke was born at Minehead, 
Somerset, England, in 1917. Following 
his secondary education and a subsequent 
job as a government auditor, he entered 
the Royal Air Force in 1941. As flight 
lieutenant, he worked with a Blind Land- 
ing Unit and also lectured and wrote on 
radar, physics, and electronics. After 
World War II Mr. Clarke studied at King’s 
College in London, where he obtained his 
B.Sc. in physics and mathematics. Since 


ARTHUR C. CLARKE (SECOND FROM RIGHT), LECTURER, WITH ARS MEMBERS AT THE WASH- 
INGTON-BALTIMORE SECTION MEETING, JUNE 6, 1952 


(Left to right: 
counsel; Mr. Clarke; 


William A. Webb, Baltimore vice-president; 
and Miss Virginia R. Erwin, 


Andrew G. Haley, ARS general 
secretary of the Washington-Baltimore 


Section.) 
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DISCUSSION PANEL AT THE WASHINGTON-BALTIMORE SECTION MEETING, JI 


6, 1952 


(Left to right: Col. H. J. Halberstadt, Deputy Director of Aeronautics and Propulsion, ARDC; 

R. W. Bonner, Head of Propulsion, Glenn L. Martin Co.; G. L. Bryan, Jr., Professor of Aero- 

nautical Engineering, The Johns Hopkins University; and John Creutz, consulting communi- 
cation engineer.) 


then he has been a full-time science writer 
and broadcaster and an energetic pro- 
ponent of space flight, astronautics, and 
astronomy. 

Following Mr. Clarke’s speech, a panel 
type discussion of the lecture was con- 
ducted by the Baltimore vice-president, 
William A. Webb of Aircraft Armaments, 
Inc. Members of the panel included Col. 
H. J. Halberstadt, deputy director of aero- 
nauties and propulsion for the Air Re- 
search and Development Command; R. 


W. Bonner, head of propulsion for the 
Glenn L. Martin Co.; G. L. Bryan, Jr., 
professor of aeronautical engineering at 
The Johns Hopkins University; and John 
Creutz, consulting communications en- 
gineer, 

This meeting was the last of the spring 
meetings of the Washington-Baltimore 
Section. The next meeting is tentatively 
scheduled to be held in Washington in 
September, at which time an active mem- 
bership drive will be undertaken. 


MODEL 408 
RECORDING OSCILLOGRAPH 


FOR VIBRATION, TEMPERATURE, STRESS, STRAIN ANALYSIS 


After conducting a na- 
tion-wide survey among 
the users of oscillograph 
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them for a number of 
years, Century is proud 
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Technical Literature Digest 


Book Reviews 


AsTRONOMY OF STELLAR ENERGY AND 
Decay, by Martin Johnson, London, 


New York, 216 pp. $3.50. - 
Reviewed by A. G. WILSON 


Dr. Martin Johnson, already noted for 
his lucid scientific and philosophical 
expositions in his books, ‘Time, Know- 
ledge, and the Nebulae” and “Science and 
the Meaning of Truth,” engages in this 
latest book in the description of the physi- 
cal problems germane to the stability 
and evolution of stars. The author has 
adopted a novel approach to an exposi- 
tory work. He has divided the book into 
two parts. The first half consists of a 
descriptive background at a popular 
level designed to familiarize the non- 
technical reader with the significant as- 
pects of the subject of stellar constitu- 
tion; the second half amplifies the de- 
_ seriptive material of the first half by 
outlining, with mathematical aid, the 


methods used by the astrophysicist to 
‘ arrive at his results. Such subjects as 


stellar models, stability criteria, and the 
size, masses, temperatures, and energy 
sources of stars are discussed. 

The book serves as an excellent intro- 
duction to the subjects. But the attempt 
to compress the mathematical elaboration 
of the descriptive material of the first 
7 part into an equal number of pages in the 
second part necessarily results in a sketchi- 

- ness which might prove frustrating to a 
beginning student. The book is further 
_ marred by a few careless statements, such 
_ _as the remark on page 83 that “there are 
7 _ only one or two telescopes in the world 

_ of sufficient power to detect a Super- 
Nova.” 


FounpaTions oF HicH SpEEpD AERO- 
DYNAMICS. Facsimiles of Nineteen Ori- 
ginal Papers, with a Bibiography, by 
G. F. Carrier, Dover Publications, Inc., 
New York, 1951, 286 pp. $1.75. 

REVIEWED By H. S. Tsien 
California Institute of Technology 


7 Out of the immense volume of literature 
on the dynamics of compressible fluid 
7 _— from 1870 to 1950, Professor Carrier 


tions of many diverse fields of knowledge. 
290 


Faber, 1950; Dover Publications, Inc., 


literature is unavailable because of security restrictions. 


selected nineteen papers which he con- 
sidered to be original and fundamental. 
The task was certainly not easy. Hence 
he should be excused for the somewhat 
arbitrary manner in making the selections. 
However, for a price of $1.75, students 
of compressible flow can be grateful to 
Professor Carrier that they can read 
the masterful works of W. Rankine, 
G. I. Taylor, Lord Rayleigh, J. Ackeret, 
H. Glauert, L. Prandtl, A. Busemann, 
and Th. von K4rmén. 


HanpBook OF SUPERSONIC AERO- 
DYNAMICS, VOLUME IV, NAVORD 
Report 1488, compiled and edited by 
C. N. Warfield, L. L. Cronvich, A. R. 
Eaton, Jr., G. M. Edelmann, and F. 
K. Hill, U. S. Govt. Printing Office, 
Washinton, D. C., 1952, 121 pp. 
$1.25 


REVIEWED By H. S. 


California Institute of Technology 


The main body of this present volume 
of “Handbook of Supersonic Aero- 
dynamics” consists of tabulated lift 
coefficients and moment coefficients for a 
two-dimensional, flat-plate airfoil in ver- 
tical and angular oscillations, computed 
under the supervision of E. C. Kennedy. 
The accuracy of these tables is very 
much higher than that of the previously 
published tables based upon the Curtiss- 
Wright Corporation Report No. P537V-28 
(April 1948). The range of independent 
variables is: frequency parameter Q, 0.01 
(0.01) 0.04 (0.02) 0.10 (0.05) 0.40 (0.10) 
1.00 (0.20) 3.00 (0.50) 5.00 (2.50) 10.00 
(5.00) 20.00; Mach number M, 1.1 (0.1) 
2.0 (0.2) 4.0 (0.5) 5.0 (1.0) 12.0. 

A brief review of the various flutter 
problems is included in the introduction. 


Out OF THE Sky, by H. H. Nininger, Univ. 
of Denver Press, 1952, viii + 336 pp., 52 
plates and 21 text figures. $5. 


Reviewed by JouHn D. BuppHUE 
Institute of Meteoritics 
University of New Mexico 


The author is the founder and director 
of the American Meteorite Museum which 
contains one of the finest private meteorite 
collections in this country, if not in the 
world. He is a recognized authority on 
meteoritics, especially as a field worker 
and by his efforts many meteorites have 


been brought to light that otherwise might 
have remained forever unknown and lost. 
He is also the author of two other books 
and more than 120 papers on meteorites. 

The book contains 33 chapters devoted 
to such subjects as Early Man and Meteo- 
rites, Typical Falls, Meteors, Composition 
and Structure of the Various Kinds of 
Meteorites, Etching and Preservation, 
Size of Meteorites, Weathering of Meteo- 
rites, Shapes and Shaping, Meteorite 
Showers, Meteorite Craters, Meteoritic 
Clouds, Origin of Meteorites, Cosmic Heat- 
ing, Planetoidal Encounters, Meteorites 
and the Moon, Tektites, and Great Meteo- 
rite Collections. 

The book appears to have been directed 
primarily to the intelligent layman al- 
though it contains a great deal of interest 
to the student, the collector, the profes- 
sional meteorist, and other persons who 
need a semitechnical knowledge of these 
objects. Some of its contents do not seem 
to have been published previously. Prac- 
tically every phase of meteoritics is treated 
in more or less detail, although the re- 
viewer would have preferred more analyt- 
ical data in the chapters dealing with 
composition, especially the mineral com- 
position. A more precise elucidation of 
the origin of the octahedral structure of 
most iron meteorites and how it is in- 
fluenced by the nickel content would prob- 
ably have been a useful addition. The 
chapter on weathering also lacks analytical 
data, but these omissions are doubtless 
intentional and due to the semipopular na- 
ture of the work. 

Two chapters are devoted to the 
author’s hypothesis that the Canyon Dia- 
blo meteorite crater was formed by the 
successive impacts of two meteorites, one 
following closely behind the _ other. 
Whether or not one is disposed to accept 
this hypothesis, the chapters contain a 
wealth of new, and relatively new, infor- 
mation about this crater. 

Other meteorite craters are treated in 
considerably less detail, but cryptovolcanic 
structures, as possible meteorite craters, 
are not mentioned until a later chapter. 
Although he does not say so, the author 


‘seems to favor the hypothesis that the 


Carolina bays are of meteoritic origin. 

One chapter describes the surface fea- 
tures of the moon and suggests that the 
craters are in fact meteorite scars. In an- 
other it is proposed that tektites are bits of 


Eprtor’s Note: This collection of references is not intended to be comprehensive, but is rather a selection of the most significant and 
stimulating papers which have come to the attention of the contributors. The readers will understand that a considerable body of 
ql We invite contributions to this department of references which have not come 
_ toour attention, as well as comment on how the department may better serve its function of providing leads to the jet propulsion applica- 
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fused lunar material blasted beyond its 
gravitational control by the impact of 
meteorites. The bulk of this chapter, how- 
ever, is a résumé of the theories proposed 
by various workers to explain these mys- 
terious objects. In the reviewer’s opinion, 
it would have been well to include a fuller 
description of their surface markings and 
the ways in which the different kinds differ 
from one another, including representative 
analyses, 

The book ends with a list of all the 
author’s books and papers on meteoritics, 
followed by a 13-page index. 

As a whole, the book is lucidly and in- 
terestingly written and fairly well docu- 
mented. It contains a great deal of infor- 
mation and, while not all meteorists will 
agree with it in its entirety, it contains few 
real errors, typographical or otherwise. 
One might mention, however, that osborn- 
ite has been shown to be titanium nitride 
(not an oxysulfide of calcium and either 
titanium or zirconium), and a theory ex- 
plaining the lamellar nature of meteorite 
rust, which the author generously credits 
to the reviewer, is actually due to D. M. 
Barringer, Sr. This book should be a 
valuable addition to the library of anyone 
interested in meteoritics or related fields 
of endeavor. 


PUNCHED Carps, edited by W. A. Casey 
and J. W. Perry, Reinhold Publishing 


Corp., New York, 1951, viii + 506 pp. . 


$10. 


Reviewed by DENsLow 
Jet Propulsion Laboratory 
California Institute of Technology 


The greater part of this book is directed 
to the application of punched-card sys- 
tems as a means of organizing technical in- 
formation. There are also sections on ap- 
plications to industrial engineering and to 
scientific computations. The book com- 
prises thirty sections, each an independent 
article by an authority on the subject. 
Short biographies of the authors precede 
the articles. 

Putting first things first, the reader is 
inspired by case histories of applications of 
punched cards to specific problems in order 
to envision their application to his own 
problem requiring organization of facts for 
interpretation. Then, assuming no prior 
knowledge, the fundamentals necessary for 
formulation of a punched-card system are 
given. Basic theory, preferred practices in 
coding, and mechanical methods are dis- 
cussed fully. An excellent bibliography is 
included. Of particular interest, an article 
by Ernst Hermann Erich Pietsch is in- 
cluded, entitled “Future Possibilities of 
Applying Mechanized Methods to Scien- 
tific and Technical Literature.” A short 
history and a discussion of problems in 
organization of the ‘Gmelin Handbuch” 
are contained, and the peculiar problems 
of the organization of philosophical litera- 
ture are discussed. 
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Mope. Jets AND Rockets For Boys, by 
Raymond F. Yates, Harper and Bros., 
1952, 108 pp. + viii, 45 illustrations. 
$2.50. 


2 Reviewed by R. C. TersBeck 
Jet Propulsion Laboratory 
California Institute of Technology 


This book is written to familiarize young 
people with the operation of rockets and 
jets and their history. The author has 
presented an authentic historical review 
of the subject and has outlined the prog- 
ress in the field during the past years. 

The fundamental principles are de- 
scribed in such a manner that the average 
youth can understand and utilize the in- 
formation to design and build many types 
of rocket and jet motors. 

The book describes in detail the con- 
struction and operation of mode] steam 
jets, a CO, cartridge jet, and simple ve- 
hicles to be driven by these motors. Many 
of the model jet and rocket motors such as 
the Dyna Jet, Jetex (solid propellant 
motor made in England), and Minijet are 
discussed together with installation and 
operation instructions for model installa- 
tions. 

The model rocket and jet enthusiast 
will find this a good reference book for 
answers to many of the operational prob- 
lems that he encounters with his hobby. 


Books 


Histoire de la Mécanique, by R. Dugan, 
Neuchatel, Switzerland, Griffon, 1950, 
649 pp. $17.50. 

An Introduction to Acoustics, by R. H. 
Randall, Addison-Wesley Press, Inc., 
Cambridge, Mass., 1951, xii + 340 pp. 
$6. 

Mechanics of Vibration by H. M. 
Hansen and P. F. Chenea, John Wiley & 
Sons, Ine., New York, 1952, 417 pp. 

Internal Combustion Engines, 6th Edi- 
tion, by L. C. Lichty, McGraw-Hill Book 
Co., Ine., New York, 1951, 598 pp. $7. 

Creep of Metals, by L. A. Rotherham, 
Inst. of Phys., London, 1951, 80 pp. + 2 
plates. 15s. 

Diffusion in Solids, Liquids, Gases, by 
W. Jost, Academic Press, Inc., New York, 
1952, 558 pp. $11. 

Rockets, Jets, Guided Missiles and 
Space Ships, by J. Coggins and F. Pratt, 
Random House, Inc., New York, 1951, 58 
pp. $1. 

Is Another World Watching?, by G. 
Heard, Harper and Bros., New York, 
1951, 183 pp. $2.35. 


a 


Jet Propulsion Engines 


Some Fundamental Considerations on 
the Propulsion of Aircraft (in Dutch), by 
J. G. Slotboom, Ingenieur, vol. 63, Oct. 
1951, pp. L.49-L.56. 


Rocket Propulsion _ 
Engines 


Liquid Films Formed by Means of 
Rotating Disks, by B. E. Dixon, A. A. W. 
Russell, and J. E. L. Swallow, Brit. J. 
Appl. Phys., vol. 3, Apr. 1952, pp. 115-119. 

The Atomization of Liquid in a Flat 
Spray, by R. G. Dorman, Brit. J. Appl. 
Phys., vol. 3, June 1952, pp. 189-192. 

JATO and Auxiliary Rocket Power 
Plants, by H. Gartmann, Weltraumfahrt, 
no. 6, 1951, pp. 134-139. 

Jet Assisted Take-Off, by W. J. Mer- 
both, Sailplane and Glider, vol. 20, Feb. 
1952, pp. 6-7. 

Bazooka Rocket Motor Testing Now 
Fully Automatic, by H. E. Stein, Jron 
Age, vol. 169, Mar. 20, 1952, pp. 100-102. 

Hot-Rod Jets, unsigned, Aviation Week, 
vol. 56, Mar. 24, 1952, p. 14. 

Determination of Rocket Motor Heat 
Transfer Coefficients by the Transient 
Method, by S. Greenfield, J. Aero. Sci., 
vol. 18, Aug. 1951, pp. 512-518. 


Heat Transfer and Fluid 
Flow 


Concerning the Cooling of a Body in an 
Air Stream Through the Injection of Cold 
Air into the Boundary Layer, by K. 8S. 
Muller, translated from German by E. 
Knuth, Zentrale fiir Wissenschaftliches 
Berichtwesen Res. Report 2020, 29 pp. 

Generalized Theory of Convective Heat 
Transfer in a Free-Molecule Flow, by A. 
K. Oppenheim, Univ. Calif. Report HE- 
150-93, Mar. 1952, 11 pp. + 9 pp. tables. 

On the Convergence of Numerical 
Solutions of the Heat-Flow Equation, by 
F. B. Hildebrand, J. Math. Phys., vol. 30, 
Apr. 1952, pp. 35-41. 

Experimental Investigation of Heat 
Transfer Through Laminar and Turbulent 
Boundary Layers on a Cooled Flat Plate 
at a Mach Number of 2.4, by E. G. Slack, 
NACA TN 2686, Apr. 1952, 13 pp. 

Liquid Metals Are Good Heat-Transfer 
Agents, by T. Trocki, Gen. Elec. Rev., vol. 
55, May 1952, pp. 22-25. 

Heat and Momentum Transfer in 
Turbulent Flow Mercury, by 8. E. Isakoff, 
Atomic Energy Comm. Report AECU- 
1199, 14 pp. 

Experimental Verification of the Theory 
of Jet Formation by Charges with Lined 
Conical Cavities, by R. J. Eichelberger 
and E. M. Pugh, J. Appl. Phys., vol. 23, 
May 1952, pp. 537-542. 

Theory of Jet Formation by Charges 
with Lined Conical Cavities, by E. M. 
Pugh, R. J. Eichelberger, and N. Rostoker, 
J. Appl. Phys., vol. 23, May 1952, pp. 
532-536. 

Heat Transfer in Nitrogen and Hydro- 
gen Sweat-Cooled Tubes, by H. L. 
Wheeler, Jr., C.1.T. Jet Propulsion Lab. 
Progress Report 20-160, Jan. 4, 1952, 27 
pp. 

Internal-Film-Cooling Experiments with 
2- and 4-Inch Smooth-Surface Tubes and 
Gas Temperatures to 2000 F, by G. R. 
Kinney, NACA RM E52B20, April. 1952, 
27 pp. 
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Note on the Application of Near-Frozen 
etlahas Flow Criteria for One-Dimensional Non- 
viscous. Expansion Through a _ Laval 
a _ Nozzle, by S. S. Penner, J. Chem. Phys., 
aa vol. 20, Feb. 1952, p. 341. 
On the Uniqueness of Gas Flows, by J. 
LL. Ericksen, J. Math. Phys., vol. 30, Apr. 
1952, pp. 63-68. 

Investigation of Flow in Duct of Vary- 
fng Cross Section with Application to De- 
sign of Ducts for Free-Flight Ground- 
Launched Model Tests, by C. H. E. 
Warren, R. E. Dudley, and P. J. Herbert, 
Aero. Res. Counc. Lond. Curr. Pap. 60, 

1951, 37 pp., Brit. Inf. Serv., New York 


n° _ Theoretical Aerodynamic Properties for 
an Inclined Flat Plate in Slip Flow, by 
«$F. Mack, Univ. Calif. Report HE-150- 


98, Mar. 1952, 14 pp. + 5 pp. tables. 
Dynamics of Pulsative Flow Through 
Jrifices, by R. C. Baird, Instruments, vol. 
5, Apr. 1952, pp. 481-486. 
Dynamical Instability in Flow Systems, 
_ by T. Baron, article from Proc. Midwest 
Conf. Fluid Dynamics, Ist conf., May 
1950, J. W. Edwards, Ann Arbor, Mich., 
1951, pp. 216-225, $10. 

Diffusion and Flow of Gases and Vap- 
ours Through Micropores. IV. Flow of 
_ Capillary Condensate, by P. C. Carman, 
Proc. Roy. Soc. Lond. (A), vol. 211, Mar. 
20, 1952, pp. 526-535. 

The Flow of Gases Through Capillaries, 
by J. A. W. Huggill, Proc. Roy. Soc. Lond. 
- (A), vol. 212, Apr. 8, 1952, pp. 123-136. 
Evaporation from Drops. II, by W. E. 
~Ranz and W. R. Marshall, Jr., Chem. Engg. 
Prog., vol. 48, Apr. 1952, pp. 173-180. 

Note on the Flow of Vapor Between 
Liquid Surfaces, by M. S. Plesset, J. 
a Chem. Phys., vol. 20, May 1952, pp. 790- 
_ Some Aspects of Fluid Flow, Edward 
ees & Co., Lond., 1951, 50s. See 


especially: Problems in the Atomisation 
of Liquids, by H. L. Green, pp. 75-86, and 
a - Boundary Layers and Skin Friction in a 

- Compressible Fluid Flowing at High 
Speeds, by A. D. Young, pp. 87-113. 
The Emission of Radiation from Dia- 
tomic Gases. III. Numerical Emissivity 
Calculations from Carbon Monoxide for 
Low Optical Densities at 300°K and 
Atmospheric Pressure, by S. S. Penner, 
M. Ostrander, and H. S. Tsien, Guggen- 
heim Jet Propulsion Center Tech. Report 
1, 1951, 8 pp. 


_ The Thermodynamics of Combustion 
Gases: General Considerations, by S. R. 
_ Brinkley, Jr., and B. Lewis, Dept. Interior 
Bur. Mines RI 4806, Apr. 1952, 61 pp. 
Combustion Tunnel Laboratory, First 
Quarterly Report, 1952, by H. W. Em- 
7 +. mons, Harvard Univ. Combustion Tunnel 
Proj. Qtrly. Report 1, Oct. 1951—Apr. 1, 


1952. 
Spectroscopic Studies of Low Pressure 
Combustion Flames, by 8S. S. Penner, M. 
; Gilbert, and D. Weber, J. Chem. Phys., 
vol. 20, Mar. 1952, pp. 1-2. 
Flame Propagation Studies Using the 
Hydrogen-Bromine Reaction, by 8. D. 
_ Cooley and R. C. Anderson, Jndus. Engg. 


Combustion 


Chem., vol. 44, June 1952, pp. 1402-1406. 


Combustion in the Mixing Zone Be- 
tween Two Parallel Streams, by F. H. 
Wright and J. L. Becker, C.1.T. Jet Pro- 
pulsion Lab. Progress Report 3-25, June 9, 
1952, 16 pp. 

Evaluation of Flame Speed at Burner 
Flame Tip, by L. C. Lichty, Ind. Engg. 
Chem., vol. 44, June 1952, pp. 1395-1398. 

Burning Velocity Determinations, by 
P. J. Wheatley, Trans. Faraday Soc., vol. 
44, Apr. 1952, pp. 338-345. 

Uber den Einfluss starker Schallwellen 
auf Fortschreitende Gasflammen in 
Rohren, by L. Ehret, U. Neubert, and H. 
Hahnemann, ZAP, vol. 4, Apr. 1952, pp. 
126-134. 

Non-Isotropic Propagation of Combus- 
tion Waves, by G. H. Markstein, J. Chem. 
Phys., vol. 20, June 1952, pp. 1051-1052. 

High Temperature Propellant Reac- 
tions, by F. E. Osborne, Aircr. Engg., vol. 
24, May 1952, pp. 134-137. 

The Combustion of Drops in a Fuel 
Spray, by G. Godsave, Nat. Gas Turbine 
Establishment, England, Memo. M-95, 
Oct. 1950, 10 pp. 

Thermal Decomposition of Nitroethane 
and 1-Nitropropane, by T. L. Cottrell, 
T. E. Graham, and T. J. Reid, Trans. 
Faraday Soc., vol. 47, 1951, pp. 1089-1092. 

Gaseous Detonations. III. Dissocia- 
tion Energies of Nitrogen and Carbon 
Monoxide, by G. B. Kistiakowsky, H. T. 
Knight, and M. E. Malin, J. Chem. Phys., 
vol. 20, May 1952, pp. 876-883. 

Reaction of Methyl Radicals with 
Nitric Oxide and the Rate of Combination 
of Methyl Radicals, by R. W. Durham 
and E. W. R. Steacie, J. Chem. Phys., 
vol. 20, Apr. 1952, pp. 582-585. 

The Heats of Combustion of the Methyl 
Substituted Hydrazines and Some Obser- 
vations on the Burning of Volatile Liquids, 
by J. G. Aston, E. J. Rock, and S. Isserow, 
J. Am. Chem. Soc., vol. 74, May 20, 1952, 
pp. 2484-2486. 

Thermal Theory of Solid-Propellant 
Ignition by Hot Wires, by D. Altman and 
A. F. Grant, Jr., C.I.T. Jet Propulsion 
Lab. Ext. Pub. 140, May 23, 1952, 9 pp. 

Steady State Burning of a Semi-Infinite 
Solid or Liquid, by W. P. Reid, NAVORD- 
R-1935, Dec. 7, 1951. 

Gasification of Solid Fuels at Elevated 
Pressures, by W. Gumz, Jnd. Engg. Chem., 
vol. 44, May 1952, pp. 1071-1074. 

Rates of Explosive Decomposition of 
Explosives, by H. Henkin and R. McGill, 
Ind. Eng. Chem., vol. 44, June 1952, pp. 
1391-1394. 

Mode of Energy Release in Combustion 
of Carbon, by J. R. Arthur and J. A. 
Bleach, Ind. Eng. Chem., vol. 44, May 
1952, pp. 1028-1034. 

Reactions of Artificial Graphite, by E. 
A. Gulbransen and K. F. Andrew, Ind. 
Eng. Chem., vol. 44, May 1952, pp. 1034— 
1051. 

Le domaine d’inflammation du méthane 
aux basses pressions, by M. Van Pee and 
G. Fally, Bull. soc. chim. Belges, vol. 61, 
Jan.-Feb. 1952, pp. 64-98. 

The Stability of Butane Air Flames at 
Low Pressures, by K. Wohl and C. Peter- 


son, United Aircr. Corp. Meteor Report 
54, Feb. 1952. 

The Oxidation of Hexane in the Cool 
Flame Region, by H. C. Bailey and R. G. 
W. Norrish, Proc. Roy. Soc. Lond. (A), 
vol. 212, May 7, 1952, pp. 311-330. 

Polyhedral Flames with Hydrogen and 
Hydrocarbons, by H. Broida and W. Kane, 
J. Chem. Phys., vol. 20, June 1952, p. 1042. 

Induced Decomposition of Acetaldehyde 
by Radicals, by R. K. Brinton and D. H. 
Volman, J. Chem. Phys., vol. 20, June 
1952, pp. 1053-1054. 

The Reaction of Methyl] Radicals with 
Some Halogenated Methanes, by F. A. 
Raal and E. W. R. Steacie, J. Chem. Phys.. 
vol. 20, Apr. 1952, pp. 578-581. 

The Azomethane Induced Oxidation 0: 
Propane, by M. D. Scheer and H. A 
Taylor, J. Chem. Phys., vol. 20, Apr. 1952 
pp. 653-657. 


Fuels, Propellants, and 
Materials 


Fluoroscope Tests Rocket Propellants, 
unsigned, Tech. Data Digest, vol. 17, May 
1952, p. 13. 

How One Base Meets Problems ot 
Handling Rocket Fuels, unsigned, Tech. 
Data Digest, vol. 17, May 1952, pp. 3-5. 

Interstate Commerce Commission 
Regulations for Transportation of Ex- 
plosives and Other Dangerous Articles by 
Land and Water in Rail Freight Service 
and by Motor Vehicle (Highway) and 
Water Including Specifications for Ship- 
ping Containers, by H. A. Campbell, 
Interstate Commerce Commission 4, April 
8, 1952, 340 pp. 

The Borohydrides, by H. C. Brown, 
Chem. Eng. News, vol. 29, 1951, pp. 5231- 
5233. 

Reaction Mechanism of the Raschig 
Synthesis of Hydrazine, by E. Wiberg and 
M. Schmidt, Z. Naturforsch., vol. 6b, 1951, 
p. 336. 

Powdered Metals; Old Art, New Sci- 
ence, by 8S. Baker, Steel, vol. 130, Apr. 28, 
1952, pp. 56-57. 

Titanium; the New Light Metal, un- 
signed, Light Metals, vol. 15, Apr. 1952, 
pp. 120-121. 

Carbon and Graphite Materials and 
Parts, by P. O’Keefe, Materials and 
Methods, vol. 35, Apr. 1952, pp. 119-134. 

New Refractories—Carbides, Borides, 
Nitrides, Silicides, by G. R. Finlay, /. 
Electrochem. Soc., vol. 99, Mar. 1952, pp. 
58c-61c. 

Determination of Particle Size in the 
Sub-Micron Range, by Y. Schubert, 
Powder Metal Bull., vol. 6, Apr. 1952, pp. 
105-109. 

Experimental Refractory Bodies of 
High-Melting Nitrides, Carbides, and 
Uranium Dioxide, by P. Chiotti, J. Am. 
Ceram. Soc., vol. 35, May 1, 1952, pp. 123 
130. 

Stabilization of Zirconia with Calcia and 
Magnesia, by P. Duwez, J. Am. Ceram. 
Soc., vol. 35, May 1, 1952, pp. 107-113. 

Vapor Deposition of Metals on Ceramic 
Particles, by J. E. Cline and J. Wulff, J. 
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Electrochem. Soc., vol. 98, Dec. 1951, pp. 
385-387. 

Ceramic Coating Saves Scarce Metals 
and Lengthens Life of High Temperature 
Parts, by H. Penton, Western Metals, vol. 
10, Apr. 1952, pp. 37-39. 


Physical-Chemical Topics 


Nucleation Phenomena, Symposium 
containing 17 papers, Ind. Eng. Chem., vol. 
44, June 1952, pp. 1270-1336. 

New Thermodynamic Concepts Bring 
Theory and Practice into Agreement, by 
N. P. Bailey, Chem. Eng., vol. 59, Mar. 
1952, pp. 150-154, 336-338. 

Thermodynamics of Irreversible Proc- 
esses, by I. Prigogine, Appl. Mechanics 
Reviews, vol. 5, May 1952, pp. 193-195. 

Experimental and Theoretical Activities 
in the Gas Properties Field, 1948 to date, 
by I. F. Weeks, Univ. Wise. CM 703, Jan. 
1952, 47 pp. 

Equation of State of Gases at High 
Temperature, by H. M. Peek and Z. W. 
Salsburg, J. Chem. Phys., vol. 20, Apr. 
1952, p. 763. 

Second Virial Coefficient of a 6-9 Gas at 
High Temperatures, by L. F. Epstein and 
C. J. Hibbert, J. Chem. Phys., vol. 20, 
Apr. 1952, pp. 752-753. 

The Second Coefficient of Viscosity of 
Liquids and Gases, by 8. M. Karim and L, 
Rosenhead, Rev. mod. Phys., vol. 24, Apr. 
1952, pp. 108-116. 

Measurement of the Vibrational Re- 
laxation Effect in COz by Means of Shock 
Tube Interferograms, by E. F. Smiley, 
Ic. H. Winkler, and Z. I. Slawsky, J. Chem. 
Phys., vol. 20, May 1952, pp. 923-924. 

Equilibrium of Solution in the System 
Nitromethane—Water, by R. M. Corelli, 
Aerotecnica, vol. 30, 1950, pp. 32-37. 

Decomposition of Nitrogen Pentoxide in 
the Presence of Nitric Oxide. II. Details 
at Low Pressures, by H. 8. Johnston and 
R. L. Perrine, J. Am. Chem. Soc., vol. 73, 
1951, pp. 4782-4786. 

Four Mechanisms Involving Nitrogen 
Pentoxide, by H. 8. Johnston, J. Am. 
Chem. Soc., vol. 73, 1951, pp. 4542-4546. 

The Thermal Decomposition of Potas- 
sium Perchlorate and Perchlorate-Halo- 
genide Mixtures: a Study in the Pyrol- 
vsis of Solids, by A. Glasner and L. 
Weidenfeld, J. Am. Chem. Soc., vol. 74, 
May 20, 1952, pp. 2467-2472. 

Thermal Decomposition of Potassium 
Chlorate and Chlorate-Chloride Mixtures, 
by A. Glasner and L. Weidenfeld, J. Am. 
Chem. Soc., vol. 74, May 20, 1952, pp. 
2464-2467. 

Photoelectric Observation of the Rate of 
Dissociation of Nitrogen Tetroxide by a 
Shock Wave, by T. Carrington and N. 
Davidson, J. Chem. Phys., vol. 19, 1951, 
p. 1313. 

Measurement of the Heat of Dissocia- 
tion of Fluorine by the Effusion Method, 
by H. Wise, J. Chem. Phys., vol. 20, May 
1952, p. 927. 

The Electron Affinity of Fluorine, by 
R. B. Bernstein and M. Metlay, J. Chem. 
Phys., vol. 19, 1951, p. 1612. 

Thermodynamic Properties of the Flu- 
orine Atom and Molecule and of Hydrogen 
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Fluoride to 5000°K, by L. G. Cole, M. 
Farber, and G. W. Elverum, Jr., J. Chem. 
Phys., vol. 20, Apr. 1952, pp. 586-590. 

Thermodynamic Properties of Chlorine 
Trifluoride, by M. D. Scheer, J. Chem. 
Phys., vol. 20, May 1952, p. 924. 

Bibliography on the Thermodynamic 
Properties of the Fluorine Compounds of 
Elements of Atomic Number 1 to 17, by 
W. H. Evans and others, Dept. Com- 
merce, Nat. Bur. Stds., Jan. 21, 1952, 17 
pp. 

Thermal Properties of Fluorine Com- 
pounds: Heat Capacity, Entropy, Heat 
Content, and Free Energy Functions of 
Diatomic Fluorine in the Ideal Gaseous 
State, by L. Haar and C. W. Beckett, 
Dept. Commerce, Nat. Bur. Stds., Feb. 1, 
1952, 16 pp. 

Thermal Properties of Fluorine Com- 
pounds: Heat Capacity, Entropy, Heat 
Content, and Free Energy Functions of 
Carbon Difluoride, Carbon Tri-Fluoride, 
Carbon Tetrafluoride and Difluoro-Acet- 
ylene in the Ideal Gaseous State, by R. M. 
Potocki and D. E. Mann, Dept. Com- 
merce, Nat. Bur. Stds., Feb. 15, 1952, 4 pp. 

Recent Determinations of the Vapor 
Pressure of Graphite, by L. Brewer, J. 
Chem. Phys., vol. 20, Apr. 1952, pp. 758- 
759. 


Instrumentation and 
Experimental 
Techniques 


Electromagnetic Levitation of Solid and 
Molten Metals, by E. C. Okress and D. M. 
Wroughton, J. Appl. Phys., vol. 23, May 
1952, pp. 545-552. 

A Linear Radio-Frequency Mass Spec- 
trometer, by P. A. Redhead, Canad. J. 
Phys., vol. 30, Jan. 1952, pp. 1-9. 

Résumé of Thermocouple Checking 
Procedures, by H. C. Quigley, /nstruments, 
vol. 25, May 1952, pp. 616-622. 

Resistance-Thermometer Measurements 
in a Low-Pressure Flame, by M. Gilbert 
and J. H. Lobdell, C.I.T. Jet Propulsion 
Lab. Ext. Pub. 131, Feb. 11, 1952, 23 pp. 

The Use of Doubly Refracting Solutions 
in the Investigation of Fluid Flow Phe- 
nomena, by B. Rosenberg, USN David 


Taylor Model Basin Report 617, Mar. | 


1952, 39 pp. 


Determination of Air Velocity by Ion 
Transit-Time Measurements, by W. C. 


Cooley and H. G. Stever, Rev. Sci. In- 
strum., vol. 23, Apr. 1952, pp. 151-154. 

An Investigation of Electromagnetic 
Flowmeters, by H. G. Elrod and R. R. 
Fouse, Trans. ASME, vol. 74, May 1952, 
pp. 589-594. 

Valve Control of Liquid Flow, by R. 
Milham and A. R. Catheron, /nstruments, 
vol. 25, May 1952, pp. 596-598. 

A Modified Budenberg Pressure Balance 
by J. F. Pearse, J. Sci. Instrum., vol. 29, 
Jan. 1952, pp. 18-19. 

A Miniature Electrical Pressure Gage 
Utilizing a Stretched Flat Diaphragm, by 
J. Patterson, NACA TN 2659, Apr. 1952, 
47 pp. 

Accelerometers for Determining Air- 
craft Flight Loads, by J. Taylor, Engineer- 


ing, vol. 173, Apr. 11, 1952, pp. 473-475. 

An Instrument to Measure and Com- 
pute the Mean and Standard Deviation of 
the Oscillations of a Recording Arm, by B. 
Hegley, J. Sct. Instrum., vol. 29, Mar. 
1952, pp. 76-81. 

Noise, Explaining Its Nature and Ori- 
gin, unsigned, Wireless World, vol. 58, 
May 1952, pp. 199-202. 

On Sound Generated Aerodynamically. 
I. General Theory, by M. J. Lighthill, 
Proc. Roy. Soc. Lond. (A), vol. 211, Mar. 
20, 1952, pp. 564-587. 

The Noise Field of a Turbo-Jet Engine, 
by H. E. von Gierke, J. Acoust. Soe. 
Amer., vol. 24, Mar. 1952, pp. 169-174. 

A Survey of the Aircraft-Noise Problem 
with Special Reference to Its Physical 
Aspects, by H. H. Hubbard, NACA TN 
2701, May 1952, 41 pp. 

Design Characteristics for Noise Con- 
trol of Jet Engine Cells, by H. C. Hardy, 
J. Acoust. Soc. Amer., vol. 24, Mar. 1952, 
pp. 185-190. 

The Design of Mobile Electronic Lab- 
oratories, by P. H. Parkin, J. Sci. In- 
strum., vol. 29, Jan. 1952, pp. 8-11. 

Tape Recording for Telemetry and 
Data Analysis, by K. B. Boothe, Tele- 
Tech., vol. XI, May 1952, pp. 44-46 and 
116. 

Guided Missile Test Center Telemeter- 
ing System, by J. B. Wynn, Electronics, 
vol. 25, May 1952, pp. 106-111. 

A High-Speed Crystal Clutch, unsigned, 
Engineer, vol. 173, Apr. 25, 1952, p. 579. 

Response Time of Magnetic Fluid 
Clutches, unsigned, Engineer, vol. 173, 
Apr. 25, 1952, p. 579. 

Missile Timers, unsigned, 
Week, vol. 56, May 19, 1952, p. 73. 

Calculating Speed of Snap-Action 
Mechanisms, by P. H. Winter, .Wachine 
Design, vol. 24, May 1952, pp. 162-164. 

A 24-Way High-Speed Rotary Switch 
for Use in Static and Airborne Strain 
Gauge Measurements, by D. Peirson, 
Aero. Res. Counce. Lond. Rep. Mem. 2232, 
1951, 9 pp., 1 diagram, 3s. 

Modern Computing Machines, by G. T. 
Hunter, J. Franklin Inst., vol. 253, June 
1952, pp. 567-583. 


Aviation 


Terrestrial Flight, 
Ballisties, and Vehicle 
Design 


Parachute Design, by T. F. Johns, 
Aero. Res. Counc. Lond. Rep. Mem. 2402, 
Dec. 1946 (published 1951), 23 pp., 11 figs. 

Interferometric Analysis of Airflow 
About Projectiles in Free Flight, by F. D. 
Bennett, W. C. Carter, and V. E. Bergdolt, 
J. Appl. Phys., vol. 23, Apr. 1952, pp. 453 
469. 

Dynamie Stability of a Missile in Roll- 
ing Flight, by R. E. Bolz, J. Aero. Sei., 
vol. 19, June 1952, pp. 395-403. 

Model Simulates Free Flight of Rotat- 
ing Missiles, unsigned, 7’ech. Data Digest, 
vol. 17, May 1952, p. 9. 

Notes on Dynamics for Aerodynami- 
cists, by J. D. Schetzer, Douglas Airer. 
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Co., Inc., Report SM-14077, Nov. 19, 
1951, 203 pp. 

Rocket Flight, by S. F. Singer, Canad. 
Aviation, vol. 25, May 1952, pp. 32-33. 

On the Calculation of Groups of Trajec- 
tories of Remote Controlled Rockets (in 
German), by E. Roth-Desmeules, ZAM P, 
vol. 2, Nov. 1951, pp. 487-489. 

Remarkable Points on a Trajectory. IV. 
Curvilinear Trajectories. Discontinuities, 
Extrema, Inflections (in French), by M. 
Garnier, Mém. Artill. fr., vol. 25, 1951, pp. 
693-741. 

Ballistics of the Upper Atmosphere; 
Atmospheric Density Profile and Varia- 
tions from the Study of Meteor Trajec- 
tories, by L. G. Jacchia, MIT Center of 
Analysis Report 2, 1948, 30 pp. 

Flight in the Aeropause, by R. Haw- 
thorne, Aviation Age, vol. 16, Dec. 1951, 
pp. 29-31. 

Radio Direction-Finding and Naviga- 
tional Aids, unsigned, Nature, vol. 169, 
May 24, 1952, p. 875. 

Navigation in War, unsigned, Flight, 
vol. 61, Feb. 29, 1952, p. 246. 

Deadly Sparrow Missile Revealed, un- 
signed, Amer. Aviation, vol. 15, Mar. 31, 
1952, p. 61. 

General Electric’s Winged Messenger; 
Hermes Al, unsigned, Aviation Week, vol. 
56, May 12, 1952, pp. 30 and 35-36. 

The Armstrong-Siddeley “‘Snarler’” Roc- 
ket Motor, unsigned, J. Brit. Interplane- 
tary Soc., vol. XI, Mar. 1952, pp. 73-75. 


V-2 Rocket Tests, unsigned, Mech. 


Eng., vol. 74, Apr. 1952, p. 329. 
Red Rockets, by D. A. Anderton, 
Aviation Week, vol. 56, Jan. 14, 1952, pp. 


37, 38, 41. oe 
Space Flight 


How Far Are We From Space Flight?, 
by F. C. Durant, Aviation Week, vol. 56 
May 26, 1952, pp. 25-35. 

Interplanetary Communications and 
Navigation, by J. P. Elliott, J. Space 
Flight, vol. 4, May 1952, pp. 1-2. 

Interplanetary Navigation, by 8S. Her- 
rick, J. Inst. Navigation, vol. 3, Sept.- 
Dec. 1951, pp. 57-59. 

Verstandigungsméglichkeiten mit Raum- 
schiffen, by G. von Schrutka-Rech- 


tenstamm, Forschungsrethe der Nord- 
westdeutchen Gesellschaft fiir Weltraum- 
forschung, no. 2, Mar. 1951, Sonthofen, 


Germany. 

Eine Methode interplanetarischer Orts- 
bestimmung, by C. Oesterwinter, Forsch- 
ungereithe der Nordwestdeutschen Gesell- 
schaft fiir Weltraumforschung, no. 2, Mar. 
1951, Sonthofen, Germany. 

Meteoric Collision Factor in Space Ship 
Design, by F. L. Whipple, Aviation Age, 
vol. 16, Dec. 1951, pp. 25-26. 

The Cost of Interplanetary Cargo 
Transportation. II, by N. J. Bowman, 
J. Space Flight, vol. 4, Apr. 1952, pp. 1-9. 

A Method of Landing a Space Ship 
Under Adverse Climatic Conditions, by 


M. Conley, 
1952, p. 6. 

Artificial Satellites, by J. Humphries, 
Aeronautics, vol. 25, Apr. 1952, pp. 62, 65- 
66, and 69. 

Generalized Interplanetary Orbits, by 
H. Preston-Thomas, J. Brit. Interplane- 
tary Soc., vol. XI, Mar. 1952, pp. 76-85. 

Ecological Aspects of Planetary Atmos- 
pheres with Special Reference to Mars, by 
H. Strughold, J. Aviation Medicine, vol. 
23, Apr. 1952, pp. 130-140. 

Medical Research in Development of 
Manned Rocket Flight, by H. Haber, 
Tech. Data Digest, vol. 17, Feb. 1952, pp. 
12-13. 

Human Experiments in Subgravity and 
Prolonged Acceleration, by E. R. Bal- 
linger, Tech. Data Digest, vol. 17, May 
1952, pp. 10-12. 

Man Without Gravity; the Physiologi- 
cal and Psychological Problems of Space- 
Flight, by L. N. Thompson, Flight, vo! 
61, Mar. 14, 1952, pp. 298-300. 


J. Space Flight, vol. 4, May 


Astrophysics and 
Aerophysics 

The Origin of the Planets, by H. Brown, 
Chem. Eng. News, vol. 30, Apr. 21, 1952, 
pp. 1622-1626. 

The Colloquium on Auroral Physics, by 
N. C. Gerson, J. Franklin Inst., vol. 253, 
Apr. 1952, pp. 331-338. 

The Role of Turbulence in the Evolu- 
tion of the Universe, by G. Gamow, Phys. 
Rev., vol. 86, Apr. 15, 1952, p. 251. 


“Here is is a better lining material 


Shown here is a group of liner and throat shapes made sel atl Carb dum C 


Silicon Carbide linings by 


“GaRBORUNDUM 


Trade Mark 


‘for uncooled rocket motors 


Exceptional resistance to the high 
temperature and severe erosive 
conditions developed in uncooled 
rocket motors is offered by spe- 
cial silicon carbide liners and 
throats. In comparative tests they 
have outlasted metals by several 
times. 

Information is available upon 
request to Dept. RS-52, Refractor- 
ies Div., The Carborundum Co., 
Perth Amboy,N.J. 
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Structure of Comets and the Formation 
of Tails, by R. A. Littleton, Nature, vol. 
169, Apr. 12, 1952, pp. 615-616. 

On the Possibility of Obtaining Radar 
Echoes from the Sun and Planets, by F. J. 
Kerr, Proc. Instn. Radio Engrs., vol. 40, 
June 1952, pp. 660-666. 

Solar Radio Noise Bursts, unsigned, 
Tech. News Bull., vol. 26, May 1952, pp. 
65-67. 

Radio Astronomy, by J. A. Ratcliffe, 
Nature, vol. 169, Mar. 1, 1952, pp. 348- 
350. 

A New Radio Interferometer and Its 
Application to the Observation of Weak 
Radio Stars, by M. Ryle, Proc. Roy. Soc. 
Lond. (A), vol. 211, Mar. 6, 1952, pp. 351- 
375. 

A Radio-Astronomical Investigation of 
Winds in the Upper Atmosphere, by A. 
Maxwell and C. G. Little, Nature, vol. 
169, May 3, 1952, pp. 746-747. 

Average Radio-Ray Refraction in the 
Lower Atmosphere, by M. Schulkin, 
Proc. Instn. Radio Engrs., vol. 40, May 
1952, pp. 554-561. 

Atmospheric Density Profile and Gradi- 
ents from Early Parts of Photographic 
Meteor Trails, by L. G. Jacchia, MIT 
Center of Analysis Report 4, 1949, 12 pp. 

Photographic Meteor Phenomena and 
Theory; Meteor Photometry, Funda- 
mental Equations and Constants, Dura- 
tions and Flares, by L. G. Jacchia, MIT 
Center of Analysis Report 3, 1949, 36 pp. 

Research in the Upper Atmosphere 
with Sounding Rockets and Earth Satel- 
lite Vehicles, by S. F. Singer, J. Brit. 
Interplanetary Soc., vol. XI, Mar. 1952, pp. 
61-73. 

The World’s High Altitude Labora- 
tories, by S. A. Korff, Physics Today, vol. 
5, May 1952, pp. 28-30. 

On the Efficiency of the Engine Driving 
the Atmospheric Circulation, by O. R. 
Wulf, J. Meteor., vol. 9, Apr. 1952, pp. 79- 


82. 
Other Classifications 


Construction of an Intermediate Reac- 
tor for Submarines, unsigned, Engineer, 
vol. 173, Apr. 25, 1952, p. 579. 

Errors in Rocket Development, 1, by 
H. Oberth, Rocketscience, vol. 6, Mar. 
1952, pp. 2-7. 

Problems of Astronautical Research, by 
E. Sanger, J. Brit. Interplanetary Soc., vol. 
XI, Mar. 1952, pp. 57-60. 

Let Punched Cards Solve That Prob- 
lem, by R. Habermann, Gen. Elec. Rev., 
vol. 55, May 1952, pp. 9-12. 

A Discussion on Friction, a group of 30 
notes and papers on friction organized 
under 3 headings; I. Friction of Metals, 
II. Friction of Non-Metals, III. Boun- 
dary and Extreme Pressure Lubrication, 
by F. P. Bowden, et al., Proc. Roy. Soc. 
Lond. (A), vol. 212, May 22, 1952, pp. 
439-520. 

Note sur l’autopropulsion avee source 
d’énergie séparée, by M. P. Blanc, Mém. 
de lartillerie francaise, vol. 25, no. 95, 
1951, pp. 103-116. 

The Role of Physics in Aeronautical De- 
velopment, by H. L. Dryden, Physics To- 
day, vol. 5, May 1952, pp. 14-24. 
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MODEL 1809 
BRIDGE CONTROL UNIT 


FOR VIBRATION AND STRESS ANALYSIS 


Designed as a companion unit to Century's famous Model 409 Oscillograph, the Model 
1809 Bridge Control Unit is the latest addition to Century’s line of industry-standard 


vibration and stress analyzing equipment. 


Packaged in a small, compact space, the unit 


contains all of the facilities necessary for use with 12 channels of resistance strain gages 
or bridge-type transducers. Where used with the Model 409 Oscillograph, it is necessary 
only to connect strain gages and power source to have a complete stress-strain measuring 


and recording system, small and rugged enough to be placed in an aircraft wing tip or — 


guided missile warhead. 


FEATURES: 

Size: 4%” x 7” x 11”. 
Weight: 10% pounds. 
Aluminum case. 


Up to 12 channels. 


Write for Bulletin CGC—307 


MODEL 
409 
OSCILLOGRAPH 


The Century Model 409 Oscillograph has been designed for recording data where space and 
weight requirements are limited. The Oscillograph has been tested to record faithfully while — 


subjected to accelerations up to 20 G's. 


FEATURES: 
Size: 5” x 615%," x 11°”, 
Weight: 13 pounds. 
Cast aluminum case. 


Paper speeds variable 2” to 6” _' 


and 2” to 24” per second. 


Write for Bulletin CGC—303 


For any resistance strain gage or bridge-type 
transducer. 
May be used with direct indicating instrument. 
Power: Control unit, 22-28 Volt D.C. 
Strain gage, 6-28 Volt D.C. 


Detachable daylight loading magazine with — 

a capacity of 356” x 100’ paper. 
2 to 14 individual channels. 
Trace identification. 


Trace viewing. 


GEOPHYSICAL CORPORATION 


TULSA, OKLAHOMA 


1505 Race St 
Philadeiphia, Pa 


3406 W. Washington Blvd. 
los Angeles 18, Calif 


238 Lafayette St 
Dayton 2, Ohio 


EXPORT OFFICE 
149 Broadway, N. City 


309 Browder St. 


Dallas, Texas 
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GROVE 
SELF-ACTUATED 


CONTROLS AND ‘RESULATORS 


High pressure liquids slid gases 


Pressure Reducing 
Back Pressure = 


Relief and By-Pass 
Manual Valves 


Special light weight units wh 
for air-borne service. 


_ Special materials for corrosive 
ae or low-temperature applications. 


GROVE CONTROLS, INC. 


6529 Hollis Street 
Emeryvil le 
ACCURACY 
SAFETY: 


California 


RINGS 
MANUFACTURED OF POLYETHYLENE, TEFLON, 
AND KELL-F TO A.N. SPECIFICATIONS 
HOSES 
FLEXIBLE HIGH PRESSURE — BRAIDED STAINLESS 
STEEL POLYETHYLENE LINED WITH A. N. END 
CONNECTIONS 

3 POLYETHYLENE SHEET, RODS, AND TUBES 

E ; LARGE STOCK AVAILABLE FOR GASKETING AND 

ae CHEMICAL CARRIAGE 

VALVES 
STAINLESS STEEL CHECK VALVES—HIGH PRESSURE 
— POLYETHYLENE OR TEFLON POPPET 

4 

ef SPECIALTIES COMPANY 


3348 East 14th Street, Los Angeles, California 


accelerometers A 

NGULAR accelerometers 
are made by Statham Lab- 
oratories for measurement in 
ranges as low as + 1.5 radians 
per second per second. 


The design permits close 
balance against linear 
acceleration effects with a 
high degree of mechanical 
shock resistance and leads to 
a damping characteristic 
relatively insensitive 

to temperature. 

The transducer element, an 
unbonded strain gage bridge, 
provides an electrical output 
proportional to applied 
angular acceleration for 
recording or telemetering in 
conventional a.c. or 

d.c. circuits. 


LABORATORIES: 
los Angeles’ 


Get the Facts about 
your Future 


Discover the 
greater opportuni- 
ties offered engi- 
neers by the 
greatest diversity of 
projects of any air- 
craft company in the 
East! Write today 
for fact-packed 
brochure. 


SEND FOR ENGINEERING BROCHURE 


THE GLENN L. MARTIN COMPANY .Y 
Personnel Dept. « Section A « Baltimore 3, Md. .Y 
Please send me your brochure describing engineering opportuni- X 
ties at Martin. >) 
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Rides a Rocket ; | 


The Naval Research Laboratory’s 
Viking rocket research at White 
Sands Proving Grounds, N. M., 
hunts facts, figures and formulas 
in the upper atmosphere. 


URTLING far into the blue, Naval Research Labor- 
H atory rockets ask questions of the earth’s upper 
atmosphere . . . flash back the answers needed to guide 
the designers of tomorrow’s piloted and pilotless super- 
altitude systems for peace or war. What are the pres- 
sures and temperatures of the earth’s atmospheric 
layers . . . the high-altitude changes in the earth’s mag- 
netic field affecting navigational instruments . . . the 
alterations in radio waves caused by the ionosphere . . . 
the effects of sun spots on communications equipment 
out beyond the filtering effects of the earth’s heavy 
atmosphere? 


Martin Viking rockets play a major role in this high- 
altitude flight research program. Last summer, the 
Viking cracked the world’s altitude record for single- 
stage rockets . . . nosing 136 miles into the heavens at 
a top speed of 4100 m.p.h. Now, an even more powerful 
Viking is being readied for launching. The Martin 
Company is proud to be a partner with the Naval Re- 
search Laboratory in these vital activities . . . helping 
to prove that America’s most valuable secret weapon is 
its scientific leadership! Tot GLenn L. Martin Com- 
PANY, Baltimore 3, Md. 


d 


AIRCRAFT 


Builders of Dependable —~s# Aircraft Since 1909 

Developers and Manufacturers of: Navy P5M-1 Marlin seaplanes « Air 
Force B-57A Canberra night intruder bombers « Air Force B-61 Matador 
pilotiess bombers * Navy P4M-1 Mercator patrol planes *« Navy KDM-1 
Plover target drones * Navy Viking high-altitude research rockets « Air Force 
XB-51 developmental tactical bomber Martin airliners Guided missiles 
Electronic fire control & radar systems « Leaders In Building Air Power to 
Guard the Peace, Air Transport to Serve It. 


SEPTEMBER-OCTOBER 1952 


‘ ‘ 
AG 
n on : Za MB 
ABA || | ZZ 
Zz g BB : g g 
| Ih: 
| 
297 


Have You Inv 


as an 
Oxidizer in Rocket Propulsion ? 


4. If you are designing rocket motors using liquid propellants, you should 
consider Nitrogen Tetroxide as oxidant. 


Here are some of the properties of Nitrogen Tetroxide which have led 


bets to its use as oxidizer in rocket research: ~ 
HIGH SPECIFIC IMPULSE: 
N2O, exceeds many other well-known oxidizers in pounds of 


wat 2 He thrust developed per pound of fuel consumed per second. | 
EASY TO HANDLE: 
= N Rey may be shipped, piped and stored in ordinary carbon steel 
equipment. It possesses high chemical stability, high density, 
low freezing point, and a reasonably low vapor pressure. me ; 


ALLIED & DYE CORPORATION 
. 


40 ies STREET, NEW YORK 6, N. Y. 7 


Technical service and development on Nitrogen Tetroxide—formerly handled by the Product Development Department, 
Solvay Process Division—are now handled by Nitrogen Division, Allied Chemical & Dye Corporation. 
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INDUS TRIAL 


SOUND CONTROL 


INDUSTRIAL 


= 
| 
> 
eg 
N,0, is available at low cost in 125-pound I.C.C. approved steel cylinders. 
| ‘ 


‘omponents for 
automatic 


Pressure Transmitters 
Precision Potentiometers 
Accelerometers 
Aeroheads 


G. M. GIANNINI & CO., INC., Pasadena 1, California 
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cONTROL @ INDUSTRIAL SOUND O 
> 
= 
” 
- wocket 
silenced by 
=) 
a Zz 
> mufflers 
* @ In the laboratory, in the test cell, on the wn 
airfield, hundreds of INDUSTRIAL SOUND 
CONTROL installations are subduing the 
noise, heat and gas velocities generated dur- 
_ ing testing of the big jets. - 
a Whatever your noise problem, ISC’s skilled 7 > 
engineering, design, and installation “know 
1°) how” — gained through years of practical 
Vv experience — can help you, and is at your 7) 
instant service. O 
a We welcome the challenge of the unusual 
z problem. For full information c 
WRITE ISC TODAY z 
Foreign Licensees: 
Cementation (Muffelite) Limited, London 
Les Travaux Souterrains, Paris 
< 
= Industrial Sound Control, Inc. 
- Industrial ound Controi,inc. - 
> 45 Granby Street, Hartford, Conn. i 
a 2119 So. Sepulveda Bivd., Los Angeles, Calif. ° 
@®@inNDUSTRIAL SOUND CONTROL 


From Precision Instrument 
Screws... 


WALTHAM makes all 
sizes and types of 
precision instrument 
screws for industrial 
and scientific use, in 
sizes from #0000-200 
pitch up. Screws #0-80 
and #1-72 are in stock. 
We can furnish cut 
thread ultra-pre- 
cision roll thread 
products, as are now 
specified for many ap- 
plications. Parts are 
produced exactly to 
your specifications 
under rigid quality con- 
trol and inspection. 
No quantity is too large 
or too small. 

Send us your blueprints 
for free estimates that are 
right and reasonable. 


Special screws, screw 
machine parts — of all 
materials — particular- 
ly stainless steels, in all 
finishes. Includes cen- 
terless grinding. We 
specialize in manufac- 
ture of unusual parts. 


Remember, an instrument is 
only as good as its fasteners. 
Waltham Fasteners make 
good instruments better. 


to Accurate 
Metal Parts... 


SKILLED WORKERS, LATEST 
EQUIPMENT, MODERN PLANT 
are your guarantee of Quality. 


Depend on WALTHAM for 
Precision * Accuracy * Uniformity 


WALTHAM 


COMPANY 


SCREW 
76 Rumford Ave., Waltham 54, Mass. 


SCREW MACHINE PRODUCTS SINCE 1893 
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valuable addition 
to your technical library 


COSTS YOU NOTHING—CAN MEAN 
REAL SAVINGS IN TIME AND MONEY 
You'll find this 20-page compilation of technical data 
and manufacturing criteria a veritable gold-mine of 
hard-to-get electrical insulation information. Com- 
plete in content, it not only encompasses the wide 
range of MYCALEX Insulation in all its various 
grades and characteristics, but includes comparative 
data on other important dielectric materials as well. 
Write today. Your copy will be forwarded promptly. 


FOR QUICK REFERENCE — CONSULT 
THE 1950 IRE YEARBOOK 
For your added convenience the entire MYCALEX 
20-page catalog appeared as a section in the 1950 
IRE Yearbook. It’s the seventh catalog in the special 
manufacturers catalog section at the back of the 
1950 issue. 


SINCE 1919 


of glass-bonded mica 


ELECTRICAL 
INSULATION 


with special section on 
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MINIATURE 
LOW-LOSS 
TUBE SOCKETS 


TABLE OF CONTENTS 


THE STORY OF MYCALEX 
Its origin and history 
TYPES OF MYCALEX 
Injection molded-Compression molded 
CHARACTERISTICS OF MYCALEX 
Comparison to other insulators 
PROPERTIES OF DIELECTRICS 
A handy comparative chart for the engineer 
MYCALEX 410: Injection Molded 
Its qualities and engineering specifications 
MYCALEX 410X 
Engineering specifications as compared 
to Mycalex 410 
DESIGN OF MOLDED INSULATORS 
For the engineer 
MYCALEX TUBE SOCKETS 
9-pin, 7-pin and subminiatures 
MYCALEX 400: Compression Molded 
Applications and uses, engineering speci- 
fication and qualities, machining and 
recommended fabricating methods 
DESIGN OF. MACHINED INSULATORS 
Design criteria for insulators 
MYCALEX K: Capacitor Dielectrics 
Uses and Engineering data 
MYCALEX ADVISORY SERVICE 
Customer service available 


MYCALEX CORPORATION OF AMERICA 


Owners of ‘MYCALEX’ Patents and Trade-Marks 


Executive Offices: 30 ROCKEFELLER PLAZA, NEW YORK 20 — Piant & General Offices: CLIFTON, N.J. 
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